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"Regularization  and  Approximation  of  a  Class  of  Evolution  Problems  in  Applied  Mathematics" 

R.  Showalter  and  G.  F.  C#ey 


The  major  effort  of  this  project  has  been  the  development  of  the  foundations  for 
regularization  techniques  related  to  conservation  equations  and  some  new  possibly  far-reaching 
contributions  to  this  ^ea.  The  approach  that  has  been  taken  is  a  departure  from  the  usual  artificial 
viscosity  type  of  strategies  which  are  produced  on  a  somewhat  adhoc  basis.  The  basic  strategy  is 
to  regularize  locally  by  a  micro- structured  parabolic  system.  A  mathematical  analysis  of  the 
regularized  equations  has  been  developed  to  support  our  approach.  Supporting  approximate 
analysis  and  numerical  experiments  have  been  made. 

The  developnieht  and  the  mathematical  foundations  of  these  microstructure  models  have 
been  primary  achievements  of  the  project.  The  relevant  nonlinear  systems  of  panial  differential 
equations  have  also  been  shown  to  provide  good  rnodels  of  diffusion  or  convection  of  fluid  or  gas 
through  a  heterogeneous  porous  medium.  Examples  include  flow  in  fissured  media,  problems 
with,  adsorption,  heat  diffusion  with  freezing-melting,  and  models  for  semiconductors.  We  have 
established  that  the  problems  are  well-posed  and  developed  the  theory  of  the  regularity  and 
dependence  of  the  solutions  on  data.  Such  information  will  aid  approximation  theory  and  the 
design  of  algorithms  to  numerically  simulate  solutions  to  these  types  of  problems. 

The  major  step  came  with  the  rather  complete  development  in  [1]  of  the  linear  case  together 
with  appropriate  convergence  and  approximation-results.  These  results  were  extended  in  part  to 
fully  nonlinear  versions  in  [2j.  Classical  versions  were  given  in  [8].  This  work  was  summarized 
in  the  review  article  [3]  and  has  been  received  with  much  interest.  Moreover,  the  applications  to 
stationary  problems  [10]  and  to  porous  rnedia  [1 1]  are  underway. 


We  have  made  related  numerical  studies  using  finite  elements  with  the  regularizing  strategy 
and  the  results  are  promising.  We  have  also  been  developing  some  related  ideas  which  are  based 
on  superconvergence  concepts  in  the  approximate  methods  used  to  solve  both  boundary  and 
evolution  PDE's  [4,5].  This  is  a  very  topical  research  area  at  present  as  far  as  post  processing 
computed  solutions  is  concerned  [13,14].  Our  approach  is  different  in  that  we  are  using  the  post¬ 
processing  strategy  to  develop  improved^  models  and  to  develop  alternative  regularization 
strategies.  This  procedure  also  is  appropriate  for  discrete  homogenization*  at  a  macro-structure 
level  and  can  be  combined  with  statistical  averaging  at  the  micro-structure  level  as  a  regularization 
strategy. 

Numerical  studies  with  regularization  techniques  have  been  applied  to  flow  calculations 
[6,12].  This  work  is  being  extended  currently  to  least-squares  finite  element  analysis  including 
local  regularization.  Our  previous  work  with  a  least-squares  finite  element  formulation  and 
parabolic  regularization.f?]  (motivated  by  [9])  confirmed  that  this  type  of  regularization  procedure 
is  applicable  and  the  numerical  results  are  positive^  Some  of  the  challenging  aspects  of  the 
formulation  arise  in  combining  the  micfbstnicture  regularization  within  the  framework  of  a 
macrostructure  Galerkin  finite  element  analysis.  A  formulation  for  this  embedding  has  been 
developed. 
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1.  Introduction 

We  shall  study  the  Gaudiy-Dirichlet  prdblem  for  degenerate  parabolic  systems 
of  ' the  form 

!(.l.l-.a)  •— a(,u).^  V  •  A(a:,  Vu)-+  /  B(x,s,yj,£^)  •  i/ds  9.  •/ , 
ot  Jr. 

(i.i.b)  y-sn. , 

(Ivl.c)  B{x,  y,V-yU)  -  u  +  fi{U {x,  y,  t)  -  u(x,  <))*  9  0  ,  y  6  Pi  . 

?Here  ft  is  a  domain  in  R"  and  for  each  yadue  of  t^  macrq-yariable  x  €  ft  is:specified 
a  domain  ftx  with- boundary  P^;  fpr  the.  micro-variable  y-  E  ftx.  Each  of  a,  6, /x  is 
a  maximal  monotone  graph.  These  graphs  are  not  nece|sarily  strictly  increasing; 
they  may  be  piecewise  constant  or.multi-valued.  The  elliptic  operators  in  (1.1. a)  and 
;(l.l.b)-are  nonlineau:  in  the  gradient  of  degree  p-^  I  >  0  and  9  —  1  >  0,  respectively, 
with  ^  -f  ^  >  ^,<so  some  .specific  degeneracy  is  also  perniitted  here.  Certain  first 
order  spatial  derivatives  can  be  added  to  (l;l.a)  and  (l.l.b)  with  no  difficulty, 
and. corresponding  problems  with  constraints,  t.e.,  variational  inequalities,  can  be 
treated  similarly.  A  particular  example  important  for  applications  is  the  linear 

This  work  was  supported  by  grants  from  tlic  National  Science  Foundation  and  tlic  Office  of  .Naval 
Research. 
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constraint 

(l.l.c')  U{x,y,t)  =  u(^x,t)  ,  y  e  ,  x  € 

which  then  replaces  (l.l.c).  The  system  (1.1)'' with  fi{s)  =  is  called  a 

regularized  micro-structure  models  and  (i.l.a),  (^l.l.b),  (l.l.c')  is  the  corresponding 
matched  micro-structure  model  in  which  (formally)  e  — >•  0.  An  example  of  such  a 
system  as  a  model-for  the  flow  of  a  fluid  (liquid  or  gas)  through  a  fractured  medium 
will  be  given  below.  In  suA  a  context,  (l.l.a)  prescribes  the  flow  on  the  global  sc^e 
of  the  fissure  system  and.(l.l.b)  gives  the  flow  the  microscale  of  the  individual 
cell  at  a  specific  point  x  in  the  fissure  system.  The  transfer  of  fluid  between  the 
cells  and  surrounding  medium  is  prescribed  by  (-l.l.c)  or  (lil.c').  A  major  objective 
is  to  accurately  model  this  fluid  exchange  between  the  cells  and  fissures. 

The  plan  of  this  paper  is  as  follows.  In  Section  2  we  shEill  give  the  precise 
description  and  resolution  of  the  stationary  problem  in  a  variational  formulation  by 
monotone  opeirators  from  Banach  spaces  to  their  duals.  In  order  to  adSieve  this  we 
describe  first  the  relevantiSdbolev  spaces,  the  continuous  direct  sums  of  these  spaces, 
and  the  distributed  trace  and  constant  functionals  which  occur  in  the  system.  The 
operators  are  monotone  functions  or  multi-valued  subgradients*  and  serve  as  models 
for  nonlinear  elliptic  equsdions  in  divergence  form.  We  develop  an  abstract  Green’s 
theorem  to  describe  the  resolution  of  the  variationeil  forrh^^as  the  sum  of  a  partial 
differential  equation  and  a  complementary  boundary  operator.  Then  sufficient  con¬ 
ditions  of  coercivity  typetare  given  to  assert  the  existence- of  generalized  solutions 
of  the  vzuriation^  equations.  In  Section  3  we  describe  the  restriction  of  our  system 
to  appropriate  products  ol  spaces.  The  Hilbert  space  case,  r  =  2,  serves  not 
oiily  as  a  convenient  starting  point  but  leads  to  the  generalized  accretive  estimates 
we  shall  need  for  the  singular  case  of  (1.1)  in  which  a  or  b  is  not  only  honlineeir  but 
multi-valued.  The  stationary  operator  for  (1.1)  is  shown  to  be  m-accretive  in  the 

space,  so  we  obtain  a  generalized  solution  in  the  sense  of  the  nonlinear  semi¬ 
group  theory  for  general  Banach  spaces.  As  an  intermediate  step  we  sliall  show  the 
special  case  of  a  =  b  =  identity  is  re.solvcd  ns  a  strong  solution  in  every  L’’  space. 
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1  <  r  <  -oo,  and  also  in  appropriate  dual -Sobolev  spaces. 

In  order  to  motivate  the  system  (1.1),  let’s  consider  the  flow  of  a  fluid  through 
a  fissured  medium.  This  is  assumed  to  be  a  structure  of  porous  and  permeable 
blocks  or  cells  which  are  s_eparated  from  each  other  by  a  highly  developed  system 
of  fissures.  The  majority  of  fluid  transport  will  occur  along  flow  paths  through  the 
fissure  system,  and  the  relative  volume  of  the  cell  structure  is  much  larger  than 
that  of  the  fissure  system.  ‘There  is  assumed  to  be  no  direct  flow  between  adjacent 
cells,  since  they  are  individually  isolated  by  thei=  fissures,  but  the  dynamics  of  the 
flux  exchanged  between  each  cell  and  its  surrounding  fissures,  is  a  major  aspect  of 
the  model.  The  distributed  micro-structure  models  that  we  develop  here  contain 
explicitly  the  locaTgeometry  of  the  cell  matrix  at  each  point  of  the  fissure  system, 
and  they  thereby  reflect  more  accurately  the  flux  exchange  on  the  micro-sczJe  of 
the  individual  cells, across -their  intricate  interface. 


Let  the  flow  regional  be  a  bounded  domain  in  R”  with  boundeiry  T  =  dCl.  Let 
and  p{x,t)  be  the  density  and  pressure,  respectively,  at  x  G  ft  and  t  >  0, 
each  being  obtained  by  averaging  over  an  appropriately  small  neighborhood  of  x. 
At  ’ach  sucli  X  let  there  be  given. a  cell  ft*,  a  bounded  domain  in  R"  with  smooth 
boundary  Fi  =  5ft*.  The  collection  of  these  ft*,  x  G  ft,  is  the  distribution  of  blocks 
or  cells  in  the  structure.  Within  each  ft*  there  is  fluid  of  density  p{x,  y,  t)  and 
pressure  p(x,  y,<),  r^pectively,  for  y  G  ftx,  t  >  0.  The  conservation  of  fluid  mass  in 


the  fissure  system  yields  the  global  diifusion  equation 


(1.2.a) 


+  —  /(x,  t)  ,  X  G  ft  , 

in  which  the  total  concentration  p  -b  ao(p)  includes  adsorption  or  capillary  effects, 
the  function  fcj-  gives  the  permeability  of  the  fissure  system  in  the  j'"  coordinate 
direction,  5(x,  <)  is;  the  density  of  mass  flow  of  fluid  into  the  cell  ft*  at  x,  and  /  is 


the  density  of  fluid  sources.  Similarly  wc  have  within  each  cell 
(i.2.b)  +  =  ,  yen,, 
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where  to  denotes  adsorption  or  capillary  effects  and  the  function  gives  the  local 
cell  permeability.  Assume  tlie.flux  across  the  cell  boundary  is  driven- by  the  pressure 
difference  mid  is  also  proportional  to  tlie  averagc  density  p  op  that  pressure  interval. 
Thus,  we  have  the  interface  condition- 

(1.2.c)  ^ pfcj  +  ApiP  - P))  3  0  >  y  €  Tx  , 

where"!/  is-the  unit  outward  normal*  oh  Fx  and  jj.  is  the  relation  between  the  flux 
across  thednterface  and  the  density-weighted  pressure  difference  as  indicated.  The 
total  mass  flow  into  the  cell  is  given  by 

In  order  to  complete  the  dynamical  system  we  need  only  to  add  a  boundary  condi¬ 
tion  on  F 'to  (1.2.a)  and  to  postulate  the  state  equation 

(1.2.e),  P  =  <s(p) 

for  the  fluid  in  the  fissure  and  cell  systems.  Heres(:)  is  a  given  monotone  function 
(or  graph)  determined  by  the  fluid.  '' 

In  order  to  place  (1.2)  in  a  more  convenient  form,  we  introduce  the  monotone 
function 


and  the  corresponding  flow  potentials  for  the  fluid  in  the  fissures  and  cells 

u  =  S(p)  ,  V  =  S(P). 


In  these  variables  with  a  change  of  notation  the  system  (1.2)  can  be  written  in  the 
form  (1.1)  together  with  boundary  conditions  on  F  for  u  or  A(Vu)  ■  i/  and  initial 
conditions  at  4  =  0  on  a(u),  b(U).  Note  that  the  average  density  on  the  pressure 
interval  p,f>  is  given  by 


•i 


As  an  alternative  to  (1.2.c),  we  could  require  that  p  =  p  on  Tx  and  this  leads  to 
(l.l.c')  iq  i)lace  of  (l.i.c)  Finally,  we  note  that  the  classical  Forclihcimcr-t3'pc 
corrections  to  the  Darcj^  Law  for  fluids  lead  to  the  casap  =  q  —  Z/2. 

Systems  pf-the  form  (1.1)  were  developed  in  [20],  [21],  [9]  in;physical  chemistry 
as  models  for  diffusion  through.a  medium  with  a  prescribed  microstructure.  Similar 
-systems  arose  in  soil  science  [4],  [13]  and  in  reservoir  models  for  fractured  media 
l[10],  [15].  By  homogenization  methods  stidi  systems  are  obtained  as  limits  of  exact 
micro-scale  models,  and  then  the  effective  coefficients  are  computed  explicitly  from 
local  material  properties  [24],  [16],  [2j.  An  existence-uniqueness  theory  for  these 
linear  problems  which  exploits  the  strong  parabolic  structure  of  the  sj'stem  was 
given  in  [23].  One  can  alternatively  eliminate  U  and  obtain  a  single  functional 
differential  equation  for  u  in  the  simpler  space  but  the  structure  of  the 

equation  then  obstructs  the  optimal  parabolic  type  results  [17].  Also  see  [12]  for  a 
nonline^  system  with  reaction-diffusion  local  effects. 

2.  The  Variational  Formulation 

We; begin  by  stating  and  reviving  the  stationary  forms  of  our  systems.  Let  ft 
be  a  bounded  domain  iii  R”  with  smooth  boundeury,  F  =  d^.  Let  1  <  p  <  oo  and 
denote  by  L'(ft)  the  space  of  p**  power-integrable  functions  on  ft,  by  L°®(ft)  the 
essentially  bounded  measurable  functions,  and  the  duality  pairing  by 

(«,  fhw  =  f  «(x)/(x)  dx  ,  u  €  L'Cft)  ,  /  e  L^'fft)  , 

Jn 

for  any  pair  of  conjugate  powers,  ,  +  ^  denote  the  space  of  in¬ 

finitely  differentiable  functions  with  compact  support  in  ft.  W’”'?(ft)  is  the  Banacli 
space  otfunclions  in  L^(ft)  for  whidi  cadi  parti.al  derirativc  up  to  order  nt  belongs 
to  L^(ft),  and  t'Fj”'*’(ft)  is  the  closure  of  Co®(ft)  inTF’”’^(ft).  See  [l]  for  information 
on  these  Sobolev  spaces.  In  addition,  we  shall  be  given  for  each  i  6  ft  a  bounded 
domain  ft^  whidi  lies  locs)l\-  on  one  side  of  its  smooth  bonndarj-  Fx-  Let  1  <  7  <  or 
and  denote  b\'  7^  :  IF* ’’(ft,)  — »  L’fFx)  the  trace  map  which  assign.*?  bomidari'  val¬ 
ues.  Let  Tx  be  the  range  of  7x;  ilds  is  a  Banaclt  space  with  the  norm  iiulu»'<vi  in* 


a 


7x  from  Since  is  smooth,  there  is  a  unit  outward  norrnal  i/x(5)  at 

each  5  €  Fj;.  Finally,  we  define  to  be  that  closed  subspace  consisting  of 

those  (p  G  with  7j;(/3  G  K.,  z.e.,  each  7r(<^)  is  constant  &.e.  on  F^.  We  shall 

denote  by  Vj,  the  gradient  on  W^’’(nr)=and  by  V  the  gradient  on 

The  essential  construction  to  be  used  below  is  an  example  of  a  continuous 
direct  sum  of  Banach  spaces.  The  special  case  that  is  adequate  for  our  purposes 
can  be  described  as  follows.  Let  S  be  a  measure  space  and  consider  the  product 
(measure)  space  Q  =  Ct  x  S,  where  has  Lebesgue  measure.  If  17  G  L^{Q)  then 
from  the  Fubinii  theorem  it  follows  that  U{x){z)  =  U(x,z),  x  E  fi,  z  G  S  defines 
U(^x)  G  L^{S)  at  a.e.  :c  G  fl,  and  for  each  $  G  L'^  XQ): 

J  (D^(x),  ^(x))  clx  =  J  U{x,z)^{x^z)dz^  dx  =  jj  . 

Q 

Thus  L^{Q)  is  naturally  identified  as  a  (closed)  subspace  of  I/’(ft,i’(5)),  the 
Bochner  3*^  integrable  (equivalence  classes  of)  functions  from  ft  to  L'^{S).  Suppose 
X  :  ft  — »  R  is  the  characteristic  function  of  a  measurable  ft*  G  ft  and  w  G  •L’(5). 
For  each  a  >  0  we  have 

{(x,  z)  Q  Q  :  X{x)w(z)  <  a}  :=  ft*  x  {z  G  S :  w^z)  <  a)  U  (ft  ~  ft*)  x  S 

and  for  a  <  0  we  delete  the  second  term.  Thus,  X  •  xv  is  measurable  on  Q.  It  follov/s 
that  each  measurable  step  function  u  =  ^XjtWj  from  Ir^(ft,L*  (S'))- is  measurable 
on  Q,  and  hence,  belongs  to  L^{Q).  This  shows  L^{Q)  is  dense  in  and  therefore 
equal  to  L^{Q,L^{S)), 

In  order  to  prescribe  a  meeisurable  family  of  cells,  {ft*,x  G  ft),  set  S'  =  R",  let 
Q'C  ft  X  R”  be  a  given  mezisurable  set,  and  set  ft^  =  {y  G  R”  ■  {x,y)  G  Q).  Each 
fti  is  measurable  in  R”  and  by  zero-extension  we  identify  L^{Q)  >  L^{Q,  x  R^), 

and  each  L^{Q,x)  ’-r*  L^(R”).  Thus  we  obtain  from  above 

L?(g)  ^  |C7  G  L^(ft,L^(R”))  :  Uix)  G  L’(ft*)  ,  a.e.  x  G  ft}\ 

We  shall  denote  the  duality  on  this  Banach  space  by 

(u, 4')l(q)  =  y) iy\ y  e i''(Q) ,  -i>  e  £»' (Q) . 
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The  state  space  Sor  our  problems  will  be  the  product  L^  {Q.)  X  L^{Q). 

Note  that  continuously  imbedded  in  uniformly  for  x  G 

It  follows  that  the  direct  sum 

Wg  =  =  {U  G  L%Q)  :  U(x)  G  ,  a.e.  re  G  H  , 

and  /  l|C/(®)llu/*.«  dx  <  00} 

Jn 

is  a  Banach  space.  We  shall  use  a  variety  of  sudi  spaces  which  can  be  constructed 
in  this  manner.  Moreover  we  shall  assume  that  each  Ctx  is  a  bounded  domain  in 
IR"  which  lies  locally  on  one  side  of  its  boundary,  Fi,  ^d  Fi  is  a  -manifold:  of 
dimension  n  -  l-.  We  assume  the  trace  maps  7i  :  — >  Ll(Tx)  sxe  uniformly 

bounded.  Thus  for  each  U  €  W,  it  follows  that  the  distributed  trace  'y{U)  defined 
by  'y{U)(x,s)  =  {■yx{U(x)){s)^  s  G  Fi,  i  G  H,  belongs  to  I-’(Fi)).  The 

distributed  trace  7  maps  Wj  onto  Tg  =  L^{Q,,Tx)  L^(F,x)). 

Next  consider  the  collection  {Wb’(fti)  :  i  G-f^}  of  Sobolev  spaces  given  above 
and  denote  by  Wi  =  Wx’*(f2*))  the  corresponding  direct  sum.  Thus  for 

each  C/  €  Wi  it  follows  that  the  distributed  trace  7(C/)  belongs  to  We 

define  to  be  the  subspace  of  those  Cf  G  Wi  for  which  7(17)  G  Wo’^(fi).  Since 
7  :  VVi  is  continuous,  Wq’^  is  complete  with  the  norm 

llt/ll^u,  =  llt7|lw,  +  hU\\w^<p  • 

This  Banach  space  x  W,  will  be  the  energy  space  for  the  regularized  prob¬ 

lem  (1.1)  and  Wq’^  will  be  the  energy  space  for  the  ephstrained  problem  in  which 
(1.1. c)  is  replaced  by  the  Diridilet  condition  (l.l.c').  Note  fiiat  is  identified 
with  the  closed  subspace  {[yU.,U]  :  U  E  of  Wo’^(f2)  x  Wg.  Finally,  we  shall 

let  Wo  denote  the  kernel  of  7,  Wo  =  {U  G  Wg  :  yU  =  Oin  Tg}. 

We  have  defined  W^’^{Q,x)  to  be  the  set  of  ru  G  Wb9(f2j.)  for  which  7xt«  is  a 
constant  multiple  of  lx,  the  constant  function  equal  to  one  on  Fx-  Thus  WI’'^{TLx) 
is  the  pre-image  by  7x  of  the  subspace  IRTi  of  Ti.  Wc  specified  the  subspace  Wj 
similarly  as  the  subspace  of  W,  obtained  as  the  pre-image  by  7  of  the  subspace 
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L^(Q,)  of  Tg.  To  be  precise,  we  denote  by  A  the  map  of  into  Tq  given  by 

Xv(^x)  =  v{x)  -  Ij,  a.e.  X  G  v  E  A  is  an  isomorphism  of  L^{Q.)  onto  a  closed 

subspace  of  Tg.  The  dual  map  A'  taking  Tg  into  is  given  by 

=  gi^^)  =  /  gx{ix)  •  v{x)dx ,  g  eT'  ve  l^{Q)  , 

Jn 

so'we  have  X'g{x)  =  a.e.  x  6  Moreover,  when  gfj;  €  (Fj,)  it  follows 

that 

gz(lx)  =  /  gxiv)  dy  , 

*/fr 

the  iritegral:of  the  indicated  boundary  functional;  Thus,  for  gr  6  L^'  (Q,  L^'  {Tx)y:C  Tg , 
X'g  E  L^'{Q,)  is  given  by 

(2:1)  X'g{x)=  /  gx{y):dy  ,  a.e;,  x^G  H  . 

dVr 

The  imbedding  A  of  into  ^  and  its  dual  map  A  will  play  an  essential  role  in 
our  system  below. 

We  consider  elliptic  differential  operators  in  divergence  form  as  realizations 
of  monotone  operators  from  Banach  spaces  to  their  duals.  Assume  we  axe  given 
A  :  Q,  X  R"  — »  R”  such  that  for  some  1  <  p  <  oo,.sri  G  L^'{Cl),  go  G  c  and 

Co  >  0 

(2.2.a)  f ).  is=continuous  in  |*G  R”  and  measurable  in  x,  and 

(2;2.b)  (ACx,  f)  -  A(x,  if),  f -  if)  >  0  , 

(2:2.c)  A(x,0  -  4  >  Colfr-  go{x) 

for  a.e.  x  G  cuid  all  ^,Tf  G  R”. 

Then  the  global  diffusion  operator  A  :  is  given  by 

^t/.(u)  =  I  A{x,fu{x))Vv{x)dx  ,  n,t;  G  . 

in 
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Thus,  each  Au  is  equivalent  to  its  restriction  to  the  distribution 

Au  =  =  -V  •  A(-,,  Vu) 

which:  specifies  the  value  of  this  nonlinear  elliptic  divergence  operator. 

In  order  to  specify  a  collection  of  local  diffusion  operators,  + 

assume  we  are  given  B  :  Q  x  R”  — >  R"  such  that  for  some  1  <  5  <  00, 

hi  €  i.Q')')  ho  €  c  and’ (Co  0 

;(2.3.a)  S(a:,t/,  ^)  is  continuous  in  (  €  R"  and  measurable  in  (x,  y)  €  <3,  and 
;|-B(a:,y,f)l  <:c|£l’“^  +7ii(i,y)  , 

:(2.3.b)!(.B(a:,  y,  f )  -  B{x,  y,  f-  ^).  >  0  , 

{2.Z.c)  B(x, y,£)  •  f  ^  Coif  1^  -  y) 

for  a.e.  (x,  y)  €  Q  and  all  f  rf  G  R”. 

Then  define  for  each  x  G 

Bxw(v)=  f  jB(x,y,  Vyto(y))Vj,t;(y)dy.  „  i<;,u  G  W- ’’(Q,).. 

./ft. 

The  elliptic  differential*  operator  on  ft*  is  given  by  the  formal  part;  of  5*,  the  dis¬ 
tribution 

Bxiu  =  =  -Vy  •  B(x,  •,  Vytw) 

in  /Also,  we  shall  .denote  by  B  :  Wj  ^  W'  the  distributed  operaiov 

constructed  from  the  collection  r  x'G  fi};by 

BU{x)  =  Bi:{Uix))  ,  a.e.  x  G  ,  C/  G  W,  , 

and  we  note  that  this  is  equivalent  to 

BU{V)  =  fB^ {U{x))V{x) dx  ,  U.VeWg. 

Jci 

The  coupling  term  in  our  system  will  be  given  as  a  monotone  graph  which  is  a 
subgradient  operator.  Thus,  assume  m  :  R— »  R’*'  is  convex  and  bounded  by 

(2.4)  m(s)  <  C'd^i’  + 1) .  sen, 


9 


hence,  continuous.  Then  by 


m(g)  =  [  [  m(;0{x,s))  dsdx  ,  g  e  (T x)) 

Jtt  Jrr 

we  obt^n  the  convex,  continuous  m  :  L^(Q,,L^(Tx))  — ♦  Assume  ^  ^  ^  so 

that  and  consider  the  linear  continuous  maps 

A  :  -»  -^'(rx))  ,  7 :  W?;  L^O^r))  . 

Then  the,  composite  function 

M[u,U]  =  mi'fU  -Xu),  u  €  W^<^{Q.)  ,  U  G  W,  , 

is  convex  and  continuous  on  x  Wq.  The  subgradients  are  directly  computed 

by  standard  results.  [11].  Specifically,  we  have  ^  €  dm{g)  if  and  only  if 

g(^x,s)  G  dm(g(x,s))  ,  a.e,  s  6  Fi  ,  a;e.  x  6 

andwc'have  (/,F]  €  if  andonly  if  /  =  -A'(/i)  in  andF  =  7-(/i) 

in  W'  for -some  n  G‘dm('yU  —  Xu). 

The  following  result  giv^  ^sufficient  conditions  for  the  stationary  rtgulgrizti 
problem  to  be  well-posed. 

Proposition  1.  Assume  1  <  p,g,  ^  and  define  the  spaces  aud  operators 

A,7  as.above.  SpeciGcally,  the  sets  {Qx  '  x  G  fi}  are  uniformly  bounded  with  smooth 
boundm:ies,  and  the  trace  maps  {■jx}  are  uniformly  bounded.  Let  thedunctions  A^  B , 
and  m  satisfy  (2.2),  (2.3)  (2.4)^  and  assume  in  addition  that 


(2.5) 

m(s)  >  coisl®  ,  s  G  E  . 

Then  for  each  pair  f.  G  W~^’^*(n),  F  G  VV'  there  exists  a  solution  of 

(2.6.a) 

u  eWo^’’’(Q)  :A(u)-y(p)  =  f  in  W-^’P\Q) 

(2.6,b) 

UeWq’.  B{u)  +'r'{n)  =  F  in  w; 

(2,G.c) 

p  G  Li' p,  Li' {Tx))  :  P  G  dm{jU  -  Au)  . 

Fjr  any  sucli  solution  wefhave 


f  IJL^x,s)ds  =  ,  a.e.  x  G  f2  , 


where  lx  denotes  the  constant  unit  function  in  W^>^(Q,x). 

Proof.  The  system  (2.6)  is  a  “pseudo-monotone  plus  subgradient”  operator  equa¬ 
tion?  of  the  form 


(2.6') 


[rx,  U]-  €  ’^(H)  X  :W,  :  for  all  [u,  V]  G  ’''(fl)  x,  W, 

Au{v)  -h  BUiy)  -f  dM[u,U]{[v,  V])  B  f{v)  +  F(y)  . 


It  remains  only  to^verify  a  coercivity  condition,  namely, 


Au{u)  +  BU{U)+M'rU->^n) 

^  ^  ll^llu^-cny+ilC^llvv,  “ 

as  llw|lwi.?(n),+J|t^llvv,  ->  +po,. 

Choose  k  =.rhax^{|y„| :  y^G  ft*,  x  G-ft}  Jmdlet  i/*  =  (^'J , . . . ,  i/”)  be  the  unit  normal 
on  F*.  For  u  G  W^’^(ft*)  we  have  by  Gauss’  Theorern 

[  (|u|’ -h  yn9lv|?“^5„t;)  =  /  5„(y„lv(y)|®)  dy. 

jQr 

=  f  x/"(s)sfi|7*u(s)|^  ds  . 

Holder’s  inequality  then-shows 

and- from  this  follows 

IMILcn,)  <  2fc|l7xv|li,(r,)  + 

by  Young’s  inequa,lity.  From  here  we  obtain 

(2.9)  coWVWl^Q^  <  +  l|Vy^lll,(Q)  ,  y  €  W,  . 
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Thus  from  the  a-pfiori  estimate 


Au{u):+  BU{U),  +  M(-yU  -  Xu)  > 

(2.1Q)  -  ||Ao||i..«) 

^coW-yU  -  ,  «  e  wl-^ia)  ,-UeW,, 

the  Poincare-type  inequality  (2:9)  smd  the  eqiiivzdehce.of  ||V«||£,;(f2)  with  the  riorm 
on  we  obtain  the  ccercivity  condition  (2.8).  Specifically,  if  (2.8)  is 

bounded  by  then  (2.10)  is  bounded  above  by 

^  +  l|Vy^^||L<(Q).+  -  Av?||L»(n,'L«(r,^^  4-  llAu||£,*(n))  , 

and  the  last  term  is  dominated  by  the  first.  This  gives  an  explicit  bound  on  each 
of  these  terms  and^  hence,  on  +  IK'ilw,* 

Finally,  we  apply  (2.6.b),  to  the  function  6  given  by  K(.'?fy)  =  v(x)Tor 
some  v'€  i^(n),  and  this  shows 


fi{^v)  =iF;v)' 


since  BU{V)  =  0,  and  thus 

/  X' fi{x)v{x).dx  =:fi{Xv)  =  fJi{jv)  =  I  '(J^('x),l)u(x)  di  . 

Jn  Jq 

The  identity  (2.'T).'now  followsTrom  (2.1); 

For  the  more  general  case  of  the  dcgenerocie  ^fofioncry  yroi/em  corresponding 
to  (1.1),  we  obtain  the  following  result. 

Corollary  1.  Let  (/? :  H— » IR.^  and  $  :  IR— »  be  convex  and  continuous,  with 

9?(0)  =  $(0)  =  0,  and  assume 

(2.11)  vj(s) < a(|sp -h  1)  ,  <i>(s) ^ c(|sp 1) ,  sen. 
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jFbr  each  pair  f  €  W  (fi),  F  £  VV^,  there  exists  a  solution  of 
(2.12.a)  u  €  ^0  4(u)  -  A'(,/ij.=  /  in  (Q) 

(2:.a2.b)  U  6  Wg  :  6  5(C7)  +  =  F  in  W' 

(2.12.c)  p  €  d7n(‘yU  —  Au)  m  jV  \yi-~  L^  (rx)),  j  3pd 

(2.12.d)  a  e.d(p{u)  in  ,  o>$(J7)  in  . 

For  any  such  solution  we  haye 

(2.13)  /  b{x\y)dy+  j  //(x,5)6'5  =■=  (F(x).  li)  ,  a.e.  x  G  ft, 

Jn*  ^Er 

Proof.  This  follows  as  above  bui  with  the  continuous  convex  function 

'?[«, =  y  ¥’(u(x))  dx+  y  /  $(CZ(x,y,)>dy  rfx+m(^j7-Au) ,  [u,  17]  €  Wi’'’(i7)x-W, 
7a  Ju.ihr 

The.subgradient  c^.be  coihputed- termwise -becauseithe  three  terms  axe. continuous 
oh  JD^(n)i;X^.((5),-  X’(rx)), -respectively. 

Remark.  The  lower  bound  (2.o)  on  m(*)  may  be  deleted  in  Corollary  1  if  such 
a  lower  estimate  is  knoivu  .to  hold  for  It  is  also  unnecessary  in  the  matched 
microstructure  model;  see  befow. 

In  order  to  prescribe  the  boundary  condition  (l.l.c)  explicitly,  we  develop 
an  appropriate  Green’s  formula  for  the  operators  Z?*.  Note  that  v/e  can  identify 
X®  (fix)  C  since  kFo’^(flx)  is  dense  in  so  it  is  rne^ingful  to 

define 

Bx  =  {«;  €  W^’^(ax)  :  e  ■ 

This  is  the  domain  for  the  c65troct  Green’s  Theorem. 

Lemma  1.  There  is  a  unique  operator  dx  :  Di  —*T^  for  which  BxV)  = 
for  all  XU  G  Dx-  That  is,  we  have 

(2.14)  Brw{u)  =iBxXu,v)i,{ar)  d-idj,xu,^xv)  ,  v  e  , 
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:for  every  ty  G  D'i- 

Proof.  The  strict  morphism  jx  of  onto  Ti  i:as  a  dual  whidi  is  an 

isomorphism  of  onto the  annihilator  in  of  the  kernel  of 

7x.  For  each  w  G  Dx,  the  difference  Sxio  —  BxW  is  in  Wq’^(Q,x  so  it  is  equal  to 
y'j.{dxw)  for  a  unique  element'  dxV)  €  Tj. 

Remark.  The  identity  (2.14)  is  a  fceneralized  decomposition  of  Bx  into  a  partial 
diffeiehtial  operator  on  9/^  tmd  a  boundary  condition  on  Fx.  If  .Tx  is  smooth,  Ux 
denotes  the  unit  outward  normal  on  Fi,  and  if  JB(a:,  *,  Vyty)  (fir)]”,  then 

w  E  Dx  and  from  the  cli^sical  Green’s  Theorem  we  obtain 

Bxw{v)-  (^SxW,v)L(a^y  =  I  B(xfS^Vyw)i/x(s)‘^v(s)ds  ,  v  G 

JVx 

Thus,  dxW-=  B(x,',V.yw)  T  ifx  is,  thedhdicated  hormad  derivative  ih  I/*'(Fx)  when 
B(x,  •,  Vyw)  is  as  smooth  as  above,-,  and.  so  we  can-  regard  dxV)  in  general  as  an 
extension  of  this  nonlinear  differential  operator  on  the  boundary, 

The  formal  part  oi  B  :  W,  -*  W'  is  the  operator  B  :  W,  -+  Wq  »jveh  by  the 
restriction  =B(17)  =  5l7|vvo-  Since  Wo  is  dense  in  L^{Q)  we  can  specify  the  domwn 

D  =  {U€Wg:B.{U)  e 

on  whicli  \ye  obtain  as  before  a  distributed  form  of  Green’s  theorem. 

Lemma  2.  There  is  a  unique  operator  d  :  D  —*  Tg  such  that 

B(c/)( V)  =  (BCD-),  +  (ay,  tK)  ,  Ue.D,v&  w, . 

Proposition  2.  Let  the  Sobolev  spaces  and  trace  operators  be  given  as  above.  IVe 
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summarize  them  in  the  following  diagrams 


U 


L^rx)  L^m 

U  U  U 

^  Ti  r, 

U  L  U-  -^L 

Wj>?(ax)  IR-lx  Wi  — -^ 

U  U  T 

>  {0}  Wo  — {0}: 

in  which ‘ji  is  the  restriction  of  j  to  Wi.  Wq  Wo  dense  in  L^(Slx),  L^(Q)> 

respectively.  Let  operator  Bx,  x  G  fl,  and  B  be  given  and  define  their  formal  parts 
Bx,  B  as  above.  Then  construct  \he  domains  Dx,  D  and  bound^y. operators  dx,  d 
as  in  Lemma  1  ajid  Lemma  2,  f&.pectively.  It  follows  that  for  any  U  €  W,, 

(a)  BU(x)  =  Bx{U(x))::ih  Wo^^flx)'  for  a.e.  x  €  fl,  and  U  e  D  if  and  only  if 
U{x)  G  Px  for  a.e.  x  G  .P  and  x  BxU(x)  belongsHq  L^  {Q); 

(b)  for  each  U  G  D, 

dU{x)  =  dx{U{x))  in  T' for  a.e.  x  6  ft 


BU  =  BU  +  7'i(^’dU,)  in  Wi  , 


andifor  eaxdi  V  €  Wi  we  have 


f  BxU(x)(V(x))  dx  =  f  BxUix)V{x)dy  dx+  f  {dxU{x),  lx)iliV){x)dx. 
Jci  JQ  dn 

Proof,  (a)  For  F  €  Wo  we  obtain  from  the  definitions  of  i?.  B  and  Bx,  respectively. 

f  BU{x)V{x)  dx  =  /  BUiVfdx  --=  /  BxU(x){y(xyfdx  =  I  BxU{x)V{x)  dx  , 
Ji'i  Jn  Jn  dn 


and  so  the  first  equality  holds  since  Wq  =  L'^  (fi,  The  characterization 


of  D  is  immediate  now. 


(b)  For  V  €  Wg  we  obtain  from  the  definitions  of  7, 5,c>i,  respectively,  and  (a) 

I  dU{'y^Vix))  dx=  f  dU('ry)dx=  (  {BU  -  BU){x)V{x)  dx 
Jo.  Jo  Jo 

=  f  dx{U{x))yxV{x)dx  . 

Jo 

Since  the  range  of  7  is  7^'  =  the  first  equality  follows.  The  second  is 

immediate  from  Lemma  2  since  on  Wi,  7  =  A  o  71  and  7'  =  7iA',  and  tlie  third 
follows  from  the  preceding  remarks. 

Corollary  2.  In  the  situation  of  Corollary  1,  /  G  (ft)  and  F  G  ((3)  if  and 
only  if  Au  €  L®* (ft)  and  B(JJ')  €  L^"{Q),  and  in  that  c^e  the  solution  satisfies 
almost  everywhere 

a(x);€  %(u(a:))  ,  d(x)  +  A«(x)+ -r  6(1, y) dy  = /(x)  +  /  F{x^y)dy  ,  x  G'Cl  , 

Jo^  JOx 

ti(s)  =  0  ,  s  €  r  , 

6(i,y)  6  5$(?7(x,  y))  b{x,  y)  +  BU.{x,y)  =  F{z,  y)  ,  y  €  ft* 

n{x,  s)  e  dm('yU(x,  s)  —  u(x))  ,  dx{U{x)){s) +  fi{x^s)  =  0  ,  5  €  F*  . 

Finally,  we  note  that  corresponding  results  for  the  stationary  matched  mi- 
crostructufe  model  aj.e  obtained  directly  by  specializing  the  system  (2.6')  to  the 
speice  W5 This  is  identified  with  {[yUiU]  :U  G  Wj’^}  as  a  subspace  of  Wq ’^(ft)  x 
Wg,  and  we  need  only  to  restrict  the  solution  [u,  17);  and  the  test  functions  [u,  V], 

V  =  7V,  to  this  subspace  to  resolve  the  matched  model.  Then  the  coupling  term 
M  does  not  occur  in  the  system;  see  the  proof  of  Proposition  1,  especially  for  the 
coercivity.  These  observations  yield  the.following  analogous  results  for  the  matched 
microstructure  model. 

Proposition  1'.  AssUxne  1  <  p,?,  ^  ^  and  define  the  spaces  and  operators 

A,  7  as  before.  Let  the  functions  A,B,  and  m  satisfy  (2.2),  (2.3)  and  (2.4).  Then 
for  each  pair  f  G  l'F“b?’’(ft),  F  G  W'l  there  exists  a  unique  solution  of 

(2.15.a)  =  /  +  in  ll'-''^'(ft) 
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(2.15.b). 

(2.15.C) 


G  Wy  :  B{U)  =  F  in  Wq 
-yU  =  Xu  in  Li{a)  C  T,  . 


Corollary  1'.  Suppose  9?,$  axe  given  as  before  and  assume  (2.11).  For  f,F  as 
above  there  exists  a  unique  solution  of 


(2.16.a);  u  €  Wo :  a  +  (6, 1)  +  A(«)  =  /  +  {F,  1)  in  W-bp'(a) 
(2.16.b)  UeWi:b  +  B(U)  =  F  in 
(2.16.C)  jU  =  Xu  in  L^(n)c% 

(2.16.d)  ae  dcp(u)  in  L^'(a):,  b€d^(U)  in  L^\Q)  . 


In  addition,  f  &  (Q.)  and  F  €  L^'{Q)  if  and  only  if  Au  G  L^’(Cl)  and  B(U)  G 

(Q),  and  iwtbatcase  the  solution  satires  almost  everywhere 

a(x)  G  d(p(u(x))  ,  a(i)  +  Au(x)  +  /  6(1,  y)  dy  =  /(x)  +  f  F(x,  y)dy  ,  x  G  ft  , 

Jstx 

«(«)  =  0  5  6  r , 

Ka;,  y)  G  9$  (tr(x,  y))  ,  6(x,  y)  +  BU (x,  y)  =  F(x,  y)  ,  y  G  ftr  , 


17(x,s)  =  u(x)  , 


5  G  r* . 


Remark.  For  the  very  special  case  of  p  =  y  >  2  and  a(u)  =  u,  b{U)  =  U  in 
the  situation  of  Proposition  1  it  follows  from  [6]  or  [19]  that  the  Cauchy-Diriclilet 
problem  for  (1.1); is  well-posed  in  the  spaceXP(0,T;  Wq ’'’(ft)  x  Wp)  with  appropriate 
initial  data  «(x,0),  U{x,y,0)  and  source  functions  /(x,i),  F(x,y,t).  A  similar 
remark  holds  in  the  case  of  Proposition  1'  for  the  matched  model  with  (l.l.c'). 
These  restrictive  assumptions  will  be  substantially  relaxed  in  the  next  section. 


Furthermore,  variational  inequalities  may  be  resolved  for  problems  correspond¬ 
ing  to  either  the  regularized  or  the  matched  microstructurc  model  by  adding  the 
indicator  function  of  a  convex  constraint  set  to  the  convex  function  'P.  Thu.s  one  can 
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hauidle  sucli  problems  with  constraints  on  the  global  variable  u,  the  local  variables 
U,  or  their  difference  Au  —  7!/  on  the  interface. 


3.  The  L’‘-Operators 

Assume  we  are  in  the  situation  of  Proposition  1.  We  define  a  relation  or  multi- 
vedued  operator  C2  on  the  Hilbert  space  x  as  follows:  C2 
if  and  only  if 

(s.i.a)  ueL-{p)r\wl'\a)zA{u)-\'pL  =  feL\a) 

(3.1.b)  U  €  L\Q)  n  W,  :  B{U)  +  7V  =  F  €  L\Q) 

for  some  /i  €  5m(7C/  —  Au)  in  L^' {CL,  L^' {Tx))- 

Thus,  G2  is  the  r«triction  of:(2.6)=to  x  L^(Q).  Note  that'  A'/i  €  L^{Cl)  hy 
(2.7). 

Lemma  3.  Jf  t :  iT-*  II  is  monotone,  Lipschitz,  and  cr(0)  =  0,  then  for  each  pair 

Caluj 9  [fj,  Fjl ,  j  =  1, 2, 

there  follows 


(/i  - /a,  +  {F1-  F2,  cr(l7i  -X/2))^2fov>  0 


LHQ) 


Proof.  Since  <j  is  Lipschitz  ^d  <r(0)  =  0,  we  have  a(ui  —  U2)  6  Wo’^(^) 
afUi  —  U2)  6  W,.  Also  the  ch^n  rule  applies  to  these  functions,. so  we  compute 

{Auj  —  .4u2,  (j(ui ,  U2))  =  /  (A(x,  Vui)  —  .4(1,  VU2))  V(ui  —  U2)a'(ui  -  ua)  dx  , 

Jsi 

{BUv- BUo,  a{Ui  -  Uo))  = 

/  [  — B(x,y.yyU2))^,j{Ui -U2)<r’{Ui -  U2)dydx  . 

Jn  in. 


IS 


Both  of  these  arc  noii-negative  because  of  (2.1.b),  (2.2.b)  and  a'  >  0.  The  remaining 
term  to  check  is 

-  /i2),  <^(W1  -  +  {l  {>M  -  /i2),  o(Ui  -  Ut)) 

=  /  /  ^)  ”  5))  -  '■fUz)  -  o-(A«i  -  Aua))  ds  dx  . 

J(i  Jr* 

Since  cr  is  a  monotone  function  and  dm  is  a  monotone  graph,  this  integrand  is 
non-negative  and  the  result  follows. 

As  a  consequence  of  Lemma  3  with  a(s)  =  s,  the  operator  C2  is  monotone  on 
the  Hilbert  space  L^(Q)  x  L^{Q).  Moreover,  we  obtain  the  following. 

Proposition  3.  The  operator  C2  is  maximal  monotone  on  x  L^(Q).  Let 

j  :  IR  ^  IR^  be  convex,  lower-senu-continuous,  and  j(0}  =  0.  If  dm  is  a  function, 
then  C2  is  also  single-valued  and 

(3.2)  (C2[ui,  tfj]  -  C2(«2,  U2],  ® 

for  any  selections  <fi  €  dj(uj  —  U2)  in  L^(Q)  and  02  €  dj(Ux  —  U2)  in  L^(Q). 


Proof.  To  show  C2  is  maximal  monotone  it  suffices  to  show  that  for  any  pair 
[/,  2^]  €  L^(fl)  X  L^(Q)  there  is  a  solution  of 

(3.3.a)  u  €  O  :  u  -f  >1(«)  -  A'(^)  =  /  in  , 

(3.3.b)  U  €  L^(Q)  nW,:U  +  B{U)  +  =  F  in  W; 

(3.3.c)  p  €  (^x))  :  p  ^  dm^jU  —  Au)  . 


The  existence  of  a  (unique)  solution  of  (3.3)  follows  as  in  Proposition  1,  but  by 
considering  the  pseudo-moiiotone  operator  [A,B\  on  the  product  space  H 

TPj’*’(f2)  xi^(Q)nVV,  and  the  convex  function,  |I|Mlit3(ft)+5liC^||^2(Q)+m(7£/— Aa). 
on  that  space. 

To  establish  the  estimate  (3.2),  we  consider  the  lowcr-scmi-continuons  convex 
function 

(3.4)  J[u,  U]  =  (^iu{x))  -i-  j{U(x,y))  dy^  dx  ,  (w,  U]  6  L-(fl)  x  lf^{Q)  . 
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The  subgradient  of  j  is  given  oh  this  product  space  by 
a  =  [<71,(72]  G  dj[u,U]  if  and  only  if 

a[v,V\  =  ^  a2{x,y)V{x,y)d‘^dx  ,  \v,V]  €1-^(0)  xL‘^{Q))  , 

where 


cri(x)  €  9j(u(a:)),  a.e,  a:  €  ,  <r2(x,y)  G  5jf  (l7(a:,  y))  ,  a.e,  '(x,y):G  Q  . 

The  Yoshida  approximation  jg  of  J  is  given  as  in  /3.4)  hut  with  j  replaced  by  jg. 
Since  the  derivative  of  je  is  Lipschitz,  monotone,  and  contains  the  origin,  it  follows 
by  Lehima  3  that  the  special  case  of  (3.2)?with jV  is  true;  Thus,  C2  is  Sj-monotohe 
[7].  and  the  desired  result  follows,  since  the  single-valued'  G2  equals  its  minimal^ 
section. 

We  define  the  realization  of  (2.6)  in  Lr(fl)  x  L'’\Q)y  1  <  r  <  00,  eis  follows.  For 
r  >  2,  C,.  is  the. restriction  of  C2  to  X  L'^{Q),  andTor  1  <  r  <  2,  Gr  is  the 

<:losure  in  L^{Q,)M  L^{Q)  of  C^. 

“Corollary  3.  The  operator  Cr  is  m-accretive  in  (Si)  X  L^(Q)  for  1  <  r  <  00. 

Proof.  Let  (/  sC2)([iijyUj])  3  [fj,Fj],  j  =  1,2,  and  assume  [fj^Fj]  €  L^{Sl)  x 
^L’'(Q)  ‘iir  >  2.  Set  j(s)  =  Is]'’,  s  G  R.  Kdm  Proposition  4.7  of  [7]  it  follows  that 

l|[ui  -U2,Ui  -U2]\\L’-^n)xL'-{Q)  ^11 1/1  -f2^Fi  p2]l|L’-(n)xL''(Q)  • 

Taking  1/2 ,  •P’2]  ^  [O1O]  we  see  that  Zr*'(fi)  x  T’^iQ)  is  invariant  under  (J  +  6:G2)"S 
and  then  the  estimate  shows  this  operator  is  a.  contraction  on  that  space.  We 
-have  Rg(I  -f  eC,.)  =  L'^(Sl)  x  L^fQ)  directly  from  the  definition  for  r  >  2,  and' for 
1  <  7*  <  2,  Rg(I  -f  eGr)  D  L‘^(Sl)  x  L^(Q),  which  is  dense,  so  the  result  follows 
easily. 

Remarks.  The  Cauchy-Dirichlet  problem  for  the  regularized  model  (1.1)  is  well- 
posed  in  L’'(Sl)  X  L''(Q)  when  a(u)  =  u,  h(U)  =  [/,  and  ?•  >  1.  This  follows 
from  Corollary  3  and  the  theory  of  evolution  equations  generated  by  777-accretivc 
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operators  in  a  uniformly  convex  Banacli  space.  For  example,  from  [18]  we  recall  the 
following: 


if  f  G  T'F^’^(0,r;X)  and  wo  £  D{Cr),  where  Cr  is  m-accretive  on  the  uniformly 
convex  Banach  space  X,  then  there  exists  a  unique  Lipschitz  function 

w  :  [0,1]  — »  X  for  which 

w' (t) +  Cr{w(t))  B  j{t)  ,  a.e.  t  £{0,T)  , 

w(t)  £  D{Cr)  for  all  i  £  [0,  T]  ,  rmd 
w(0)  =  Wo  • 

See  [3]^for  details  (Theorem  III.2.3)  and  references;  By  applying  this  result  to  the 
operator  given  in  X  =  L^(Q,)  x  1  <  r  <  op,  we  obtain  a  generalized 

strong  solution  to(<)'  =.[it(<)',X/'(t)]' of  the  system 

+  Au(x,  <)  +  — 1’^ --  dy  =  /(i,  0  +  F{x,y,  t)  dy  ,  a:  €  ft  ,  f  €  (0,  T)  , 

u{s,  i)  =  0  ,  s  €  r  ,  ■ 

+ =  F{x,y,t)  ,  y€&.  , 

p,{x,s,t)  £  dm{U{x,  s,t)  -  u{x,t))  ,  dxU{x,  s,t)j+  fi{x,  s,t)  =  0  ,  s  e  Tj;  , 
u(fCjO)  =  Uo{x)  ,  U{x,y,0)  =  Uo{x,y)  . 

The  restrictions  on  the  data  f(t)  =  [f{t),  F(t)]  and  u>o  =  [wo>  Uq]  can  be  considerably 
relaxed  in  the  Hilbert  space  case  r  =  2  [7). 

By  applying  Proposition  1'  similarly,  it  follows  that  corresponding  results  for 
the  matched  model  are. obtained.  Thus  one  obtains  a  generalized  strong  solution  in 
L''(ft)  X  L’'(Q),  1  <  r  <  oo,  of  the  system 
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+  Au{x,  t)  +  j_  =  /(®>  0  +  J  y^i)di ,  X  eo, ,  te  (o,  t) 

u(^$,t)  =  0'  5  €  F  , 

+  BUix,y,t)  =  ni,v,()  ,  V  6  ft.  , 

U:{x,S,t)  =  u{x^t)  y  5  G  Fi  , 

u(a:,0)  =  uo(a:)  ,  U(x,y,0)  =  Ui{x,y)  . 

This  follows  as  above  from  the  analog  of  Proposition  3  and  Corollary  3. 

We  return  to  consider  the  fully  nonlinear  model  (1.1).  The  generator  of  this 
evolution  system  will  be  ;obtained  by  closing  up  the  composition  of  C2  with  the 
inv^se  of  [d(p,  9#]  in  L^(p)  x  Thus,  we  begin  with  the  foilOwing. 

Definition.  C(a,  6]  3  if  Catu,!/]  3  (/,F|  and  a  €  dip{u)  m  L‘^{Cl\  b.  e  d$(U) 
in  L^{Q)  ioT  some  pair  [u,f7]  as  im(3.1). 

Lemma  4.  The  operator  C  is  accretive  on  L^{Q,)  x  {Q)-if  either  dm  is  a  function 

or  if  both  d(fi  and  9$  are  functions,  < 

Proof.  Let  e  5>  0  and  suppose  that  (/  +  eC)[ay,  6j]  3  [fj,  P!j]  for  j  =  1, 2.  Thus  we. 
have  eC2[uj,C/y]  9  [/j  —  aj,Fj  —  hj],  aj  €  d(p{uj),  bj  €  d^{Uj)  as  above.  First  we 
choose  <7(5)  =:  sgh^ (s),  the  Yoshida’ approximation  of  the-maximal  monotone  sgn"^, 
apply  Lemma  3  and  obtain 

(oi  -  a2,sgnj‘(ui  -  +  (61  -  62,sgn+(C/i  -  U2)):,^^^Qy 

<  l|(/l  -  /2)‘^lll,»(0)  +  “  ■^2)^1|l>(Q)  ' 

If  9^  and  9$  are  functions, 4hen 

(ai'‘-  a2)sgn3;(a,  -  U2):=  ;(ai  -  02^'^  , 

{K-h2)sffitiU, , 

so  letting  6  -+  0  gives 

(3.5)  ll(a,-d2)+|U,(„,+||((.,-62)+||i.(Q)  <  l|(/i-/2)+||2,.(n)+ll(0-^2)+lli,.w)  • 


22 


The  same:holds^or  negative  parts,  so  it  follows  that  (I+eC)“^  is  an  order-preserving 
contraction  with  respect  to  B\Q,)  x  L^(Q)  for  each  e  >  0. 

Next  we  suppose  drn  is  a  function.  Choose  J(s)  =  5"^,  sO  that  dj  =  sgn^,  and 
then  set 

cri(x)  =  sgn^:(ttj  -  U2  +  ai  -  aa)  €  sgn'^(ui  -  ua)  H  sgn'^Xaj  -  ca)  , 

<^2(^1  y)  =  sgn^(l7i  -  U2  +  bi  -  6a)  6  sgn'^(17i  -  U2)  D  sgn'^(6i  -  6a)  . 

Proposition  3  applies  here  to  give  (30).  A  similar  estimate  for  negative  parts  yields 
the  result. 

Although  G  is  not  accretive  on  V  for  1  <  r,  we  can  obtain  estimates  when 
the  graphs  dipyd^  aremot  too  dissimilar. 

Corollary  4.  If  (/  -f  eC)[a,  6]  3  [fyF]  with  e  >  0,  then 

||a'^|Uoo(n)  <rnax(ao(^*),|l/'‘‘|lL.»(n))  , 

(3.6) 

l]^'‘’lli:.~(Q)'^'hiax(6o(fc),||F^||L~  , 

where  =  max(a^^(ir/+l|L~ ),  6^^(i|i^^^ 

Remarks.  Here  gq  is  the  minimal:  section  (5v?)o,  is  the  minimed  sectioii  of 
(9^)“^,  and  bo,bQ^  are  defined  similarly  from  5$.  Specifically,  we  obtain  an  ex¬ 
plicit  a-priori  bound  on  l|a'*';||L«>(n)  and  ;{|6+|1£,oo(q)  when  ||/■*^|| £«>(«)  ^  Ryid'p) 
2uid  ||i^'^|lL«>(Q);  €  Rg{d^),  By  siinifar  estimates  for  negative  parts^  we  obtain  ex¬ 
plicit  esUmates  on  ||a||x~(h)  and  ||6||2ioo(Q)  for  any  pair  /  €  i®®(ft),  F  G  L°°{Qr) 
if  Rg(d(p)  =  R,  and  i2^(5$);  =  11  or  (trivially)  if  both  Rg{d(p)  and  Rg(d^)  are 
bounded  in  R.  Finally  we  note  that  in  the^special  case  (p  =  $,  we  obtain  - 

max(||a+||Lco(fi),  ||6+I|£,«o(q))  <  max(|I/+|Uco(n),  l|i^'*'|U~(Q))  • 

Proof.  By  the  choice  of  k  >?0.we  have 
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for  some  pair  Subtract  these  from  the  operator  equation,  multiply  by  -either 

sgn^£(u  -  k)  ,  sgnjp  -  k) 

or  by 

sgnj(a  -  +  it  -  k)  ,  sgn^  (6  -  £2  +  U  -^k)  , 

depending  on  whether  and  are  functions  or  dm  is  a  function,  respectivelj'. 
Apply  Lemma  3  and  let  6  — »  0  or  apply  Proposition  3,  respectively,  to  obtain 

il(a  -  «i)+lll.-(fl)  +  IK*’  -  «2)+|U.(«)  <:||(/'*'  +  IKf”^  -  lli.(«)  . 

The  right  side  is  zero,  so  tliAresult  follows. 

Proposition  4.  (Moser)  Let  (it,  C/)  €  WQ’^(f2)  x  VVj  be  a  solution  to 

A(u)  -  X' 3  f  in  , 

B(Uj  +  7'/x  BT  in  W;  , 

(i  €  .dth('yU  —  Au). . 

1)  If{f,F).e  L^'(Q,)  x  L^\Q)  withr^^>  J,  and 

(2.2.c')  A(x,f )'-f  >  coiei'’  -  ffoXx) 

where; go  ^  L^'{Q),  then  L°®(ft). 

2)  If,  additionally,  F  €  L°°[ayl/’ (Qi)]  with,  t'  > 

(2.3.c')  B{x,  y,^)  •  I* >  col^l’  -7io(x,  y) 

where  IiQ  €  [(2 j  L-  (fli)],  and  m  satisfies  the  growth  condition  (2.5)  and  m(O)  =  0, 

thenUeL°°{Q). 

Proof. 

1)  Estimate  (2.7)  of  Proposition  1  shows  that.  X'/i  6  L^' (Q.),  so  that 

^(tx)  =  /-AV  =  /eL"'(fi)  . 
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Lemma  3  of  [22]:  can  how  be  used  to  conclude  u  € 

2)  Define  U  —  U  —  ul.  Since  B(JJ)  =  follov  s  that 

B{U)  +  7V  =  in  W'  ,  /i  €  dm{^U)  , 

and  for  almost  every  x  £  ^1,  and  every  V  €  W, 

(*)  /  B{x:r,VyU{x))  ■VyV(:x)+  t  (x(xy^Vix)  =  f  F(x,-)V{x)  , 

Jcir  Jr^x 

with  iJ.{x)?E.  dmi^Uix)).  We  now  use  Moser  iteration  with  (*)  to  conclude  |lt^(®)l|L~(n*)  ^ 
C,  where  C  is  to  be  chosen  independently  of  x  6.  Q>. 

If  U^x)  6  (r  =  q  suffices  for  the  first  iterate),  define  s  =  1 

(7+:^  =  1).  Let  if  €  C^(IR-)  satisfy  H{s)  =  |s|*  if  [s]  <,so,  if  affine  for  Js]  >  sq?  and 
define  G^s)  =  d|.  Since  if  has  line^  growth,  it  follows  that  G(  17.)  €  W,. 

Substituting  G(I7)  for  V  in  (.♦,)  gives 

f  Kt)7G(K(x))  =.  /  F(^,-)Q{U(,x)): . 

The  first  term  of  the  above  is  bouhded  below  using  (2.3-c').  To  estirnate  the  second 
term,  use 

(i)  >  7n(Cf  )  (as  m(0)  =  0),.and 

(ii)  sgn(U)  =  sgn(G(Cf))  (so  that  G(U)/U  >  0  when  Cf  ^  0), 
to  get 

IxGiU)  =  fiU  G(U)IU  >  m{U)  G{U)iU 

>  colif|«  G{U)!U  =  co\U\>-'^\G{U)\  . 

Co  f  \V:yU\<^G'iU)  +  co  f  Fr'lG(Cf)l  <  /  ■FGiU)  +  hoG'{U). 

JVx  Jcix 

The  first  term  inay  be  written  as  |yyif(i7)|^  which,  using  the  Sobolev  embedding 
theorem,  is  bounded  below  by 

cicwmr  , 
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where  e  >  0  can  be  chosen  zn’bitr^iJy  small  (see  equation  (2;9^).  The  right  hand 
side  is  bounded  using  Holder’s  inequality. 

c(£)iiff(K)r^  +/  ipr>iG(£?)i-£|H(g)i« 

L^-na.)  •'r. 

when  So  ooj  |C/|*,  .and  |W|*“^|Ct(17;)|  ^  i7(r)|t/’.|*^,  where  7i(r)  =  ^s®  = 

^(1  +  H  e  is  chosen  as  e  =  min it  follows  that 

ll^llLM("/n-o(n,); $  maxflvlL?||L-(nx)]i  •. 

The  result  now.  followsiby  iteration  of  the  above  estimate. 

'Theorem  1.  Assume  the  hypotheses  of  Proposition  1,  Corollary  1,  Lemma  4,  and 
Propostioii  4.  Also,  assume  that  Rg{d(p)  and  Rg(d^)  are  both  bounded  or  that 
both  are  equal  to  K,.  Then  C,  the  closure  of  C  in.L\(Q,)  x  L\Q),  As  m-accretiye. 

Proof.  Let  /  €  L'^(fi)>and'F*€  L°^(Q^.  CbroUaxy  1  asserts  there  is  a  elution  of 
(2.12).  ILthe  graphs  5^.and  5$  have  bounded  range,  then  a  €  L°°(f2),  b  € 
and  it  follows  from  Proposition  4  that  u  €  L^(A)  ^d  U  €  I?{Q).  This  shows 
C2[u,  17]  9  [a  —  f,b  —  J?],  so  (J  +  C)([a,6])  9  [/,^].  Thus,  Rgfl  +  C)  is  dense  in 
and,  hence,  equal  to  (12)  x  L^{Q). 

If  therranges  of  593  and  5$  equal  II,  then  by'Corollafy  4  aiiy  solution  satisfies 
:i]alU~(n)  <  ■K'  ,  ||6||£,«(Q)  < /v  , 

where  K  depends  on  f  and  i^.  Replace  9^,5$  by  the  appropriately  truncated 
The  solution  with  these  tfuncated  graphs  then  is  a  solution  of  the  equa¬ 
tion  with  the  original  graphs,  so  we  are  done. 

Corollary  5.  Under  the  hypotheses  of  Theorem  1,  problem  (1.1 )  has  a  unique 
generalized  solution  {a^b)  6  C[0,r;L^(12)  x  L^{Q)],  provided  the  data,  satisfy 
if,  F)  €  (0,  T;  V  (fi)  X  {Q)\,  and  (a(0),  6(0))  €  ^Q. 
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Thip  follows<ffom  the  Grandall- Liggett  Theci:em'[8],  which  is  proved  by  showing 
that  the  step  functions,  iConstructed- rrom  solutions -to  the  differencing 

sclieme 

(r  =  ^ )y  converge  uniforrnly  when  the  operator  C  is  mraccretive.  Benilan  [5]  proves 
that  these  generalized  solutions  are  unique. 

All  of  our  results  hold  for  the  matched  microstfucture  noodeLprobleni.  Specif- 
ica,lly,  Lemma  4  and  Corollary  4  are  obtzuned  from  Proposition  3,  and  Proposi¬ 
tion  4  is  actually  simpler  for  the  matclied  problem.  The  analogs  of  Theorem  1  and 
Corollary  5  show  that  the  matched  problem  (l.l.a),  (l.l.b),  (l.l.c')  h^  a  unique 
generalized  solution  (a, -6)  €  C(Q,T;ii^(fi)  x  L^(Q)]. 

Themext  theorem  shows  that  if  the  data  is  further  restricted,  the  generalized 
solutions  will  satisfy  the  partial  differential  equation  (1.1).  The  following  notation 
is  used, 

r(T)  =  L%T-,r{ayxr{Q)]  i<r<oo, 

V  =  w^’p{a)  xWg , 

V{T)  =  LP  [0,T;  Wi’^(n)]  x  , 

V(r)  =  [0,T;  x  wj]  . 

Theorem  2.  Assume  the  hypotheses  of  Theorem  1  and  in  addition  that  {f,F)  € 
L-(T)  n  V(T)'  and  (a(0),6(0))  G  D{C)  f)V'.  Then  the  generalized  solutions  of 
Corollary  5  satisfy 

(3.8.a)  (a,6)eV(T),  {u,U)eV{T), 

(3:8.b)  !(<.,  6)  +  (Aiu)  -  V/i,  B(,U)  +  7V)  ='(X  F)  in  V(f )'  , 

(3iS.c)  (a,  6)  G  (cld(«),  5^(1/))  ,  p  Gdm{Xu  —  '/U'). 

Proof.  The  results  of  Grange  and  Mingot  (14)  show  that  the  step  functions  {(F’.  b-^') 
and  generated  from  the  differencing  scheme  (3.7)  converge  weakly  in  V(T) 
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and  V('r),  respectively.  Mor6,  ver,  equation  (3.8)  will  be  satisfied  in  the  limit,  pro¬ 
vided!  the  weak  limits  (o,  6).and  (tt,  17)  satisfy  (a,  !))  G  (cl^(u),  5‘I>(C/)).  To  establish 
this  inclusion,  let  (UiV)  G  V(T)  and  (5,6)  G  (d^(v),di*(V)).  The  growth  condi¬ 
tions  on  ^  and  $  guarantee  that  (a^ ,  b^)  and  (a,b)  €  V(Ty  are  functions,  so  it  is 
possible  to  define  (a^  —  5,  6^  —  6),  to  be  the  pmr  of  functions  truncated  above  and 
below  by  ±s  (s  >  0).  This  pair  of  functions  is  bounded  in  L°°(T)iand  converges  in 
L^(r)-to  (a  —  5,  6  —  6)a,  so  converges  \n  I/i^T)  for  1  <  r  <  oo.  If  r  >  max(p',  5'),  it 
follows, that  L'^(T)  C  V(T)',  so  the  sequence  (a^  —  dy  b^  —  6),  converges  strongly 
in  V(T)'.  The  mOnotonicitj'  of  9^  and  9$  imply 

0  <  ((a^  -  a,  6^  -  b),  ,  (u^  -v,U^-  V)^  . 

Passing  to  the  liiriit  as  N  —*  00,  and  then  letting  s  r-+  00  yields 

6  <  ((a  -  a,  b  -  b)  ,  („  -  V,  U  -  Vfj  ,  (5,  b)  G  (9^(u),  9$(F))  . 

Since  (9^(*),  9$(^))  is  maximally  monotone, it  follows  that  (a,  6)  G  (5^(u),  9$(C7)).  I 

Finally,  we  note  that  the  corresponding  solution  of  the  matched  problem  sat¬ 
isfies 

(3.8.a')  (q,,6)  G  V(r)  ,,  {lyU,  ff)-G  V(T)  , 

(3.8.b')'  ^^(a,b)-+{A{-tU)MU))  =(.f,F)  in  V,(T)' , 

(3.8.c')  (a,  b)  €  (av’C-rC^).  .  U€Wo, 

where  the  space  Vo(r)  is  givenhy 

Vo(T)  =  [U  G  T^.(0,r;  Wi] :  7(i^)  G  T'’[0,T;  TFj-'’(a)] } 
with  the  appropriate  norm  for  which  €  V(T)  for  each  U  G  Vo(T). 
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The  work- described  below^wM  done  with  N.Jv  WjE^ngton  of  Garnegier 
Mellon  University.  We  ^tablish  that  the  Cauchy-Dirichlet  problem  for 
degenerate  parabolic  systems^of' the  fonri 

(T.a)  ^a(«)  -V-  A(x,  Vti)  + /  B{x,s,V^U)-Vds  ^  f  ,  *  , 

Jtm 

(l.b)  |5(K).-V,.-B(ii!,,V,B)3J’,  y  m., 

(i.c)  9  0  }  y.  €,r,. , 

is  well-posed,  Here-ft  is  a  domain  in  R”  and  for  each  value  of  the  macro- 
variable  X  €  «  is  specified  a  doinain  Q,  with  boundary  F*  for  the  micro- 
variable  y  €  ilx’  E<^h  of  a,  6,  ft  is  a  inaximal  monotone  graph.  These 
^aphs  are  not  nec^sarily  strictly  inerting;  they  may  be  piecewise  con¬ 
stant  or  multi-valued.  The  elHptic  operators  in  (i.a)  and  (l.b)  are  non¬ 
linear  in  the  gradient  of  degree  p  —  1  >  0  ^d  y  —  1  >  0,  respectively, 
with  ^  -b  ^  so  some  specific  degeneracy  is  also  permitted  here.  Cer¬ 
tain  first  order  spatial  derivatives  can  be  added  to  (I.a)  and  (l.b)  with  no 
difficulty,  and  corresponding  problems  with  constraints,  t.c.,  variational 
inequ^ities,  caii  be  treated  similarly.  A  particular  example  important  for 

*  This  work  was  supported  by  grants  from  the  National  Science  Foundation  and  the 
OflSce  of  N aval  Research . 
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applications  is  the  linear  constraint 

(lx')  U(xyy,i)  =  w(x,#)  ,  y  e  Fi  ,  T  €  fi 

which  thenvreplaces;(l.c).  The;system;(*i);with  /i(.s)  =  is  called 

a  regulaTizei  micro -structure  models  and  (l.a),  (l.b),  (l.c')  is  the  cone- 
sponding  motcAcd  micro-structure  model  in  which  (formally)  e  — ♦  0.  An 
example  of  such  a  system  is.  a  model  for  the  flow  of  a  fluid  (liquid  or  gas) 
through  a  fractured  medium.  In  such  a  context,  (l;a)  prescribes  the  flow 
on  the  global  scale  of  the  fissure  system  and  (l.b)  gives  the  flow  on  the 
microscale  of  the  individual  ceU  at  a  Specific  point^z  in  the  &sure  systOT. 
The  transfer  of  fluid-  betwera  the  cells  and  surrounding  mediiun  is-pre- 
scribed  by  (l.c)  or  ^l.c').  A  major  objective  is  -to  accuratdy  modebthis 
fluid  exchsmge  betwc^  the  cellSzand  fissures. 

Systems  of  the  fom  (1)  were  developed^in  [10],  [11],  [3]  in  phya^ 
chemistry  as  models  for  difliision  throu^  a  mediuin  with  a  prescril^ 
microstructure.  Similar  systems  arose  in  soil  science  [2],  [6]  and  in  reservoir 
models  for  fractured-  media  [4],  [7].  By  homogenization  methods  such 
systems  are  obtained  as  limits  of  exact  micro-scale  models,.^d  then  the 
effective  coefficients  are  computed  e^licitly  froin  local  material  properties 
[13],  [8],  [1].  An  exiStence-tiniqueness  theory  for  these  linear  probl^s 
which  exploits  the  Strong  parabolic  structure  of  the  system  was  given  in 
[12].  One  can  alternatively  eliminate  and  obtain  a  sin^e  functional 
differentia  equation  for  u  in  the  simpler  space  X*(ft),  but  the  structure  of 
the  equation  then  obstructs  the  optimal  parabolic  type  results  [9].  Also 
see  [5]  for  a  nonlinear  system  with  .reaction-diffusion  local  effects. 

We  begin  by  stating  and  resolving  the  stationary  forms  of  our  systems. 
Denote  by  7*  :  W’^’’(Dx)  -*  Xr*(r*)  the  trace  map  which  assigns  boundary 
values.  Let  T*  be  the  range  of  7*;  this  is  a  Banach  space  with  the  norm 
induced  by  7*  from  Since T*  is  smooth,  there  is  a  unit  outward 
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normal  at  each  .‘f  €  Tr-  We  shall  denote  by  Vy  the  gradient  on 

'Trb?(ni)  and  hy.  V  the  gradient  on  Tr*'^(n). 

In  order  to  prescribe  a  measurable ;fcimily  of  cells,  €  fi},  let 

Q  C  fi  X  FL”  be  a  given  measurable  set,  and  set  ft*  =  {y  GTC*  :  (x,  jr)  € 
-Q}.  Each;  f2x  is  mes^urable  in  R"  and  by  zero-extension  we  identify 
»  L^(ft  X  R”)  and  each  L’(f2x)  ►  Jj’(R’*).  Thus  we  obtain  from 

above 

L^(Q)  ^  {ff  €  L’(ft,i«(R"))  :  U{x)  €  ,  a.e.  x  €  q}  . 

We  shall  denote  the  duality  on  this  Banach  space  by 

W  *)l,(Q)  =  </i ,  U  6  ,  «  6  e'lQ)  . 

The  state  space  for  our  problems  will  be  the  product  (ft)  x  (Q). 

Note  that  W^'*(ft*)  is  continuously  imb^ded  in  L*(ft*),  tiniformly 
for  x:€  SI.  It  follows  that  the  direct  sum 

W,  s  L»(ft,  W^>»(ftx))  =  {If  €  L*(Q) :  U(x)  €  W’*>»(ft.)  ,  a.e.  x  €  ft  , 

and  /  ||tf(x)||^i,,  dx  <  oo}. 

Jn 

is  a  Banach  space.  We  use  a  variety  ^  such  spac^  which  cw  be  con¬ 
structed  in  this  manner.  The  Banach  space  Wj ’*'(ft)  x  will  be  the  en- 
eryy  space  for  the  regidarized  problem  (^1).  We  ^uine  the  trace  maps  7*  : 
W*’*(fti)  L^(rx)  sxe  uniformly  bounded.  Thus  for  each  If  €  VV,  it  fol¬ 
lows  that  the  distributed  trace  7(1/)  defined  by  7(lf)(x,s)  =  (7r(tf(x))(s), 
5  €  r*,  X  6  ft,  belongs  to  L*(ft,X*(r*)).  The  distribute  trace  7  maps 
W,  onto  r,  =  X«(ft,T*)  ^  X*(ftiX«(Fi)). 

Denote  by  A  the  map  of  X*(ft)  into  7^  given  by  At;(x)  =  t;(x)  •  T*,  a.e. 
X  €  ft,  V  €  X*(ft);  A  is  an  isomorphism  of  X*(ft)  onto  a  clo^d  subspae  of 
Tg.  Tne  dual  map  A'  taking  7^^  into  X’'(ft)  is  given  by 

yg(v)  =^(Au)  =  f  ^*(1*)  ■  v{x)dx  ,  g  e  Tg  ,  u  €  X’(ft)  , 


o 


so  we  have  .A?fif(x)  =  a.e.  x  €  fi.. 

We  consider  elliptic  differential  operators  iri.divergehce  form  as  real- 
izations  of  monotone  operators  from  Banach  spaces  to  their  duals.  The 
global  diffusion  operator  — is  given-by 

>tu(u)  =  f  A{x.yVu{x))Vv{x)dx  ,  €  Wp’^(n)  . 

Jn 

Similarly' define  for  each  x  €  ft 

Bxtv{v)=  I  S(i,y,V,ii>(y))V,u(y)dy  ,  w,v  . 

Ja, 

Also,  we  shall  denote  by  B  :  Wq  VV'  the  dwtn’&nted;  operator  con¬ 
structed  from  the  collection  {B*  :  x  €  ft}  by 

BU{x)  =  Bx{Uix))  ,  a;e.  x  €  ft  ,  U  €  , 

and  we  note  that  this  is  equivalent  to 

BU{V)f^  /  Bx(U{x))V{x)dx  y  U,V€  Wj  .. 

Jit 

The  cpupling  term  in  omrsystem  wU  be  given  as  a  monotone  graph 
which  is  a  subgrrwiient  operator.  Thus,  assume  m  :  R  — »  R"^  is  convex, 
continuous  and  bounded  by  G(ls|*  -j- 1),.  Then  by 

•Tn((3)  =  /  /  m{g{x,  s))  ds  dx  ,  g  €  I’  (F* )) 

JotJt, 

we  obtain  the  convex,  continupusim  :  i^(ft,L*(Fi))  R"^., 

For  the  case  of  the  degenerate  stationary  problem  correspondingrto' 
(1),  we  obtain  theTollowing  result. 

Let  tpr:  R  — >  R"^  and  $  :  R  — ♦  R"*"  be  cpnyex  and  cdntinuousj  with 
V?(0)  =  $(0)  =  0,  and  assume 

+  ,  $(s)  <  C(|s|«  +  1)  ,  s=€R  . 
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‘For  each  pair  /  G  (fi),  F  G  W', >there  exists  a  solution  of 

(2.a)  u  G  Wo^ ’'’(«) :  q  +  A{u)  -  A'(//)  =  /  in 
(2.b)  V  G  W,  :  6  +  B{U.)  +  =  F  in  W''. 

(2.c)  ’fi  ^  dmi.'yU -r  Xu)  in  L®  (fi,L’,(F.t))  j  and 

(2.d)  a  G  av?(ti)  in  b  em^U)  in'  i«'(Q)  . 

For  any  such  sdlutioh^we  have 

I  b{x,y)dy  +  I  /i(x,s)d5  ==  (,F(®),1*)  a.e.  x  G  ft  . 

JClg 

Next  we  restrict  our  system  to  appropriate  products  of  L’’  spaces. 
The  Hilbert  space  case,  r  =  2,  serves  hot  only  as  a  convenient  starting 
point  blit  leads  fo^he;  generalized  accretive  estimates  we  sh^  n^d  for  the 
singular  case  of:(l)  imwhich  a  or  6  is  not  only  nonlinear-but  miilti-valued. 
The.statiohaxy  operator  for  (1)  is  shown  to  be  m-accretive  in  the  space, 
S.0  we  obtain  n  gtnerdlizei  soluiion  in  the  sense  of  the  npnlineM  serhigroup 
theory  for  general  Banaxh  spaces.  As  an  intermediate  step  we  show  the 
special  case  of  d  =  b  —  identity  is  resolved  as  a  strong  solution  in  every 
i*"  space,  1  <  r  <  oo,  and  also  in  appropriate  dual  Sobolev  spaces. 
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NODAL  SUPERCONVERGENCE  AND  SOLUTION  ENHANCEMENT  FOR  A 
CLASS  OF  FINITE-ELEMENT  AND  FINITE-DIFFERENCE  METHODS* 

R.  J.  MACKINNONt  ANO  G.  F.  CAREYt 

Abstract.  A  class  of  finite-element  methods  for  elliptic  problems  is  shown  to  exhibit  nodal  superconver¬ 
gence  in  the  approximate  solution,  and  some  equivalence  properties  to  familiar  fiiiite-diflerence  operators 
are  demonstrated.  The  supe'rconvergence  property  is  expIoited:in  a  Taylor  scries  analysis  to  demonstrate 
Gauss-point  superconvergence  for  the  derivatives  of  the  approximation.  A  post-processing  formula  for  the 
derivative  at  the  nodes  is  constructed  and  shown  to  exhibit  superconvergence.  The  nodal  superconyergence 
property  can  be  exploited- recursively  to  further  enhance  the  finite-element  or-finite-difference  solution. 
Supporting  numerical  studies  are  given. 

Key  words,  finite  element,  finite  difference,  superconvcrgence,  post-processing 

AMSfMOS)  subject  classificatiohs.  6S-L60,  6S-N30 

1.  Introduction.  In  this  note  we  consider  a  Galerkin  finite-element  approximation 
of  the  Dirichlet  problem  for  the  equation  Lu  =fm  CL.  Here  0  is  a  union  of  rectangular 
subdomains,  L  is  a  secgrid-order  elliptic  differential  operator  with  smooth  coefficients, 
arid  «  is;assumed  to  be  sufficiently  smooth;  By  introducing  an  appropriate  integration 
rule  for  element  quadrature  we  show  that  therGalerkin  approximation  u^,  defined  on 
a  square  mesh  of  piecewise  bilinear  elements^  is  equivalent  to  a  familiar  finite-difference 
approxiination  of  m.  Discrete  uniform  error  estimates  for  this  difference  approximation 
are  known  (Bramble  arid  Hubbard  [1]).  These  estimates  iiriply  that  difference  quotients 
of  the  error  have  the  same  order  of  convergence  as  the  error  itself;  i.e.,^0(/l^)  for  the 
bilinear  element.  It  follows  that  th[s  nodalrsupercqnvergerice  property  holds  for  the 
standard  Galerkin  approximation  with  higher-order  (or  full)  integration.  We  use  this 
result  to  prove  new  superconvergence  results  arid  show  how  simple  and  accurate 
supercoriVergent  postTprocessing  formulas  for  the  solution  and  derivatives  can  be 
derived  using  Taylor  series  expansions.  Although  the  fdrmulatipn  and*  analysis  pre¬ 
sented  here  is  for  problems  in  two  dimerisibns,  the  results  apply,  to  problems  in  one 
dimension  as  well,  and  extend  directly  to  three  dimensions. 

2,  Formulation  and  analysis. 

2.1.  Nodal  solution  superconvergence.  Consider  the  boundary  value  problem 

in  the  unit  square  fl  =  (0,  l)x(0,  l)"with  Dirichlet  data 

(2)  u  =  g  on  dCl. 

Here  we  assuirie  that.a,'b,  c,  and  /  are  smooth,  and  L  is  uniformly  elliptic  in  fl. 

The  Galerkin  finite-element  approximation  to  (1)  is  defined  to  be  m,,  ^  satisfying 
the  essential  boundary,  condition,  and  such-that 

(3)  B{ui„Vh)  =  (:f,.Vh) 
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for  all  VhS  Ho^  where  (•,•)  is  the  L^(fl)  inner  product  and  is  the 

bilinear  functional 

(4)  B(w,  !/>)=  [a,w^i/>^  +  a2Wyif>y-i-(b,w^  +  b2Wy-i^cw)i;/]dxdy. 

Jn 

Now  consider  a  uniform,  partition  of  ft  into  square  elements  of  size  /i  and  take 
H-(ft)  to  be  spanhedf by  C“  piecewise-bilinear'functions  defined  on  this  partition. 
Approximating  the  integrals  in  (3)  by  a*  suitable  integration  rule  applied  over  each 
element,  we  get  the  approximation  «/,)  =  (/ U/,) a  for  all  u;,  s  The  resulting 

algebraic  system  is 

(5)  B*u,  =  ffc 

where  the  precise  forms  of  B^  and  depend  on  the  particular  integration  rule  used. 

For  clarity  of  exposition,  let  us  first  consider  the  case  where  coefficients  a,  b,  and 
c  are  constants,  and  at(2  x  2)  trapezoidalUntegration  rule  is  used  to  evaluate  integrals 
in  (3).  Accordingly,  evaluating  the  coefficients  in  (5),  for  typical  interior  node  point  i 
at  (.Xi,y,)  whii  test  function  0*1,  we  obtain 

Bh{Uh,VM)  =  -{ai[Uhix,+h,yt)-2Uh(x,,  y,)  +  «*(x,  -  h,  y,)] 

+  Oiluhixt,  yi  +  h)-2uh(.Xi,yi)  +  Uhix,,  y.  -  /i)]}: 

+  h  Y  [wa(x,+:A,  y,)  -  u^ix,  -  h,  y,)] 

+  h  j[Uhix,,  yj+h)-  u^{x,,  y,  -/))]  +  cu^ix,,  yi)h^ 

(f,.VM)h^fixi,yi)hl 

For-this  case;. we ;see  from  (6);and  (7)-tfiat  (5)  is  equivaleht^to  the  five-point  central 
difference  approximation  to  ( 1 ). 

Bramble.:and  Hubbard  [1]  Kave  shown  that,  for  a  solution  u  of  (l);having  bounded 
fifth  derivatives,  the  gridpoint  error  e,  -  u(x,,  y,)  -  Uh{x,,  y,)  for  the  five-point  difference 
approximation  satisfies 

(8)  e,  =  4){x,,  y,)h^+  B(x,,  y,,  h)h^ 

where  <(>  has  Lipschitz  continuous  second  derivatives,  and  R  is  uniformly  bounded  in 
x^yi,  and  A.  It  follows  from  (8)  that  the  solution  to  this  finite-element  problem  (5) 
has  gridpoint  errors  oForder  0{h^).  It  should  be  emphasized  that  this  estimate  is  a 
gridpoint  result  for  the  discrete  problem,  and  is  of  the  same  order  as  the  global 
estimate  usually  encountered  in  finite-element  theory. 

Remark.  If  the  domain  discretization  error  is  zero  (as  assumed  here)j,.  then  it 
follows  directly  from  the  proof  in  Bramble  and  Hubbard  that  max,,. |B(X|,  y,,  h)\  g  Ch, 
constant  C,  sb  the  final  term  in  (8)  is  actually  0{h*). 

In  (8)  0.is  the  solution  to  the  auxiliary  problem 

L4>  =  t(u)  in  ft, 

(9)  '  ’ 

0=0  on  5ft 

with  T  the  truncation  error.  In  particular, 

(10)  Bh(u,  VH,)  =  h^Lu(Xi,y,)  +  h*T,{u)  +  0(h^) 


(6) 


and 

(7) 
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where  T((u)  denotes  t(«)  evaluated  at  interior  gridpdint  (x.-.y,).  For  trapezoidal 
integratioii  t(«)  corresponding  to  (6)  is 

(11) .  t(m)  =  [aiU„xx+a2Mwy”2(6,u«,+ 

According  to  (10),  the  discrete  approximation  given  in  (6)  and  (7)  has  local  truncation 
errors  of  order  ~Oih*).  On  dividing  by  h^,  we  see  that  the  differential-  operator  is 
approximated  to  a  local  accuracy  of  order  0(/i^).  Even  if  a  more  accurate  quadrature 
scheme  is  used  for  integrating  (3),  0{h*)  truncation  errors  remain.  Their  precise  forms 
depend  on  the  integration  rule  used.  It  follows  that  the  estimate  in  (8)  gives  the  best 
possible  rate  for  the  nodal  solution  error  irrespective  Of  the  increase  in  quadrature 
accuracy; 

This  conclusion  also  holds  for  the  case  of  smooth  variable  coefficients^  since  their 
variations  only  introduce  0(/j'*)  truncation  errors.  (See  the  Appendix  for  an  example.) 

2.2.  Derivative  superconvergence  points.  Consider  first  the  problem  of  deriyative 
calculation  from  the  bilinear  finite-element  nodal  interpolant  u,  of  u.  Simply  differen¬ 
tiating  the  expansion  on  element  D,,  we  have 

4 

(12)  M,x(x,  y)  =  I  y) 

Jml 


where  u^.are  the  interpolated  nodal  values  for  are  the  element  basis  functions, 
and  X,  y  is  an  arbitrary  point  in  the  element. 

Next,  we  introduce  Tay!or  serie!..expansions-for  M;  =  u(x^-,yy)  about  x  =  jc,  y  =  y 
with  5j=  X;-x,  SJ =.yj-y  to  obtain 

Uj  =  u(xj,yj)=:u{x,p)  +  uiixJ)SfA-Uy<^xJ)S^ 

-b.«xx(x,y)(«;)V2!-+M^,y)(5j’)V2!H-Mx.(x,y)s;5/^  ■ . 


Using  (13)  in  the  right  side  of  (12)  and  regrouping  terms,  we'havc' 

«x=  S  ^JxiSj?}uxx 

"  [^i. 

where  for  notational  convenience  u  =  u(x,  y),  =  Ujc(x,  y),  tfex  =  </';x(^,  y).  and  so  on. 

A  similar  expression  holds- for  Uy. 

Now  the  derivative  of  the  approximate  solution  Uh,at  x,  y  in  D,  is 


(15) 


W/ix-  Z  ^hj^jx- 


Subtracting  (15)  from  (14)  yields  the  error  In  the  derivative 


(16) 


:=  I  e/yx-  Z  fe(5y)^  Wx, 
y-i  L2!y-i  J 


"[^i-  '^x6;6j]:f3x,+o(/.^). 

When  we  introduce  nodal  estimate  (8)  for  ej  and  use  the  fact  that  derivatives  and 
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hence  difference  quotients  are  bounded,  the;first  term  on  the  right  in  (16)  satisfies 

Z  ^jxCj  =0(h^). 

j-t 

This  implies  that_(l6)  will  be  an  0(A-)'approximation,  provided  the  femainingifirst- 
order  tehns  are  zero  or  coliectively  cancel.  On  examinatipn^  we  find  that  coefficients 
of  Uxy  and'Uyy  are  zero  for  all  x,  j>  in  but  the  coefficient  of  u,*  is  zero  for  all  /  with 
jc  =  (xi  4-X2)/2.  Therefore,  U/ix  is  supercpnvergent  on  the  //he  bisecting  the  horizontal 
sides  of^Oei  Similarly,  u^y  is  Qih^)  along  the  line-bisecting  vertical  sides  of  ft,.  Hence, 
the  centroid  (Gauss  point)  is  the  superconvergent  point  for  both  u**  and  Uhy. 

2.3.  Nodal  derivative  extraction.  Now  consider  the  calculation  of  derivatives  (flux 
components,  stresses)  at  interior  node  . point  x,,  y,.  For  a  solution  m  of  (1)  having 
sufficient  smoothness  in  the  interior  of  ft.  Bramble  and  Hubbard  [1]  prove  the  following 
estimate  for  the  equivalent  finite-difference  scheme: 

(17)  \Dle(x,,y,)\^c„[\eU,  +  0{h^)] 

where  is  ah  nth  order  difference  quotient  having  O(h-)  truncation  error,  c„,'is  a 
constant  independent  of  /i,  and  le|n„  =  max;f,_y,  |e(x;,y,)|.  For  the  problem  considered 
here  we  have,  according  to  (8), 


Thus,  (17)  becomes 


|e|n*  =  0(/i^). 


\Dle{:X„y,)\^eW. 


This  result  can  npw  be  used^to  derive  a  supercbrfvergentiapproximatipn  for.  the^flux 
components  a,  Ui,-.a2Uy(and'>hehce  derivatives  if  desired)  atinode  point  x,,yi. 
A.Taylor  series  expansipn.fof  u{xi-TH,y,)  about  (X(,y()  yields 

Q  ‘h 

(19)  a,  u,(x,,.;',)  = [m(x,,  y,)  -  u(x,  -  h,  y,)] +-  a,  Uy^{Xi,  y,)  +  O(h^). 

Replacing  a, u,„c(x,,  y,)  in  (19)  using  differential  equation  (1); and  then  introducing;the 
following  difference  formulasTor  a,x,  u^,  u,.,  and  (a^Wy)^ 

au(xr,yi)  =  — - p - ^+0(/i), 

,  ^  u{x,,y,)-u{;x,-h,y,)  ^ 

Ux(:Xi,y,)  =  - - +  Oih), 

^  uix,,y,+h)-uix,,y,-h)  ,  ^,,2^ 

Uyixi^yi)  =-^- - ^ - -r~+  0{.h% 

,  /  ia2ixi,y,  +  h)  +  a2(,x,,yi))\.  ,  .  , 

{a2Uy(x,,y,))y  =  (u(x,-,  y, + /i)  -  u{x„  y,)] 


ia2ix,,y,)  +  a2{x,,y,-h))^  ,  ,  , 

- rn - Mxr,  yi) -uiXi,  y,  -/))]  +  0(/)^), 
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we  obtain 


(ai{x,,y,)  +  ai{x,-H,y,)) 
2h 


[u{x,,y,)-u{x,-h,yi)] 


(20) 


+ 


{:a2(.xi,yt+'h)+a2(x,,yi)) 

'4/1 


[w(x,,  y,  +  ft) -m(x, ,>>,)] 


ia2{xi,y,)  +  a2{x,,yi-h)) 
Ah 


[U{x,,y,)-u{x,,yi-h)]i 


+^[h(x,,  y,)  -  u{x,  ^  ft,  y,)] + Y  [«(x,,  y,  +  ft)  -  u(x,,  y,  -  ft)] 
2  4 

+  c|u(x,,  yd~^fixi,  y,)  +  d(ft^). 


Note  that  (20)  is  ah  O(ft^)  difference  formula  involving  nodal  values  of  the  exact 
solution  u.  Oh  introducing  the  finite-element  approximation  for  u  on  the  right  in 
(20)i.we  define  the  approximation  for  aiUx(xi,y,) 


aiUx{x,,y,) 


iai{xi,yi).+  ai{x,-h,y,)): 

'2ft 


[Uh{Xi,yi)-Uh\Xi-h,yi)] 


(21) 


(a2(j^(.yi+ft)-+a2(yi.y()) 

4ft 

(a2(Jfi.yf)tQ2(x„y,,rft)) 

4ft 


[Mh(x,,y,  +  ft)-u;,(x,,y()]= 


[wa(^C/,  y/)  -  UH(Xi,  y,  -  ft)]; 


+Y  [“*(*'»  3’/)  ^  Uhixi  •r  ft,  y/)]H-^  [Hfc(jc,,  y,+-ft)  -  u^(3c,,  y,^  ft)] 

+jUH'^x,,y,)^^fix,,y,). 


Subtracting  (21)  from=(20)  and  using  (18),  we  find  that  (21)  is  a  superconvergent  C*'  ■) 
flux  apprpxi|[iation.:(For  a  related  study  of  derivative  approximations  see  MacKinnon 
and  Carey;[5].) 

Finally, ilet  us  use  this  result  to  analyze  a  finite-element  projection  technique  for 
flux  post-processing.  This  technique 'is  based  on  the  integration-by-parts  procedure  in 
the  finite-element  integral  statement,  from  which  we  define  the  projcction.relatiohship 
for  a,«J: 

(22)  a,M*Uwd5=  (ajM*xt'wx+a2«A;-t'My+(b'Vufc-+CM*-/)uJdxd'y 

Js,  Jn, 

where  si  is  defined  by  element  sides  connecting  gridpoints  (x„  yi  -  ft),  (x,-,  y/),  (X|,  y* + ft ) 
and  ft,  represents  the  two-elemeht  patch  defined  by  gridpoints  (x,-  ft,  y,j,  (X{,y,), 
(xj,  yi + ft ),  (Xj  -  ft,  y/ + ft),  (X(,  yi  -  ft),  (x,  ~h,yi-  ft).  This  approach  has  been  examined 
in  one  dimension  by  Wheeler  [6],  Dupont  [4],  and  Carey  [2].  In  two-dimensional 
numerical  teM  cases  the  method  has  been  demonstrated  to  yield  an  0(ft*)  approximation 
to  the  nodaf  flux  when  n,Ux.is  assumed  to  be  piecewise-cohstant  over  Sj  (Carey,  Chow, 
and!Seager:[3]).  Indeed,  iFwe  integrate  (21)  usihg  the  trapezoidal  rule,  as  described 
in  the  Appendix,  then  the  resulting  discrete  formula  for  a,uj  is  identical  to  (21).  This 
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then  confirms  the  observed  numerical-convergence  rate  of  If  a  more  accurate 

quadrature  scheme  is  used  to  evaluate  (22),  the  resulting  difference  formula  . is,  in 
general,  different  from  (21).  However,  a  simple  Taylor  series  analysis  confirms  that 
the  formula  is  0(/i^)  accurate. 

3.  Nodal  solution  enhancement.  In  this  section  we  apply  the  results  obtained  in 
the  foregoing  analysis  and  formulate  a  new  scheme  to  compute  an  accurate  approxima¬ 
tion  e*  tp  the  gridpoint  error  e,  =  u(x,,  y,)  -  Uh(xi,  yi).  This  approxiniation  may/^then 
be  used  to  improve  u*  and,  moreover,  increase  the  asymptotic  rate  of  convergence  of 
u*  and  its  derivatives.  We  point  out  that  although  the  formulation  presented  here  is 
for  problems  in  two-dimensions,  it  includes  the  one-dimensional  case  by  simply  setting 
y  derivatives  equal  to  zero. 

In  the  interest  of  clarity,  we  restrict  our  analysis  to  the  constant  coefficient  case 
described  by  (6)-^(  11),  The  extension  to  other  cases  involving  different  quadrature 
schemes  and  variable  coefficients  is  straightforward  in  view  of  our  previous  results. 
First  recall  (10)  and  (11): 

(10)  Bfc(u,  Vm)  =  h^Lu(x,, yi)  +  h*T,(u)  +  0{h^), 

(11)  T(u)  =  -n[a  +  OiUyyyy-Kbi 

^XXX  +  biUyyy)y. 

Now  from  estjmate  (10,  since  is  0(/j^)  accurate  at  the  node  points,  any  nth-order 
difference  quotient  D)|?of  U;,  also  converges  to  the  exact  value  D"u  ara  rate  of  P(/i^). 
^erefqre,  at  any  interior  node  point./,  t(«)  cah;be  rewritten  using  difference  quotients 
DhUh  as 

n(«)  =  -  i^[aiDtx/M/.  +  a2PvW/.-2(^DL/MA  +  fc2Dj,^(Mj]-t-0(/i^) 

^  ^rM  +  Oih% 

(Note  that  since  the  fifth  derivatives  of  w- are  assumed  bounded,  then  t,(m)  atmode 
.points  on  bdundary.afl  can  also  be  approximated  to  0{h^)  accuracyiby  simply  using 
ah  0(/i^)  extrapolation  to  the  boundary.) 

Next,  interpolate  the  nodal  values  T/(M)-ih  the  piecewise-bilinear  basis  as 

(24)  t(u)=  I  Tj{u)^jix,y)  +  0{h^) 

j-i 

where  N  is  the  number  of  node  points. 

Introducing  (23)  in  (24),  we  have 

(25)  r(u)=Z  -ri.jMi/'jix,  y)  +  O(b-). 

Replacing  7-,(u)  in  (10)  using  (25),  we  have 

(26)  B^(u,  Vhi)  =  h^Lu(Xi, >'()+./i'‘t,,,(wJ-f  0(/)*). 

Using  (26)  in  place  of  (10),  the  estimate  (8)  now  has  the  form 

(27)  e,  =  (l>*{x,,  y!)h-+;R*{Xi,  y,,  b)h^ 
where  </>*  satisfies  the  auxiliary  problem 

i-4>*  =  fhiUh)=I.  Ujiuh)<pj{x,  y)  in  II, 
y-i 

<f>*  =  0  on. 512. 


(28) 
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The  objective  now  is  to.construct  a  finite-element  approximation  </>*  to  (;>*  in  (28), 
and  then  use  this  approximation  in.the  leading  term  On  the  right  of  (28)ito  obtain  an 
accurate  correction  to  the  nodal  solution. 

First,  let  us  assume  that  we  have  already  computediu^  from. (5)  using  LU  fac¬ 
torization  and  haversaved  the  computed  matrix  factors.  The  Galerkin  finite-element 
approximation  to  (28)  is  asTollows.  Find  <^*6  //''  such  that 

(29)  B{<f>t,Vh)  =  {rH{uH),Vh) 

for  all  u*e  H*c//o*  As  before,  if  we  evaluate  integrals  in  (29)  using  the  (2x2) 
trapezbidaf  integration  rule  we  get 

(30)  = 

where  t/,i(uj,)  is  defined  in  (23)  and  i)/,()  is  analogous  to, (6). 

Since  the  matrix  factorization  of  B/.  is  already^given  from  the  previous  calculation 
of  U/,,  the  approximate  function  in  (30)  can  be  computed  efficiently  once  T;.i(M,,) 
are  computed  (for  the  Dirichlet  problem  t,.i(u*) -is  needed  at. interior  points  only); 
t/,i(«/,)Ts  easily  computed  using  ohe-dimetisional  difference  formulas.  For  example, 
we  may  write 

(c+l  W" 

D"fU=  I  ^%(^)«(§)  +  0(/i''"-"). 

(31)  t  . 

=  DZfU  +  p(/.'‘-"'-"),  k^n.  S  =  x,y 

where  u(fy)^are  node  point  values  of  function  m;  and  are  Lagrange  polynomial 
shape  functions  oFdegree-fc  In  particular,:for  a  second-order  (0.(/i^))  approximation 
to  d%/dx*,  n-4  and  k  =  l  Noteahat  (31)  canrbe  usedUo  approximate  d‘*«/dx‘‘  at 
interior  nodes  near  the  boundary.  For  this  case  (31)  is  simply  arone-sided  difference 
formula  involving  interior  node  point  values  of  li- only. 

Solution  <^jf  from  (30)  will  approximate  4>*  with  accuracy  0(/i'’)  at  all-  node 
points,  where  p  depends  on  the  smoothness  of  solution  h  to  (1).  Note-that  t  in  (24) 
is  C*  in  view  of  the  assumptions  on  m.  Moreover,  f,,  in  (28)  is  C”  by  construction  so 
<f)*s:C^  and'p^  1.  Replacing  in  (28)  with  ^*,  we  have 

,  ,  ef  =  cf>t(x,,ydh^+R*{x,,y,,h)h^+OW-^'’) 

(32) 

=  4>U^+OXh^^'’) 

since  R*  is  0{h). 

^is  importaritiresult  implies  that  we  can  compute  node  point  errors  e*  haying 
at  lept  0{h^)  accuracy,  and  ©(/i"^)  accuracy  (p  =  2)  for  sufficiently  smooth  solution. 
An  immediate  consequence  of  this  result  is  that  we  can  also  increase  the  accuracy  of 
our  approximation  Uh  (and  its  derivatives,  if  desired)  from  0(/j3)  to  at  least  0(/i^). 
That  is,  the  enhanced  gridpoint  value  obtained  by  adding  the  nodal  correction  ef 
becomes 

(33)  uf,!=UM+ef. 

The  solution  enhancement  procedure  may  be  summarized  as  follows: 

(1)  Solve  the  finite-element  problem  using  sparse  LLf  factorization  and 

save  matrix  factors. 

(2)  “Process”  approximationrufc  and  form  associated  vector  Xhiuh).  Then,  using 
matrix  factors  of  B*,  solve  auxiliary  problem  Bft<j)*  =  T,,(MA). 
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(3)  Compute  approximate  node  point  error  correction 

ef  =  y,)h^ 

and  hence  the  “enhanced”  solution 

uf  =  u„+-ef. 

3|1.  Numerical;exampl^.  Numerical  test  studies  have  been  made  to  demonstrate 
the  effectiveness  of  the  nodal,  enhancemerit  post-processing  procedure.  (For  results 
related  to  the  application  of  post-processing  derivative  formula  (22)  and  related 
formulas,  we  refer  the  reader  to  Carey  [2],  Carey,  Chow,  and  Seager  [3],  and  MacKin- 
mon  and  Carey  [5].) 

In  the  first  test  case  we  consider  the  two-point  boundary  value  problem 

-:u^+u^  +  u=f,  0<x<\, 

(34) 

w(0)=.m(1)  =  0 

where /  is  constructed  such  that  the  analytic  solution  is;M  =  x(l  -x)(l  +  x)*. 

We  take  a  sequence  of  uniform  mesh  refinements  with  /i  =  5,  -re,  and  55,  Numeri¬ 
cal  integration  is'performed  using  the  trapezoidal  rule,=and  derivatives  u^xx  and  Uxxxx 
in  f  are  approximated  to  order  0(/i^)  and  P(/i^),  respectively,  by  six-point  difference 
formulas.  A  six-point  formula  for  Uxxx  was  used  because  it  is  computationally  convenient 
to  simply  differentiate  this  formula  and  use^the  result  to  approximate  Uxxxx- 

Node  point  errors  £,,  Ef  for  approximations  and  enhancement  uf  are  presented 
in  Table  1.  Note  the  substantial  increase  in  accuracy  and  asymptotic  rates  of  conver- 
gence  afforded  by  the  enhancement  procedure. 

Next  we  examine  three  approximatiohf  to  m*  at  x=  1.  These  approximatioris.are: 
the  standard  0(h)  derivative -approximation  u^x;  the  post-processed  derivative  u* 
given-by  (21);  and'the  enhanced  derivatjvc'denoted  by.uj*  and  given  by  an  0(/i^) 
sixrpdirit  difference  formula  operating  on  enhanced  solution  u*.  Results  are  presented 
in  Table  2.  Approximations  u^  and  u**  ate  O(h^)  and  Oih*)  accurate  as  predicted. 


Table  1 

Node  point  errors  £(x,),  E*(x,)  for  the  case  b  =  l. 


H 

£(0.2) 

£•(0.2) 

£(0.6) 

£•(0.6) 

£(0.8) 

£•(0.8) 

V 

0.085004 

"0.016524 

0.161265 

0.200329 

_0.162612 

0.015906 

is 

0.021283 

0.204E-3 

0.050081 

0.724E-3 

.0.040634 

0.389E-3 

is 

0.005323 

0.27E-4 

QQR:  1 

0.012520 

.0.61  E-4 

0.010157 

0.46E-4 

ii 

0.001330 

0.18E-5 

0,0031301 

■tX-llliOM 

0.3IE-5 

-(?(/.•) 

~0(/.») 

-Oih*) 

~om 

~0{h*) 

~0(/<*) 

~0{h^) 

Table  2 

Derivatives  and  u**at  1.  The  exact  derivative  is  m,(1)  =  -3^0. 


h 

uj* 

—m 

-14.303480 

-34.933828 

-31.476970 

17.69652 

2.933828 

0.5230300 

A 

-22.034701 

-32.136437 

-31.982434 

9.965299 

0.7364370-- 

0.0175660 

is 

-26.716385 

-32.184295 

-31.999583 

5.283615 

0.1842950 

0.417E-3 

is 

-29.280084 

-32.046085 

-31.999975 

2.719916 

~0{h) 

0.0460850 

~0(h^) 

0.25E-4 

~0(h*) 
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As  a  two-dimensional  test  problem  we  take  the  example 

(w**+u^^)  =  {2-A2x^){y-i)  -A2y\x'^-x’)  in  ft  =  (0, 1)  x  (0, 1) 

with 

(35)  «  =  0  on  aft. 

The  analytic  solution  is  the  polynomia} 

(36)  M  =  (x^-x’)(y-/). 

Node  point  results  for  a  sequence  of  calculations  bn  uniformly  refined  meshes  of 
=  and  p  are  given  in  Tables  3-5.  Agajn,  the  observed  rates  of  convergence 
coifobofate  our  analysis. 

Conclusion.  By  a  suitable  choice  of  quadrature  rule  the  finite-elemejit  approxima¬ 
tion  for  a  two-dimensional  elliptic  problem  has  been  related  to  a  familiar  finite- 
difference  approximation.  Nodal  supercbnvergence  ofthe  solution  then  follows  from 
an  estimate  of  finite-difference  theory.  Moreover,  any  nth  ordef  differehce  approxima¬ 
tion  having  Taylor  series  truncation  error  of  P(/»^)  at  a  node  point  converges  to  the 
exact  value  at  a  rate  of=0(/i^).  Therefore,  accurate  derivative  extraction:formulas  can 
be  derived  directly  using  Taylor  series  ideas. 


Table  3 

Node  point  errors ‘E(Xi)i  E*{x,y  along  x  =  0.8. 


h 

£(0.8,'0.4) 

£*(0.8, 0.4) 

£(0.8, 0.6)' 

£*(0.8, 0,6) 

£(0.8, 0.8) 

£*(0.8, 0.8) 

i 

0.012138 

0.004074 

0.023837 

0.005127 

■0.032477 

0.006333 

IHI 

0.006824 

A 

HmAjXHJxB 

0.40E-^-4 

0.006087 

0.80E-4 

0.007625 

0.93E-4 

3^ 

0;2E-6 

0.680E-3 

0.3E-6 

0.4E-6 

-O(h^) 

=  ~0(/i*) 

~0(fc») 

~-0(/.«) 

~0(A?)  . 

~0(A*) 

~0(A*) 

^0(A*) 

Table  4 

Node  point  errors  £(*(),  E*ix,)  along  y  =  0.8. 


h 

£*(0.4, 0.8) 

£(0.6, 0.8) 

£’(0.6, 0.8) 

i 

0.007514 

0.003668 

0.017644 

.  0.003680 

0.028428 

0.005004 

A 

:  0.001982 

0.21E-4 

0.004496 

0.22E-4 

0.0O7233 

0.48E-4 

Mil  II 1  1  f 

0.2E-6 

0,502E-3 

0.4E-6 

0.807E-3 

0.5E-6 

~0[h^) 

~0(A") 

-0(A=) 

-0(A*) 

~0(A') 

~0(h*) 

Table  5 

Derivatives  and  (*;>')  =  (0.8, 0.8).  The  exact  derivative  is  u^(0.8, 0.8)  =  -0.1387215. 


h 

i 

0.282010 

-O.miST 

-0.155165 

0.420731 

0;267035 

0.016443 

A 

0.137307 

-0.208948 

-0.139683, 

0.276028 

0.070227 

0.961E-3 

is 

-0.0310355 

-0.146723 

-0.138726 

0.107686 

0.4E-5 

~0(A) 

~0{h') 

-0(A*) 
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We  emphasize  that  since  this  Taylor  series  approach  relies  only  on  elementary 
analysis  concepts,  it  is  straightforward  to  understand  and  implement.  Furthermore, 
although  this  is  not  taken  up  here,  the  method  can  be  easily  applied  to  higher-order 
elements  and  problems  in  three  dimensions.  Also,  derivatives  can  be  extracted  from 
finite-difference;splutions  ih’tlie  same  manner. 

Finally,  using  the  truncation  error  in  an  auxiliary^problem  the  nodal  superconver- 
gcnce  property=can  be  further  expjoited.to  enhance  the  gridpoint  solution  accuracy. 
These  results  are  of  practical  significance  in  solution  and  derivative  post-processing 
and  also  for  a  posteriori  error  analysis  in  conjunction  with  adaptive  refinement.  The 
adaptive  refinement  aspects  will  be  taken  up  in  future  studies. 

Appendix.  Trapezoidal  rule  and  variable  coefficients.  For  the  case  of  variable 
coefficients  and  trapezoidal’integration,  we  have  from  (3)  at  interior  gridpoint  i 

(Ifl) 


where 


B, 


,  V  \ai(xi+  h,yi)  +  atixs,yi)  ^  , 

,(«*,  t'/.f)  =  -| - ^ +  h,  yi)  -  Uh{Xi,  yi)] 


a\(Xi,yi)  +  at(Xi-;h,yi) 

2h 


[Uh(x,iy,)-u^{x,-h,y,)'] 


(1.2) 


j  a2(Xi,yi+Jt)+a2{Xi,yi) 

2h 


[Uh{x„yi  +  h)-Ukix,,y,)] 


a2{Xiiyi)+a2ix,,y,~h) 

2h 


[uh{x,,yi)-UH{x,,y, 


+  h^ix,iyi)[Uh(x,+Hiy,)--  Uh{xi-^h,y,)] 

+  h~ix,,yi)[uH(xi,  yi +h)-  Uh(x,,y,  -h)]+ c(x,,  yi)Uh(x,,y,)P 


and 

(1-3)  fu=fixi,y,)h\ 

For  a  smooth  function  H',.ahd  using  Taylor’s  theorem, 

(1.4)  B,iw,v,,)  =  h^Lw{x„  y,)  +  h\,{w)  +  0[h^) 

where 

(15)  =  Hxwx  +  a|x»V«x  +  OlxxWxx  +  d|xxx>*’.v 

+  a2Wyyyy  +  a2yW^yy+  QiyyWyy  +  OlyyyWy  -  2(b  f  tiWyyy)]. 

From  (1.4)  and  (1.5)  we  see  that  variable  coefficients  produce  additional  p(/i'‘) 
terms  in  r.  Hence,  the  accuracy  of  (1.1)  remains  0{h^)  and  the  results  demonstrated 
for  the  constaht;coefficient  case  extend  directly  to  the  case  of  variable  coefficients. 
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Superconyergence  error  estimates  are  established  for  a  class  of  semilinear^problems 
dehned  by  a  linear  elliptic  operator  with  a  nonlinear  forcing  term.  The  analysis  is 
for  rectangular  biquadratic  elements,  and  we  prove  supefconvergence  of  the  de¬ 
rivative  component  along  associated  lines  through  the  Gauss  points.  Derivative 
postprocessing  formula  and  formulas  for  integrals  are  also  considered  and  similar 
superconvergen^eestimatesproven. 


1.  INTRODUCTION 

In  recent  years,  much  effort  has  been  ft^used  on  superconvergehce  phe¬ 
nomena  in  the  numerical  solution  of  differential  equations.  In  the  context 
of  finite  element  methods,  it  was  observed  that  under  certain  conditions 
the  convergence  rate  at  specific  points  or  along  specific  lines  is  higher  than 
the  global  convergence  rate  of  the  approximations.  These  phenomena  may 
occur  naturally  or  “artificially"  from  postprocessing  procedures.  For  linear 
problems,  an;  extensive  review  on  superconvergehce  may  be  found  in  the 
paper  by  Krizek  and  Neittaanmalci  [1].  For  nonlinear  problems,  few  results 
are  available,  even  though  numerical, evidence  strongly  suggests  that  most 
of  the  results  for  the  linear  problems  may  be  carried  over  to  reasonably  well 
behaved  Nonlinear  problems.  (See  Ghow  and  Lazarbv  [2],  Wheeler  and 
Carey  [3],  and  Carey  et  al.  [4].) 

In  this  paper,  we  identify  a  class  of  semilinear  problems  for  which 
many  superconvergence  results  rernaih  valid.  In  particular,  by  restricting 
oufselvesjto  biquadratic  elements  for  simplicity,  we  establish  superconver¬ 
gence  of  approximate  derivative  components  alongJines  through  the  Gauss 
points  and  superconvergence  of  the  gradient  at  the  Gauss  points.  We  also 
ptoversuperconvergence  of  the.  boundary  flux  and  linear  functionals  ob¬ 
tained  from  certain  postprocessing  procedures. 

The  standard  notational  cohvenri’ons  for  Sobolev  and  Hilbert  spaces  are 
employed  (e.g.,  see  Showalter  [15]).  For  a  given  domain  ft,,  the  Sobolev 
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space  k  =  {0, ±1,±2,,..}, p;£  [1, ooj,  is  equipped  wlth  the  usual 

norm  |lil)i.p.n  and  seminorm  l-li.nn.  We  shall  omit  the  index;/?  whenip  =  2. 
and  write  //*  for  if'*' ^  and  //  for //d  (ft),  the  subspace  of //‘(0)= with  ele¬ 
ments  vanishing  on  the  boundary  of  ft.  Since  the  distinction;  should  be 
clear  from  the  context,  we  use  (•,•)  to  denote  both  the;£’(ftj.;innef;prpduct 
and  the //“'(ft)  x  //d(ft)  duality  pairing. 

For  the  present  superconvergence  analysis,  we  restrict  the  domain  ft 
throughout  to  be  a  subdomain  of  consisting  of  a  union  oFrectangjes  with 
sides  parallel  to  the  coordinate  axes  and;  let  5ft  denote  ;t he  boundary.  We 
specifically  consider  the  class  of  semilinear  bdundaryrvalue  pifobjems  de¬ 
fined  by 

U  .  =  m  m 

,  U  =  0,  -6h:5ft 

where  (i)  a',/3  are  uniformly  Lipschitz  continuous  in  ft;  (ii);;a,  /?;.are  uni¬ 
formly  bounded  above  and  bejow  by  some  positive  const3nts  ci  .ahd;Co,  rcr 
spectively;  (iii)  f(u)  is  a  mohotbnically  decreasing  Lipschitl  continuous 
function  of  «  £  R  such  that  for  all  Hiv  in  fi, 

[/(w)  -  /(v),w  -  v]->  y\\w- .v||E  (2): 

with  y  >  0;  and  (iv)  for  all  v,w  in  //,  there  exists  a  constant'5  >  0  such  that 
Colw  -  v|t  -  y\w  -  v|ld  ^  5|w  vlf.  (3) 

Using  the  Green-Gauss  identities,  the  corresponding-weak  statenient  of 
eqs.  (1)  and  (2)  becomes:  Find  m  £  //  such  that 

=  [/(«), v]  for  all  v  £  //,  (4) 

where  a{u,v)  =  f-  p^:^)dxdy.  Under  the  assumptions  (i)-(iv) 

following  eq.  (1),  it  is  easy  to  verify  that  the  energy  functiphaJ  associated; 
with  a(M,v)  is;  strictly  convex,  and.thus  the  existence  andiurdq'ieness  of  the 
weak  solution  to  the  variationarproblem  (4)  is  guaranteed: 

To  construct  the  finite  element  spaces,  we  first  partition ;ft  with  rectan¬ 
gular  elements  whose  sides  are  parallel  to  the  coofdinate?axes.  AVeiassunie 
that  this  discretization  r*  is  regular  [6],  with  the  parameter  /i  denoting  the 
maximum  diameter  of  the  elements  in  Th,  and  satisfies,  the  uhvefse  hy¬ 
pothesis  hmin  s  Ch  (throughout  this  paper,  C  is  taken  aS;a  generic  cohctant 
and  is  independent  of  /i).  Introducing  a  local  polynomiaKbasis  of  degree  ^ 
on  each  element  and  imposing  continuity  across  element  boupdaries,  we 
define  the  standard  C°  Lagrange  piecewise-polynomial  finite  element 
space  5*  C  //‘(ft).  Let  5n  denote  the  subspace  of  functions  m *5'' tha  van¬ 
ish  on  5ft.  The  finite  element  approximation  «/.  £  Sn  of  eq..(4)  is  obtained 
by  solving 

ai,Uh,Vh)  =  (/(«/,), va)  for  all  v/,  in  5o=. 


(5) 
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Using  the  idea  of  arirelliptic  projection  operator  P  on  W  defined-by 

a{ii  -  Pu,Vh)  =  0  for  alhv/,  in  So,  (6) 

Noor  and  Whiteman  [7]  showed  that  the  global  energy  error  estimates  for 
the.ifinlte'element  approximation  of  degree  k  defined  by  eq,=(5)  are  simijar 
to  those  obtained  from  standard  finite  element  theory  for  linear  problems, 
i.e., 

||H  -  u/,||i  <  C/I'‘||«||*+I .  (7) 


These  global  estimates  are  optimal.  However,  for  linear  problems,  it  has 
been  observed  in  practice  that  the  finite  element  solutions  may  achieve  a 
higher  order  of  accuracy  0{h'‘*')  in  the  gradient  at  certain  Gauss  points. 
For  linear  problems  iii  one  dimension  this  problem,  as  well  as  nodal  conver¬ 
gence  of  the  solution,  has  been  extensively  investigated.  In  two  dimensions 
this  property  for  the  derivative  was  shown  to  hold  by  Ogahesjan  and 
Ruhovec‘[8]  and>Andreev  and  Lazarov  [9]  for  linear  and  quadratic'triangu- 
iar  dements,  respectively,  and  by  Zlamal' [10,11]  and'Lesamt  and  Zlamal 
'il2i:for  quadrilateral  elements. 

In  numerical  studies  by^Ewing  and  WHeeler.[13]  arid  Ewing  et  a|.  [14],  it 
was  observed -that  the.  respective  component  derivatives  exhibit  supercori- 
vergence  behavior  along  Gaussian  lines.  MacKinnon  and  Carey  [15]  used  a 
Taylor-series  analysis  to  prove  this  result  and^  further  superconVergehce 
properties^  Superconvergence,  along  Gauss  lines  has  aisd  recently  been 
demonstrated  by  Ewing  et  al.  [16];fdr  the’mixedmiethdd.  It  is  of  interest. to 
seelf  tht..je  naturally  occurring  superconvergence  pheripmena  are  also  pres- 
ent  for  the  semilinear  probleriis  (1)!^(2).  GHdw  and  Lazarov  [2]  and-Wheeler. 
and  Carey  [3]  have  performed  related  studies  for  the  orie-difnensidnal  non¬ 
linear  problem. Here,  we  consider  the  two-dimensional  case  and  biquadratic 
elements  (k  -  2). 

Appropriate  postprocessing  techniques  can  also  be  introduced  to  yield 
superconvergent  results.  For  example,  after  the  finite;  element  approxima¬ 
tion  is  obtained,  an  approximate  flux  function  may  be  computed  from  this 
solution  using  the -integration  by  parts  formula  (e.g.,  see  Carey  et  al.  [4]). 
More  specifically,  recall’ that-the  true  normal  flux  q  is  defined  on- 90  by 


du  ^du 


where  (ni,/Z2)^  is  the  dutward  normal  on-50,  arid,  by  virtue  of  the  Green- 
Gauss  identities,  that  ^  satisfies  the  equation 


(q,v)  =  a(M;v)  -  (/(m),v)  for  alliv  E  //'(O), 


where 


{q,v)  =  qvds 
Jin 


Thus,  an  approximate  flux  q/,  may  be  qbtainecLby  solving 

=  a{iih,Vh)  -  {f{iih),Vh)  for  all  Vh  G  S^. 


(8) 


(9) 


T 
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Note  that:fpr  convenience  we  have  neglected  the  effect  of  numericar  inte¬ 
gration,  in.  eq.  (9).  The  approximate  boundary  flux  thus  computed  is  super- 
convergent.  This  result  has  been  proved  for  linear  triangles  by  Lazarov 
et  al.  [17]  and  will  be  shown  tochold  for  quadrilateral  elements  in  the  next 
section.  A  related  result  is  given  by  Douglas  et  al;  [18].  Averaging  deriva¬ 
tives  between  elements  can  also  produce  siiperconvergent  results  (see  [15]). 

Another  postprocessing  application  arises  in  the  evaluation  of  integrals 
such  as 

F  =  F{u)  =  P  uiftdx,  (10) 

Jn 

where  u  is  the  weak  solution  of.eq.  (4)  and  i/r  is  a  sufficiently  smooth  func¬ 
tion.  An  obvious  approximation  to  eq.  (10)  is  then 

Fh  =  F{uh)  =  i  Uhtj/dx.  -(11) 

Jn 

With  the  aid  of  a  negative  norm  estimate^  we  show  that  Fh  is:a  superconver- 
:gent  apprpximation>to  F.  (For  related  work,;see  [20].) 


2.  SUPERCONVERGENCE  At  GAUSS  POINTS  AND  ALONG 
GAUSSIAN  LINES 

As  most  of  the  calculations  are  carrieci  out  over  the  reference  element 
e  =  [-1,1]  [“l».i]>  let  us  first  considej  the' transformation  of  an  elenient 

e  to  e  and  the  corresppnding.change  of  the  bilinear  form  Suppose.e  is 
an  element  inT*  with  center  (xo,yo)-  Kthe  side  pafallel  to  theA;  axis  fs  of 
-length  /i<^and  the  side  parallel  to  the  y',axis  is  of  length  ke,  then  a  ppint 
(x,y)  in  e  is  relateditp  corresponding  point  (f,  tj)  in  e  via  thedinear  map 

X  =  ATo  +  y  =  y(^,v)  =  yo-+  (12) 


Any  function  v  =  v(x,y)  defined  on  e  may  then  be-transforined  as  a  fiinc- 
;tipn  V  =  ’v[x{^,T]),  y{^,ri)]  and  analyzed -on  e.  Due  tp  the  restrictions  im¬ 
posed  ondhe  discretization  offl,  the  Jacobian  is  constant  on  each  element 
with  ="^htke.  Mpfeover,  since  there  are  np  mixed  derivative  ternis-in 
the  original  differential  equation  (1),  we  see  that 


/  V  f  f  du  dv  ,  , 

a{u,v)  =  2,  :  aT"7-  + 

JJf  dx  dx  dyrdy 

^Ikc  h,  ff:  -du  dv  \ 

<  \hc  J  Je  d^  d^  ke  j  j«-  drj  dr]  } 


(13) 


Let  G  denote  the  set  of  all  images  of  Gauss  points  (±\^/3,  ±V5/3) 
throughnhe  mapping  (12)  to  the  rectangles  in  the  discretization.  Wevnow 
show  that,  under  appropriate'  regularity  assumptions  on  «,  superconver- 
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igence  in  the  finite  element  approximation  for  the  components  of  the  gradi¬ 
ent  occurs  along  coordinate  lines  through  these=  image  points. 

Lemmafi.  Let  u  G  then  for  aiiy  v;,  G  5* 

|a(M  -  u,,Vh)\  ^  Ch^u\i  +  |m|4)|va1i  ^  (14) 

where  «/ is  the  finite  element  interpolant  of  m  in  5o. 


Proof.  The  proof  follows  closely  that  of  Zlanial  [10].'  Note  that  we  have 
not  imposed  any  boundary  condition  on  v*  and  only  seminorms  appear  on 
the  right-hand  side. 

We  see  from  eq;  (13).that  we  may  estimate  aiu  -  Ui,vh)  term  by  term.  By 
evaluating  the  coefficient  a  at  the  centroid  and;talcing  this  as  constant  on 
an  element,  we  may  simplify  the  iiitegral  terms.  This  implies  that  on  setfing 
a°  =  a(0,O),  we  have  to  estimate 


.  a(M  -  u,)  dvh . 
“'-'If— = 


-  a®) — - zjd^d'n  +  a°Li(uh. 


(15) 


where 


Uu)  = 


d{u  -j  jii)  dVh 


d^dr]. 


Now  q  is  uniformly  =Lipschitz  contmuous;  so  |a  a®l  ;<  |a|i,oo,;  and 
-hence  we  may  bound  the  first  term  on-the  right-hand  side  of  eq.  (15)  by 
G|a|i,»,c|M  -  M/|i;«iv/i|i.e.  From  standard  approximation  theory,  we  have 
|w  -  M/|if<  ^  C|m|3,>';  applying  this  result  and  transforming  back  to  the 
original  element.e,  we  have  the  bound 


(a  -  a®) 


„„  a(«  -  «;)  dVH 


d^drj  <jCh'^\u\i,,\yh\u. 


(16) 


To  estimate  Li(m)  we  observe  that  for  v/,  fixed,  Lj  is  a  bounded  linear 
functional  on  u  E  H\e),  so  |L|(«)|  <  cj|M||4,e-|vA|i,(.  Moreover,  note  that.if  m 
were  a  cubic  polynomial  on  e,  then  since  dvh/d^h  a  linear  polynomial  of 
we  would  have  Li{u)  =«0.  (We  only  heed  to  check  the  cases  u  =  and 
H  =  77^).  Then,  according  to  the  Bramble-Hilbert  lemma  [6],  we  may  re¬ 
place  the  norm  on  «  by.  a  seminprm  to  yield 

|Li(m)[  <  C\u\A.;\h\u;  ^  Ch^\uU,\vh\u-  (17) 


Combining  estimates  (16)  and  (17)  in  eq.  (15), 


M 


,  a(M  e  M,)  dvf, 


a- 


<  CP(|m|3,<  +  |«|4,«)|Va|i.<. 


(18) 
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In  a  similar  fashion  we  may  obtain  the  same  bound  for  terms  involving  77 
derivatives.  Thus  we  have 

\a(U  -  M/,V,,)|  <  E  C/i^(|m|3,,  +  |t/|4.^)|V/,|u 
<  C/t^(|M|3  +  lMl4)iv/.|i 

as  stated.  I 

From  Lemma  i,  we  may  derive  various  global  superconvergence  results 
relating  the  finite  element  approximation  h/,,  the  interpoiant  /i/  of  u  in  5o, 
and  the  elliptic  projection  Pu  of  u  in  So- 

Lemma  2.  Let‘t/ £ //‘‘(H)  be  the  solution  of  eq.  (4).  Then 

||m/  -  Ma||i  <  C/i^)|m1|4-  (19) 

and 

||F«  :£  C/j%/1|4.  (20) 

Proof.  Frpm-'eqs.  (4)  and  (5), .we  nave  for  al.lfy/,  in  So 

|a(MA  m,,va)|  :£  \a{u  m,,v/,)|  +  |(/(ma)  -  f(u),Vh)\. 

Using  Lemma -i  and  the  Lipschitz  property  pf(^  we  get 

|a(MA  -  iihVh)\  <  C/t^(|M|3  +  |m14)|va|i  +  CIIh  -  Va||o||v/.||o 

<  C/ji’||t<l|i!|vA||i ,  (21) 

where  we  have  assumed  the  global  L^-error  estimate  ||h  -  ha|'1o  <  C/r^||//||3 
for  the  finite  elethent  approximation.  This  is  the  standard  fesultior  the  lin¬ 
ear  problem  and  also  holds  for  the  semilinear  problem,  as  seen  in  Section  4 
using  negative  norm  estimates. 

Finally,  noting  that  iiu  —  «/  £  So,  and  that  fl(-,  •)  is  //-elliptic,  we  set 
Vh  =  Uh  —  ui  to  obtain  the  estimate  in  eq.  (19). 

To  see  that  the  estimate  (20)  holds,  we  first  show  that  ||/’//  -  «/||i  ^ 
Gh\\u\i  +  |m|4).  Now  a(M  -  ui,Vh)  =  a{Pu  -  ui,vo)  igT  all  v/,  in  So,  so  by 
Lemma  1  we  haye  |a(/*M  -  H/,Vfr)|  <  Ch%\u\3  +  [M|4)|vA|i..Againf observing 
that  Pu  -  Ut  is  in  SS,  we  set  v*  =  Pu  -  «/  and  make  use  of  the  H-elliptic- 
ity  of  a{-,  •)  to  obtain  the  desired  bound.  The  estimate  (20)  then  follovvs  eas¬ 
ily  by  applyingthe  triangle  inequality.  I 

The  above  lemma  uiaicaies  that  both  the  finite  element  approximation, 
of  «  and  the  elliptic  project  ion  of  m  are  exceptionally  close  to  the  inter- 
polant  oft/.  As  we  shall  see  in.  the  next  theorems,  the  derivatives  of  the 
interpolant  are  exceptionally  close  to  those  of  the, solution  at  Gauss  points 
arid  along  Gauss  lines.  This  can  be  used  to  then  show  the  derivative  super¬ 
convergence  property  of  the  finite  element  approximation. 


SEMiLINEAR  PROBLEMS  251 


Theoremsl.  Let  w  G  then.the  arithmetic  mean 

^=77-2  |V«(xj)  -  Vh/,(x^)|3 
iVc  J=1 


(22) 


where  Na  denotes  the  total  number  of  Gauss  points  XgrE  G,  g  =  1, 
2,...,Ng,  is  bounded  by 


A-<M 


or,  equivalently, 

Ng  ■ 

MA)'(Xg)|  <  C/l>||4:. 

Proof.  (See  also  [10].)  For  Gauss  point  (x^),  we  have 


(23) 


(24) 


^(«  T  U,)  _  2_  d(u  -;»/)  {Q 
dx  'hi 

where  ==  .(±\^/3,  ±v5/3):  Now  Z.2(m)  =  {d/dO{u:~  m/)(^)  is  a  linear 
iunctionai  of  u  G;H\e)  mth  \Li{i^)\  S  |m  ■?  m/|i;o»,;  <  C||m||4;c**  In  particu¬ 
lar,  if  M  were  a  cubic  pdlynomiai-, -then  Liiu)  would  vanish  identically, 
which  implies  on  applying  the  Bramble-Hilbert  lemma  that 


and 


|L2(m)|  ^  C1m|4,>*  S  Ch?|«|4,o 


[— («  -  h,):(xj) 


(25) 

(26) 


A  similar  bound  follows  for  {d/dy)  (,u  -  ui)  (Xg),  arid  combining  these 
results 

|V(h  -h,)(x,)|  <a^|«|4,,  (27) 

at  any  Gauss  point  x^  £  G. 

Summiiig  over  all  the  Gauss  pojnts  and  taking  the  average, 

—  1  |V(«  -  «;)  (x^)|  <  Nd'Cfi^-  2  |m|4,, 

g»\  t 

<  C/r'l«|4Nc''' 

<  C/r^(rr|4, 

since  =  5:,//“/iek,d^£/77  <  Cli^Nc. 

Recall  that  for  v*  belonging  to  5^  a  finite  dimensional  space,  we  have 
the  inverse  inequality  \vh\o.'«,t  ^  Ch~^\vh\o,c-  Setting  v*  =  (d/dx)(,Ui  -  ut,) 
and  va  =  (d/dy){ui  —  m*)  respectively,  we  can  use  this  inequality  to  obtain 
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No  ^ 

2;|v(m/  -  ma)(Xj)|  <  2  |v(m/  - 

«=i-  t 

<  c/i-‘5  lv(M/  -  ma)1o.. 

-e 

<  CNg\\ui  -  MaIi  . 

Thus  applying  Lemma  2,  we  have 

1  Nc 

TT  .E  |V(«/  -  ma)  (x,)|  <  CA^|iM||4 .  (30) 

iVc  g»i 

The  desired  bound  on  A  is  then  obtained  by  writing  u  -  Uh  =  «  - 
ui  +  {u/  -  Uh)  in  eq.  (22),and  applying  eqs.  (28)  and  (30).  The  second  esti¬ 
mate  (24)  follows  e^ily  from  the;  inverse  hypothesis  on  the  triangulation..;| 

Note  that  for /*«iwe  have  a  similar  result, d.e., 

1  Na 

^  B|V(«  -  Pm)(x,)|  <  Ch^(|M|3  +  |«|4):.  (M). 

This  estimate  may'be  obtained’by  noting  that  eq.  (29)  holds  with  ma  re- 
placedfby /?«  and  that  a  bound  for  ||Pm  -  M/||ids  available. 

The;fesu|t  in  Theorem  1  maybe  strengthened  with  the  introduction  of-a 
discrete  L^-horm,.as,in  ref.  [ll]:.The  seminprm  |v|a  is  defined  by 


where  =  (±v5/3;  ±vf/3). 

This  semiriorm  is  equivalent'uniformly  on  5*  to  the  seminorm  hi  i 


c"*|va|i  <  1va|a  ^.;c|va|i  for.all  va  in  5*.  (33) 


Since  for  Vh  E  So,  the  Poincare  inequality  implies  that  (va|i  isia  norm  on  5^ 
we  seedhat  |va[a  is  also  a  norm  otv^So-  With  this  in  mind,  we  write  the  semi- 
norm:|ijA  as;||i|A. 


Theorem  2.  Let  «j£  then 


||h  -  i<A|ift  ^  C/l^j|«|i4 


(34) 


and 


|1«  -  Pu\]/:^  Ch^(\uU  +  |«l4)  . 


(35) 
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Proof.  We  first-estimate  ||m  -  U/||a.  From  the  propf-of  Theorem  1,  we 
have  |L2(m)|  <  C|M|4;t*,  and  similarly,  for  the  77  derivative  term,  so 


<  Ch^\u\^ . 

Next  observe  that  («/  -  Uh)  E  .So.  Applying,eq.  (33)  and  Lemma  2,  we 
obtain 

||u/  -  unWh  S  Ch%U\\t 
and 

||M;  -T  Pu\\h  <  Ch^{\u\^  +  |u|4).  (37) 

The  triangle  inequality  may  how  be  used  to  show  the -desired  estimates. 
This  is  a^strongef  fesult  thhn  Theorem  1  since /I  <  C||m  -  | 

Next  we  turn  to  investigate  derivative  supwconvergence  along  special 
lines.  More  specifically,  we  show  that,  as  in  the  linear  problem^[15]  and 
also  the  dinear  mixed  finite  element  method- [16],  superconvergence  of 
dUh/dx-and  aMA/ay  . occurs ^aipng  lines. that  pass  through- the  Gauss  points 
and  parallel  to  the  y  axis  and  x  axis,  respectively. 

We  introduce- the  following  seminpf ms  for  sufficiently  smooth  v: 

f  A  f‘  ao 

[vlt  =  S  2  \  —{sf^k)'^ds  , 

I  ^  k  K  K\  dr) 

where  the  second  sum  is  performed- over  all  the  Gauss;points  lying  on  a 
straight.line  parallel  to  one  of  the  coordinate  axes,  w*  are  the  weights,  and 
^k  are  the  Gauss  .points  on  a  one-dimensional  reference  element.  In  our 
case  iVk  =  1  and  =  ±(V3/3). 

Theorem  3.  Let «  £  //‘‘(ft),  then 

[«  -  ma],-  <  C/i^|1m||4,  /  =  1,2  (39) 

and 

[m  -  Pm]/ <  C/i^(|m|3  +  |m|4),  i  =  l,2-  (40) 

Proof.  Again  we  estimate  [m  -  m/]i  and  [u/ -  ma]i  and  then  apply 
the  triangle  inequality.  Let;(^A,j),  —1  <  5  <  1,  be  fixed  and  set  L3(m)  = 
(a/a^)(u  -  Ui)((k,s).  This  is  clearly  a  linear  functional. of  m  G  H*(e)  with 
|L3(m)|  ^  [m  -  M/|i.»,e  <  C||m||4,;.  If  u  is  a  quadratic  polynomial,  m/  =  u  and 
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L3(m)  ==0.  Ifrii  =  then  M/  =  ^and  Zi3(«)'=  (3^^  -  l)|(a.j)  =  0,  since 
are  roots  of  the  Legendre  polynomial  3^^  -  1.  If  tr  =  77^  ui  =  ■jj  and 
Liifi)  =  0.  Thus,  L3(m)  vanishes  identically  for  cubic  polynomials  «,  and 
once  again  ajjplying  the  Bfamble-Hilbert=;}emma,  we  get 

|L3(m)|  :S  e\u\,.i.  (41) 

Thus, 

[m  -  ui]i  =  js  2  I  ^[L3{M)]^d5| 
f  r'  l‘« 

<\2c\ju\hdsj  ^^2) 

^'Chy\,. 

Observe  thabfor  vA  G  5*,  [va]i  =  H3va/3^||o.  since  the  two-point  Gauss  rule 
is  exact  fonquadratic  polyriomials.  As  ur-^  Uh  G  5o,  we-Have  [«/  ^  uh%  = 
\\(d/dx)(ui  rr  :UA)]||o  ^  \\ui  -  Ma||i,  SO  applymg.Lemma  2  wc  get  [u/  -  Uh]i  ^ 
G/i^|u|l4  and -hence  [m  «*]!:£  CA^||M||4.In:a  similar  fashion  we  may  show 
[u  -  P«]i  ^^GA?(|u|3  +  \u\i). 

To  provie-the  estimate  in-.the  [r]:  sehiinprm,  we  fix  -1  <  j  2  T,  set 
La^u)  =  (d/d^)  (u  -r  u))  {s,  §)  and  proceed; as  before.  I- 

Remark;  From  Theorem  3  we  may  deduce  that  at  the  Gauss  points  along 
the  sides  of'the  elements,  the  tangential  derivatives  areisuperconvergent. 
Note  also  that  even  though  we  have  restricted  our  anaiysis;to  quadratic  ele¬ 
ments  in  two  space  dimensions,  these  results  may  readily  be  extended  to  el¬ 
ements  of  higher  degree  and'in  different  dimensions  (see  also  Zlamal  [19]). 

3.  APPROXIMATE  BOUNDARY  FLUX;  CALCULATIONS 

Let  us  now  consider  the  problem  of  estimating  the  error  of  the  approxi¬ 
mate  boundary  flux  procedure  (9).  Subtracting  eq.  (9)  from  eq.  (8)  we  have 
forallvftG.S'' 

{q  -  qh,VH)  =  q(u  -  Uh,Vh)  -  {/(«)  -  .  (43) 

Using  the  global  L’-error  estimate  ijn  -  u'aIo  s  C/i^||h|3  and  the  Lipschitz 
property  of/,  we  get 

!(/(«)  -  f{Uh),Vh)\  <-C/l^j|M||3||VA|lo.  (44) 

In  eq.  (43),  we  write  «(«  -  Uh,Vh)  =  a{ii  -  iii,Vh)  +  «(«/  -  iih,Vh)  and 
bound  a(t/  -  h;,va)  using  Lemma  1.  Then  applying  the  Cauchy-Schwarz 
inequality  to  a{H/  -  Uh.Vh)  and  invoking  Lemma  2, 
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Thus, 

K<?  -  <7A,Vft)  I  S  C/z^lluyvAlli .  (45) 

Let  5^denote  the  space  of  functions  obtained  by  restricting  lunciions  in 
5*  to  the^bOundary  aft.  This  is  precisely  the  function  space  to  vvhich  the 
approximate  flux  qh  belongs.  As  S»  is  a  subspace  of  L\diX),  for  each 
9  G  L\dCi)  we  may  define  the  projected  function  in  S*  by  the  relation 


{q  -  Rq,Vh)  =  0‘  for  all  vA  in  5*.  (46) 

Setting^  -  qh  =  q  -  Rq  +  Rq  -  qn  in  eq.  (45)  and  applying  eq.  (46)  with 

Vh  =  Rq  -  q';,,  we  obtain 

1|/?^  -  ^  C/i"|1m114|1va|1i  ,  (47) 

Where  v*  now  represents  the  extensiph  of  thejboundary  function  Rq  -  qh  to 
a  fuhctiomthat  is  defined  over  ft  and  is  in  5*. 

In  view  of  the  inverse  inequality  [18] 

|V,1,  <  G/r'«llv,|lo.an,  (48) 

we  may  combine  eqs.  (47)  and  (48) -to  obtain 

l|/?P~PAllo.afl.sC/t'^"'l«ll4.  (49) 

It  is  not  difficult  to  see  that  for  the>f lux  function  q, 

Iq  -  I?pio,ah  =  inf|vA  -  pljaan,  v*  G  5a 

<  Ch^%qhn.,n  (50) 

^  Ch^ul. 

Thus,  we  may  combine  these  results  to  obtain  the  following; 


Theorem  4.  For  u  G  H‘’(fi),  the  error  of  the  approximate  flux  computed 
from  (9)  may  be  estimated  by 

'  h  -  9Allo.an  S  C/i^'^||h1|4  .  (51) 


Remarks.  This  result  is  not  entirely  satisfactory  since  in  computations^ 
the  rate  of  convergence  oTthe  fluX  approximation  has  been  observed  to  be 
0(h^)  for  biquadratic  elements,  indicating  that  the  above  result  is  probably 
suboptimal.  Note,  however,  that  it  is  still  a  superconvergence  result  since 
the  global  rate  of  convergence  of  is  expected  to  be  0{h^'^)  only. 

4.  EVALy^ON  OF  INTEGRALS 

In  many  applications,  integrals  pf  the  form  F(m)  =  or  G(u)  = 

JnVu  •  '^dx  must  be  evaluated  wheje  u  is  ihe  weak  solution  of  eq.  (4) 
and  \\i  and  4^  are  sufficiently  smooth  functions.  To  estimate  the  error 
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|/^(w) r-  F(Mfc)|  it  is  convenient  to  have  at  our  disposal  negative  norm  esti¬ 
mates  fore  =  u  —  Uh. 

For  ah  integer  ^  >  0,  we  have  by  definition 

|Hl-,  =  sup^,  (52) 

Thus,.for:^  E  //*(fl)\{0};,  we  seek  to  establish.a  bound  for  (e,  <^).  Before  do¬ 
ing  so,  let  us  recall  that  if  u.is  the  weak  solution  of  eq.  (4)  and  the  finite 
elemenbapproximation  in  eq.  (5),  then  by  the  mean  value  theorem  we  have 
the  following  orthogonality  relation: 

a{u  ^  Uh,Vh)  -  f  {uH),Vh) 

=  a(u  -  UhiVh)=-r  i  f  [  —[tu  +  (t  -  /)«/,] dmu  — Ma),va 
W  Jo  du  j 

=  a{e,Vh)  -  {hix)e,vi,) 

=  0  for  all  v*  imSo . 


Let  E  //‘(ft),  be  the  datfrih  the  auxiliary  problem 

Lijf  -  g(x)ilf  =  <f>  iii^ft,  (54) 

t/f  =  0  oh  3ft,  (55) 

with  L'the  linear  elliptic  Operator  in  eq.'(l)  arid  g(x)  -  (3//3m)[«(x)],  m  be¬ 
ing  the  weak  solution  of  eq,.(4).  Then 

ieA)  -  [e,  Lilf  -  gix)il/] 

=  u(e,#-  (eg(x);iff). 

Let>A  be  an  arbitrary  element  in  5o.  Setting  v/,  =  i/f/,  in  the  orthogonality 
relation‘(53)  and  subtracting  from  eq.  (56), 

=  aie, iff  -  iftn)  -  ([|(x)  -  /i(x)]e,i^)  +  ih{x)e,il/h  -  <fr) 

(57) 

#  Glejil^lr  -  ^i,\i  +  |([g(x)  -  h{x)]e,ilf)\  +  \{h{x)e,^  -  'Ph)]- 


Now 

l5(-»^)  “  ^'C-^)!  =  f  ?(“)  -  +  (1  -  0«a]  dt 

•  Jo\dlt  du 

-  \  l^(“*)  (1  “  ^)^'=  I"  ~  "aI- 

JQ  OU 

So  underUhe  additional  assumption  that  df/du  is  uniformly  Lipschitz  con¬ 
tinuous,  wc  Have  |g(x)  -  /j(x)j  ^  C\u  —  Uh\  and,  thus,  by  the  imbedding 
theorem  for  dimension  h  <  3, 

-  /i(x)]e,^)|  <  Cf(e^.^)|  <  C|H|g.4Wo  <  C|le||?|[^io.  (59) 
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Also,  as  |/i(a:)|  ^  C  fpr  all  j:  due  to  the  uniform  Lipschitz  bound  on  df/du, 

\(hix)e,tl/—  il/h)\  ^  C\e\i\il/  -  (60) 

Next,  leti^A  be  the  finite  element  solution  oftthe  variational  problem 

aitl/h,Vh)  -  [8(x)>l/h,Vh]  —  ((f>,Vh)  :fpr  all  Vft  in  So 

associated  with  the  linear  auxiliary  problem  (54),  (55).  Assuming  that  u  is 
sufficiently  smooth  that  ij/  G  HoiCl),  we  have  the  error  bound 

(61) 

From  eqs.  (57)-(61), 

If  5*  consists  of  elements  of  degree  A:,  then  5  +  1  <  A:  and 
(e, <!>);<  +  CA^1|Klk%,Wo. 

Thus  for  0  ^  5  <  A:  -  1,  we  have 

y.,  +  ch^\\u\u.  (62) 

In  particulafi  the  L*-norm  estimate  is  0(/i**'),  in  accordance  with  the  corre¬ 
sponding  estimate  for  linear  problems.  We  summarize  the  result  as  follows: 

Theorem  5.  Let  m  be  the  solution  of  eq.  (4)  with  m  G  //**'(ft).-Let  5*  be 
the;space  of  finite  elements  pf  degree-A:.  If,  in  addition  to  the  properties 
(i)-=§(iv)  listed  for  the  boundary  value  problem  (1),  (2),  the  nonlinear  forcing 
term/is  such  that  cl//dH  is  uniformly  Lipschitz  continuous  in  R,  then  for 
0  ^  ^  A:  -  1,  we  have 

|le||_,  =  0(/i‘*'*').  I 

Returnjng  to  the  problem  of  estimating  the  error  in  approximating  the 
integrarf(M)‘_by  F(ma),  we  have: 

Corollary  1.  Let  i/r  E  f/’(^l)  and  let  S'"  be  a  finite  clement  space  of  de¬ 
gree  A:.  Let  F(m)  =  / niulfdx.  For  0  <  s  <  k  -  I,  v/c  have 

\F{iO-  Filin)] 

Proof.  |F(m)  -  F(ma)|  =  1(22  -  Uh,<li)]  ^  C|[i2  -  Applying 

Theorem  5,  we  have  the  desired  result.  I 

Remark.  The  approximation  F(m*)  is  supcrconvcrgcnt  whenever  5  S  I. 
Note  that  if  i^r  is  very  smooth,  then  s  =  k  -  l  and  we  have  the  maximum 
rate  |F(m)  -  F(ma)  <  Oih^). 
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Corollary  2.  Let  £  [//**' and:let  5*  be  a  finite  element  space  of 
degree  k.  Let  G(«)  =  /jj  Vh  •  ^^dx.  For  0  <  s  <  k  -  I, 

\G(u)  -  G(u,)  I  <  C/«?"^"'|l^jU.I!'4*.  •  (63) 


Proof,  .‘\ppiying  integration  by  parts, 

G(tO  -  G(ma)  =  f  V(«  -  tih)  -  yf,dx 
Ja 


=  I  (m  -  iik)yt-  hds  -  [  (h  -  Uh)y  -yfdx, 

-'in  Jn 


where  n  is  the  outward  unit  normal  on  dfl.  Thus, 

lG(i/)  —  G(m/,)|  ^  [[m  —  ■  hll,+i/2.«i 


(64) 


<  Cil«  -  UA|I-,:nWUi- 

The  last  inequality  is  obtained  using  the  imbedding.'//^*'(ft)‘-»-  H**''^(dCl). 
Now,  by  applying  Theorem  s,  we  obtain.the  estimate  (62).  I 


5.  CONCLUSIONS 

Under  appropriate  assumptions  on  the  solution  u  and  nonlinear  function 
f{u)  v/e  determine  error  estimates  involving  the  finite  element  approxima¬ 
tion  to  a  class  of  second-order  elliptic-^milinear  problems.  In  particular, 
the  respective  first  derivatives  are  shown  to  remain:Superconvergent  along 
Gauss  lines  in  a  discretization  of  rectangular  elements^  The  analysis  is  ex¬ 
tended  to  include  postproccKing  formulas  for  the  boundary  flux  and  evalu¬ 
ation  of  integrals.  Thus,  one  can  extend  the  superconvergence  theory  to 
this  class  of  semiiinear  problems  under  the  stated  assumptions. 

This  research  has  been  supported  ip-part  by  the  Texas  Advanced  Tech¬ 
nology  Program,  the  National  Science  Foundation,  and  the  Office  of 
Naval  Research. 
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1.  Introduictioni 

We  consider^boUndary-value  problemsifor  degenerate  dliptic  systems- of  the 
form 

(i.l.a)  a{u)  “  V  •  A(x,  Vu);+  [  fi(x,s,  VyU),*  uds  9  /  ,  x  €  Q  , 

Jri  '  '  '  ' 

(i.l.c)’.  .B(a;,.y,  VyU.)  •  ,.  y.^€  Tx 

Here  0,  is  a  domain  in  R"  and  for  each  value  of  the  macro- variable  x  €  0  is 
specified  a  domain  fl,  with  boundary  P*  for  the  micro- variable  y  G  ft* .  Each  of 
a,  6,  /I  is  a  maximal  monotone  graph.  These  graphs  are  not  necessarily  strictly 
increasing;  they  may  be  piecewisercphstant  or  multi-'s^lued.  The  elliptic  operators 
•in  (I.l.a)  and  (l.l.b)  are  nonlinear  in;  the  ^adient  of  degree  p  —  1  >  0  and 
g  -  1  >  0,  respectively,  with  i -j-  ^  >  i,  so  some  specific  degeneracy  is  also 
permitted  here.  Certain  first  order  spatial  derivatives  can  be  added  to  (I.l.a) 
and  (l.l.b)  with  no  difficulty,  ^dscorresponding  problems  with  constraints,  i.e.. 
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variational  inequalities,  can  be  treated  similarly.  A  particular  example  important 
for  applications  is  the  iinear  constraint 

{l.l.c')  U,{x.^y)  =^.u{x),i  y.€’r*  ,  i 

which  then.replaces  (l.'l.c).  The  system  (-1.1^  with /i(s)  =  is  cedled  a  re^- 

ulafized  micTo~struciufe  mddely&nd  (l.l.a),  (l;l.b')v(r.T.c')  is=the  corresponding 
matched  micr/o- structure, .model  iniwhich  (formally)  e  0. 

The  time-depeh^eht  form  of  such  a  system  arises  as  a  model  for  the  flow 
of  a  fluid' (liquid  or  gM)  through  ^  fractured  medium.  This  is  assumed  to  be  a 
structure  of  porous  and  permeable  blocks  or  cells  which  are  separated  from  each 
other  by  a  highly  developed  system  of  fissures.  The  majority,  of^fluid  tr^sport 
will.pccur.8Jong  flow,paths  through.,the  fissure  , system, iand.the  relative  volume  of 
the  Celhstnicture  is  much  larger  than  that  of  the  fissure  system.  There  is  assumed 
to  be  no  direct  flow  between  adjacent  cells,  . since  they  axe  individually  isolated 
by  the- fissures,  but  the  dynamics  of  .the  flux  exchanged:  between  each  cell  and 
its  surrounding  fissxires  is  a  major  aspect  of  themodel,  The  distributed  micros 
structure  models  that  we  develop  here  contain  explicitly  the  local  geometry  of 
the  cell  inatrix  at  each  point  of'the  fissure  system,  and  they  thereby  r^ect  more 
accurately  the  flux  exChjmge  on  the  imcro-scale  of  the  individual  cells  across  their 
intricate  interface.  In  such^a  context,  .(l.l.a)  prescribes  the'flow  on  the, global 
scale  of  the  fissure  system  and  (l.l.b)  pves  the  flow  on  the  microscale  of  the 
individual  cell  at  a  specific  point  x  in  the  fissure  system.  The  transfer  of  fluid 
between  the  cells  and  surrounding  mediuni  is  prescribed  by  (l.l.c)  or  (l.l.c').  A 
major  objective'is  to .accurately  model  this  fluid  exchange  between  the  cells  and  ' 
fissures. 

The  plan  of  this  paper  is  as  follows.  In  Section  2  we  shall  give  the  precise 
description  and  resolution  of  the  stationary  problem  in  a  variatipned  forniulation 
by  monotone  operators  from  Bamach  spaces  to  their  duals.  In  order  to  achieve 
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this  we  describe  first -the  relevant  Sobolev  spaces,  the  continuous- direct  sums  of 
these  spaces, .and  the  distributed  trace  and  constant  functionals  which  occur  in 
the  system.  The  operators  are  monotone  functions  or  multi-'ralued  subgradients 
and  serve  asimodels  for  nonlinear  elliptic  equatioris  in  divergence  form.  In  Sec¬ 
tion  3  we  develop  an  abstract  Green’s  theorem  to.  describe  the  resolution  of  the 
variational  form  as  the  sum  of  a  partial;differential  equation  :Md  a  complemen¬ 
tary  boundary, operator.  Then  sufficient  conditions  of  coercivity  type  axe  given 
to  assert  the  existence  of  generalized  solutions  of  the  variational  equations. 

Systems  of  the  fom  (1.1)  , were  developed  in  Rosen- (1952),  Rosen  and  Winshe 
(1950),  piesler  and  Wilhelm  (1953)  in  physical  chernistry  as  models  for  diffusion 
through  a  medium  with  a  prescribed  microstructure.  JSimilar  systerns  arose  in  soil 
science  from  Barker  (1985),  v^  Genuchten  arid  Dalton, (1986)  and  in  reservoir 
models  for  fractured  media  in  Douglas,  et  al.  (1987),  Hon:iung  ( 1988).  By  ho¬ 
mogenization  methods  .such  systems,  wre  obtained  ^  limits  of  exact  micro-scale 
models,  and  then  the  effective  coefficieiits  are  cornputed  explicitly  from  local 
material  properties  in  Vogt  (1982),  Hornung  and  Jager  (to  appear),  ArBogast, 
Douglas  and  Homuiig  (1990),  An  eidstence-uniquenesS  theory  for  these  lineeir 
probleihs  whicli  exploits  the  strong  parabolic  structure  of  the  system  .was  given 
in  ShoWalter  and  Walkihgton  (1991).  One  can  alternatively  eliminate  U  and  ob¬ 
tain  a  single  functional  differentiai^equatipn  for  u  in  the  simpler  space  L^(n),,but 
the  structure  of  the  equation  then  obstructs  the  optiinal  parabolic  type  results; 
see  Horhung:and  bhowalter  (1990);  Also  see  Friedman  and  Tzavaras  (1987)  for 
a  nonlinear  system  w'ith  reaction-diffusion  loceil  effects. 

2v  The  Variational  Eormulatidn 

We  shall  resolve  oi  r  systems  as  monotone  operator  equations.  Let  ft  be  a 
bounded  domain  in  with  smooth  Boundary,  F  =  5ft.  Let  1  <  p  <  oo  and 
denote  by  X'’(ft)  the  space  oLp'^  power-integrable  functions  on  ft,  by  X‘^(ft)  the 
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essentially  bounded;measurable  functions^  . and- the  duality  pairing  % 

J{1 

for  any  pair  of  conjugate  powers,  p,  +  ^  ^  denote  the  space  of 

infinitely  differentiable  functions  with  compact  support  in  ft.  W*”’P(ft)  is  the 
Banach  space  of  functions  in  L^(ft)  for  which  each  partial  derivative  up  to  order 
m  belongs  to  P’(ft),  and  is  the  closure  of  C^(ft)  inW"*'^(ft).  See 

Adams  (1975)  tar  information  on  these  Sobolev  spaces.  In  addition,  we  shall  be 
given  for  each  z  e  ft‘ a  boxinded' domain  ft*  which  lies  locally  on  one  side  of  its 
smdoth  boundary  f*.  ^t  1  <  q  <  oo  and  denote  by  7*  :  W^’*(ft*)  -+  2/*(F*) 
"the  iraic  inap  which  assigns  boundary  values.  Let  T*  be  the  range  of  7*;  this  is  a 
Banach  space  with  the  norm  induced  by  7*  . from  W^’*(ft*).  Since.r*  is  smooth, 
there  is  a  unit- outward  normal  i/*(s)  at  each  s  € T*.  Finally,  we  define  IVj’*(fti:) 
to  ‘^be  that  closed  subspace  consisting  of  those  <p  €  W^’*(ft*)  with  7*^  G  K, 
t/c.,  each  %(v?)  is  cqnsijmt  a.e.  onTi.  We  shall  denote  by  Vy  the  gradient  on 
W^’(ft*)  and%  V  thcii^adient  ph  W’^'*’(ft)i 

In  order  to  prescribe  a  iheasurable  family-  of  cells,  {ft* ,  z  €  St) ,  set  S  —  jl" , 
let  Q  C  ft  X  R"  be  a  giyen  measurable  set,  and  set  ft*  =  {y  G  R"  :  (z,y)  G  (?}. 
Each  ft*  IS tneasurable  in  R"  and  by  zero-extension  we  identify  If’(Q)  i’(ft  x 
R”)  and  each  I’(ft*)  ^  Ir*(R").  Thus  we  obtain 

mQ)  S  {C7  G  L’(ftfI’(R"))  €  L’(ft*) ,  a.e.  z  eft}  . 

We  shall  denote  the  duality  on  this- Banach  space  by 

/.li  U{x,  »)$(*,  y)dy^dx,  17  £  i»(e)  ,  $  e  i«'  (Q)  . 


The  state  space  {01  our  problems  will  be  the  product  r^(ft)  x  L}  (:Q). 

Note  that  W^>®(ftx)  is  continuously  imbedded  in  L’(ftr),  uniformly  for  z  G 
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SI..  It  follows-that  the  direct  sum 

Wg  =  {U  €  L^(Q}:.V(x)  e  W^’\a,)  ,  a.e.  a:  €  SI  , 

and  f  111/(1)11^1,, *dx  <  00} 

J(i 

is  a  Banach  space.  We  shall  use  a  wiety  of  such  spaces  which  can  be  constructed 
in>  tKisJihanner.  Moreover  we  shall’ i^sume  that  each  is  a  bounded  domsdh  in 
H”  which  lies  locally  on  one  side  of  its  boundary,  F*,  and  F*  is  a  C^-manifold 
of  dimension  n  —  1.  We  assume  the  ;trace  maps  7^  :  )  — ♦  L®(F*  );are  uni- 

forrnly  bounded.  Thus  for  each  t/  €  W,  it  follows  that  the  distributed  irace  %{U) 
defined  by  7(f/)(a;,s)  =;(7*(l7(i))(s),  s  €  Fi,  x  €  fl,  belongs  to  L«(ft,X?(Fx)). 
The  distributed^trace  7  maps  W,  onto  7^  =  L^{Sl,Tx)  L^{Sl,L^(Tx)). 

Next  consider  the  collection  ^  €  fi}  of  Sobolev  spaces  given 

above  and  denote  by  Wi  =  i*(fli  Wx‘*(ft*))ithe  corresponding  direct  sum.  Thus 
for  each-l/  €  Wi  it.follows  that  the  distributed  tr£tce  7(^)  belongs  to  L^{Sl)..  We 
define  to  be  the  subspace  of  those  U  €  Wi  for  which  yiU)  €  Wo  ’^Cfi). 
Since  7  :  Wi  ^  L^(Si)  is  continuous^  is  complete  with  themofm 

This  Banach  space  Wo  ’^(fi)  x  VV,  will  be  the  energy  space  for  the  regularized 
problem  (M)  ffid  will  be  the  energy  space  for  the  constrained  probleih 
in  which  (l.l.c)  is  replaced  by  the  Dirichlet  condition  (l.l.c').  Note  that  Wq 
is  identified  with  the  closed  subspace  {(7l7,  C/]  :  1/  €  Wp’^^  of  .W’o  ’^(fi)  X  W,. 
Finally,  we  shall  let  Wo  denote  the  kernel  of  7,  Wo  =  {U  €  Wj  :  =  0  in  7^}; 

We  have  defined  to  be  the  set  of  w  €  W^’®(fir)  for  which  7xU> 

is  a  constant  inultiple  of  1*,  the  constant  function  equeJ  to  one  on  F*.  Thus 
;Wj  ’^(fix,)  is  the  pre-image  by  7*  of  the  subspace  IR -lx  of  Tx.  We  specified  the 
subspace  Wi  similarly  ^  the  subspace  of  Wg  obtained  as  the  pre-image  by  7  of 
the  subspace  L^(Sl)  of  Tg.  To  be  precise,  we  denote  by  A  the  map  of  L^(Sll)  into 
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Tq  given  by  At>(x)  =  u(x)  •  1,.  a.e.  x  €  fi,  u  €  A  is  an  isomorphism  of 

L’(Q)  onto  a  closed  subspace  of  Tq.  The  dual  map  A'  t^ing  into  L'^' {€l)  is 
given  by 

A'p(t')  =  =  f  !7i(lx)  •  v(ar)dx  ,  g  e  ,  v  €  , 

Jn 

;So  lye  have  A'ji(x)  r=  5*.('lx))  €  ft.-  Moreover,  when  follows 

that 

:  5x(lx)'?=  I  9x{y)4y 
Jr, 

the  irite^U  of  the  indicated  boundary  functional.  Thus,  for  g  E  I»®*(Sl,Ir®*(r*)) 
C  Tg,  \‘g  e  L^'  iQ)  is  given  by 

(2.1).  .:A'^(a:)i=  f  gxiy)4y  ,  a:e..  x  g 

Jr, 

The  imbedding  A  of  L®(fl)  into  7j  imd  its  dual  map  A  wilUplay  ^  essential  role 
in  opr  system  below. 

We  construct  the  elliptic  differentid  operators.ih  divergehce-form  as  realiza¬ 
tions  of  monotone  operators  from  Banach  spaces  to  their  duals.  Assume  we  are 
given  A  :  Q  x  R"  -*  R"  such  that  for  some  1  <.p  <.po,/5i  E  L^' (Q),  go  E 
c  and  Co  >  0 

(2;2.a)  A(x, |*)  is  continuous  in  R”  and  measurabteintx,  and 

+ ji(i) ,  ■ 

■(2:2.b)  (A{i,f)-i{a:,7),|'-0>O, 

(i2x)  X(i,f)-f>C(,#-5(,(i) 

for  a.e.  X  €  ft  and  all  ^,17  €  R". 

Then  the  global  diffusion  operator  A  :  WQ’^(ft)  W”bP:(.fl).is  given  by 
Au(v)  =  f  A(x,  Vu(x))Vu(x)<ix  ,  u,  v  E  Wj’^(ft)  . 

Ja 
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Thus,  each  An  is  equivalent  to  its  restriction  to  C^;(=Q),  the  distribution 

Au  =  1^00 (ft)  =  '  -A(:,  Vu); 

which  specifies- the  vEdue  of' this  nonlmear  elliptic  divergence  operator. 

Assume  we  io-e  given  B  :  Q  x  IR"  — ♦  R”  such  that  for  soihe  1  <  g  <  oo, 
Til  €  ho  €  and  Co  >  0  ■ 

(2i3.a)  is  continuous  in  ^  €  R  wid  me^urable  ih  (.x,yj)  €  Q,  and 

d-  hi(x,y)  , 

(2-3.b)  (f  (i,  y,  f )  rr  .B(x,  y,  0,  f  ^ 

(2h3.c)  l(x,y„|‘)^f>.co|f|^^/io(x,y) 

for  a.e.  (x,y)  €  Q  and  all  y  €  R". 

Then  define  for  each  X' 6  Bi  : 

Biw(v)  =  I  B(x,y,  Vyto(y)jVyv(y)dy  ,  to,  u  €  . 

Ja, 

The  elliptic  differentijad  operator  pn  is, given  by  the  formal  part  the 

distribution 

BrW~  ^xto|poo(„-)  =  '  B{x,  •,  Vyto) 

rn  Wg'^mx)'-  Also,  we.shall  denote  by  B  :  Wj  — >  W'  the  distributed  Operator 
constructed  froHirthe  collection  {0*  :  x  e  ft};  by 

BUix)  =  BriU{x))  ,  a.er  x  €  ft  ,  U  €  >V,  , 

and  we  note  that  this  is  equivalent  tp 

BU{V):~  f  B:c(Uix))Vix)dx  ,  U,V  €  W,  . 

Ju 

The  exchange  term  in  our  system  w'ill  be  given  as  a  monotone  graph  -which 
is  a  subgradient  operator.  Thus,  assume  m  ::R-+  R’^  is  epnvex  and  bounded  by 

(2.4)  m(s)  <  CCIsj’  + 1)  j  5  €  R  , 
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hence^.  continuous.  Then  by 


M9)  =1  /  m{g(:^xjys))dsdx  ,  ^  €,X’(fi,T’(rr)) 

JnJi\ 

we  obt|un,the  convex,. continuous?m  :  L%^,  L^{yx)).  -^/R'^ .  Assume  ^  ^  ^ 

so  that  consider  the  linear  continuous  maps 

A:  Wo^’''(Q)-*L^(il,L^rx))  ,  7  :  W, .D’(Q,L«(rx))  .  ^ 

Then  the  compositeTunction 

Mlu,U]  =  m{-yU  Au) ,  u  €  Wi’^(fl)  ,-U  €  W,  , 

is  convex  and  continuous  on  Wo  ’-(fl)  x  Wg.  The  subgradients  are  directly  com¬ 
puted  by  standard  results  .Ekland  and  T<^^  (1976).  Specifically,  we  have 
p  €  9m(S')  if  and  only  if 

^(x,s).'€'9m(flr(xj  5))  „  a^e.  .s?€  F*  ,  a.e.  x  €  , 

and  we  have  [/,F]:  €  5M[u,l7j;  if  and.  only  if  /  —  =-A'(^)  in  W"bi»':(Jiy  and 
F  =  in  VVg  for  some  €  -  Au)i. 

The  following  r^ult  gives  sufficient  conditions  for  the  re^u/arizcd  problem  to 
be  well-posed. 

Theorem  1.  Assume  1  <p,q,  ^  ^  >  ^,  And  define  the  spaces  and  operators 

c 

A,  7  as-above.  Specifically,  tbe^sets  {fix  '  a  €  ft)  are  unifoTmly  bounded  with 
smooth  boundaries,  and  the  trace  maps  {%}  ate  uniformly  bounded’.  Let  the 
functighs  A,B,  and  m  satisfy  (2.2),  (2.3)  (2.4),  and  assume  in  addition  that 

(2.5)  m(s)*>  cojs}’  ,  s  €  R  . 

Then-for  each  p^r  /  €  (fl),  F  €  W'  there  exists  a  solution  of 

(2;6.a)  u  €  :  A(u)  -  A'(m)  =  /  in  W-»’'(ft) 

(2;6.b)  U€  yVg  :  B(U^;+7W=  F  in  Wj 

(2.6.c)  p  €  L^'  (ft,  F*'  (Fx))  :  p  €  dm(7U  -  Au)  . 
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(2.6') 


For  any  such  solution  n'e  have 

(2.7^  .7  /i(i,  5)  (f5  =  (F(x),l‘x )  ,  a!.e.  X  €  ft  , 

where  1*  denotes  the  constant  dhit  function  In 

Proof.  The  system  (2.6)  is  a  “pseudo-monotone  plus  subgradient”  operator 
equation  of  the;  form 

TVi’»'(ft)  xyVg  :  for  all  [v,V]  €  l^’^(ft);x  W, 
vlu(u)  -t-  BU(y)  -f-  aM[u,  I7]([u,  VJ)  3  /(v)  -f  F{V)  . 

It  remains;only  to  verify  a  coercivity  condition,  neunely, 

as  irw||»v».>(n)  +  IlKlIw,.  -♦  +00  . 

Choose  k  =  max{lj/hl. :  y  €  ftx,  ®  €  ft}  and  let  i/g  =  (i/J, . . , ,  i/^)  be  the  unit 
normal  oh  F*.  For  v  €  W*’?(ft*)  we  have  by  Gauss’  Theorem 

/  (kl*  +  yn?|vl* -^nv)'=  /  5n(ynk(y)|’)dy 

/ft,  ‘  /ft, 

’  =  /  *'?(«)«niT*v(3)k'd5  . 

/r. 

Holder’s  inequality  then  shpws 

lklllf(ft,)  ^  9^lkllx,»(n,)ll^nvIlL?(ft,)  1 

and  from  this  follows 

?Ml7x»Hi,,r.)  +  (2t)’(9  - 1)’’’  lia»vlll,(„,, 

by  Young’s  inequality.  From  here  we  obtain 

(2.9)  collVll’.,,,,  < :ll7l'IIL,(fl,x.(r,)):+  ,  V  6  W,  . 
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Thus  from  the  a-pnon  estimate 

Au{u)  ^  BUii:)  +  M'yU  -  Xu)  > 

(2-10)  co||Vt.||J,,„,-%o||p(n)  +  co|feUlli,,e)-|l«ollt.(«^ 

,  +  colht; -- ,  t, 

the  Ppincare-type  inequality  (2.9)  and  the  equivalence  of  ||Vullj[,F(n)  with  the 
norm  pn  .  we  can  obtzdn  the  coerciyity  condition  (2.8).  Specifically,  if 

(2.8)  isibounded  by  K,  then  (2.10)  is  bounded  above  by 

■K’(i|Mllvy».^(n)  +  |l^y^llL«(Q)  +  j|7^.IU*(p.L*(r.))) 

^  ■^(ll“llvy‘;''(n)  +  IP|^lU«{Q)  +  iW  ~  +  ll'^wlli»(n))  , 

and  thedast  term  is  dominated  by  the  first.  This  gives  an  explicit  bound  on  each 
of  these  terms  and^thence,jon  +'|j,i^l]w,* 

.Ein^y,  we  apply  (2,6.b)  to,  theTfimctipn  V  €  W,  pven  by  .K(x^;y)  =  ,v(x) 
for  some' V  €  X*(ft),  and  this  shows 

p(7v)'=.(F,y) 

since  BU{V)  =  0,  Md  thus 

f  X'ji{x)v{x)dx  =./i(Au)  =  /<(7t;)  =  /  (F(x),l)u(x)dx  . 

Ju  Jh 

The  identity  (2.7)-now  follows  from  (2.1). 

The  more  general  case  (1.1)  of  a  monotone  pointwise  perturbation  is  easily 
handled  likewise. 

Corollary  1.  (p  :  R  — ♦  R^  and  i  :  R  R"^  be  convex  and  continuous, 
with  !p(0)  =  $(0)  =  0,  and  assume 


(2.11) 


V>(i)  <  C(|5i!  +  1)  .  $(s)  <  C(|s|’ +  1)  , 


5  €  R  . 


For  each  p^r  f  G  there  exists  a  solution  of 

(2.12;aK  u  €  d  +  A(u)V- A'(/i)  =  /  in 

(2.12.b)  UeW,  :  b  +  B(U)  +  7'(p)  =  F  in 
(2.12.c)  p  G  drh^^U  —  \u)  in  (flfL^  fTg))  ,  and 

(2.124)  a  G  d<p(u)  in  p'(n)  i  b  G 

For  any  such  solutiqnfwe  have 

(2.13)  J  K^yy)dy  +  J  p(x,s)ds  =  (F{x),lx)  ,  a.e.  x  G  tl  . 

Jqs  Jr. 

Proof.  Tnis  follows  M  above  but  with  the  continuous  convex;  function 

^[u,U^,=  f  (p(u(x)),dx  -^  I  [ 

Ja  JnJh, 

+  miiU  -  Au> ,  [u,  U]  G  H^:^^(ft)  x  >V,  . 

The  subgradient  can.be  computed  tenhwiseibecause  the  three  terms  are  contin¬ 
uous  on  I-* (ft),  F^iQ),  and  L’(ft,L’(r*)),  r^pectively. 

ReThark.  The  lower,  bound  (2.5)  on  m(*)‘  may  be  deleted  in  Corolla^  1  if  such 
a  lower  estimate  is  known  to  hold  for  It  is  also  uimecessary  in  the  matched 
microstructufe  model;  see  below. 

3.  The  Green’s  Formula 

In  order  to  prescribe  the  boxmdmy  condition  (1.1  .c)  explicitly,  we  develop 
an  appropriate  Green’s  formula  for  the  operators  Bx-  Note  that  we  cmi  identify 
P'(P-z)  C  (ft*)  since  Wp’*(ft*;);is  dense  in. (ft*),  so  it  is  mezmingful  to 

define 

Dx  =  {wG  W*-»(ftx)  :  .B*u»  G  i’’(ft*)}  . 

This  is  the  domain  for  the  absirdci  Green’s  Theorem. 


Lemma  1.  There  is  a  unique  operator  5,  :  D,  — ♦  for  vhich  Brtv  =  BrW  + 

7j5xU>  for  ^1  u<  €  Dr-  Thai  is,  we  have 

(3.1)  Bxw{v)  =  +  (5xtt>,7*u)  »  t>i€  , 

for  every  w  £  Dx- 

Proof.  The  strict  jnorphism  7*  of  onto  T,  h^  a  dual  7^,  which  is  an 

isomorphism  of  T'  onto  the  annihilator  in  W^**(ftx)^  of  the  kernel  of 

7x.  For  each  w  €  Dx»  the  difference  BxU>  —  B*ti»  is  in  Wo'’'(n  x)"*",  SO  it  is  equal 
to  7i(9xtw)  for  a  unique  element  5xto  €  T,. 

Remark.  The  identity  (3.1)  is  a  generalized  decomposition  ot  Bx  into  a  partial 
differenti^. operator  on  fix  and  a  boundary  condition  on  F*.  If  F,  is  smooth, 
Vx  denotes  the  unit  outwardmormal  on-F*,  and  if  *B(x,  *,  V^ta)  £  [Wr^’’*(ftx)]’*» 
then  w  E  Dx  and  from  the  classical  Grera’s  Theorem  we  obtain 

B,u;(v-)-^'(Bxto,v)L(n.)  =  /  ■B{x,sy%w)ux{s)'iv{$)ds,  ,  u  €  F*’*(ftx)  • 

Jr. 

Thus,  dxW  =  B(i,  •,  Vyti;)  •  i/x  is  the  indicated  normal  derivative  inX*"'(Fx)  v’hen 
B(x,*,  Vyw)  is  as  smooth  as  above,  wd  so  we  can  regiurd  dxtu  in  general  as  an 
extension  of  this  nonlinear  differential  operator  on  the  boimdary. 

The  formal  part  of  B  ;  W,  — »  W'  is  the  operator  B  :  rVq  — ♦  Wq  given  by 
the  restriction  jB(l7)  ^  Since  Wo  is  dense  m-Ll{Q)  we  can" specify  the 

domain 

D  =  {U  €  W,  :  BiU)  E  l^(Q)} 

on  which*we  obtain  as  before  a  distributed  form  of  Green’s  theorem. 

Lemma  2.  There  is  a  unique  operator  d:  D  —*T^  such  that 

B{U)(V)  =  V)t(Q)  +idU,yV) ,  UeD,V€  W, . 


Theorem  2.  Let  the  Sobolev  spaces- anddrace  operators  be  given  as  above.  We 
sununcLrize  them  indhe  following  diagrams 


L^m 

LOffl^L^Tx)) 

U.  : 

U 

u 

U 

W^'^flx) 

Tx 

w,. 

u; 

A»tl 

u 

w^‘^(dx) 

- — ►  Rlx 

Wi 

x*(n) 

^  ‘ 

U 

u 

T 

Wi'^iflx) 

— .  {0} 

Wo 

—  {0} 

-  in  n^hich  71  is  the: restriction  of  7  to  Wi.  ^0  dense  in  L^(Clx), 

L^{Q),  respectively.  Let  operator  Bx,  x  €  (l,  and  B  be  given  and  de£he  their 
formal  parteBx ,  Bj  as  above.  Then  construct  the  domains  £>x  j/D  ancf  boundary 
operators  dx,  d  as. in  Lemma  1  tmd  L^ma  2,  respectively.  It  follows  that  for 
any  U..€  W,, 

(a)  BU(x):-  Bx(U(x))  in  for  a.e.  x  e  ft,  and  Ue  D  if  and  only  if 

U{x)  £  Dx  for  a.e.  x  €  il  and  x  h-*  BxU{x)  belongs  to  L*  (Q); 

(b)  for  each  U  €  B>, 

dUix)  =  dxpix))  in  T^fdra.e.  xGft 

and 

BU  =  BU  +  'ri{X’dU)  in  >V,'  , 
and  for  each  V  €  Wi  we  have, 

I  BxU{x){V{x))dx  =  f  BxUix)\qx)dydx+  f  (a,t?(x),lx)(7i  V)(x)dx  . 

Jfi  Jq  Jn 

Proof,  (a)  For  V  €  Wo  we  obtain  from  the  definitions  of  B  and  Bx,  respec¬ 
tively, 

f  BU{x)V{x)dx  =  /  BU{V)dx  =  f  5*C/(x)^(x)) :dx  =  /  BxUix)V{x)dx  , 
Jn  Jn  Jci  Jci 
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aii'dtso  the  first  equality  holds  since  ’’(Hx)').  The  characteriza¬ 

tion  of  D  is  immediate  now. 

(b)  For 'V  €  W,  we  obtain  from  the  definitions  of  7,  respectively,  aiid 

(a) 

f  dU{'rxV{x))  dx=  f  -dUi'rVydx^  f  {BU  -  BU){x)V{x)dx 

Jq  Jq  J(i 

=  f  dx(Uix))'fj;V(x)dx  . 

JQ 

Since  the  range  of  7  is  =  L*  (fi,r^),  the  first  equality  follows.  The  second  is 
immediate  from  Lemma  2  since  on  Wi,  7  —  A  o  71  and  7'  =  7|  A',  and  the  third 
follows  from  the  preceding  remarks. 

Cbrpllary  2.  Ip  the  situation  of  Corollary  1,  f  E.L^'  fff) 

only  if  Au  €  {0).  and  L9'{Q),  saiddn  that  case  the  solution  satisfies 

almost  everywhere 

a(x)  €  5v?(u(i))  ,  a(x)  +  Xu(x) -t-  K^^y)dy  f(x)  +  /  F{x,y)dy  ,  x  €  fi , 

JQ,  JQ, 

u(s)  =  0,  seT, 

6(x,  y)  €  d^{U{x,y))  ,  6(x,y)  -f  BU{x,y)  =  F(x,y)  ,  y  €  ft*  , 

p{x\s)  €  dm(y[/(x,s)  -  u(x))  ,  d*(i7(x))(s)  -f  p(x,s)  =  0  ,  s  € T*  . 


Finally,  we  note  that  corresponding  results  for  the  matched  microstructure 
model  are  obtained  directly  by  specializing  the  system'(2.6')  to  the  space  Wq’^. 
This  is  identified  with  d7C/,  C/]  :  U  €  Wq’^}  as  a  subspace  of  PFo’^(ft)  x  W,, 
^d  we  need  only  to  restrict  the  isolution  [u,  U]  and  the  test  functions  [v,V], 
u  =  7V,  to  this  subspace  to  resolve  the  matched  model.  Then  the  exchange  term 
M  'does  not  occur  in  the  system;  see  the  proof  of  Proposition  1,  especially  for 
the  coercivity.  These  observations  yield  the  following  analogous  results  for  the 
matched  microstructure  model. 
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Theorem  1'..  Assume  lt<  p,  g,  ^  ^  and  define  the  spaces- and  operators 

as  before.  Let  the  functions  A,  B,  and  m  satisfyi(2.2),  {2.Z)  and  {2.4).  Then 
for  each- pair  /  €  F'€  W{  there  exists' a  unique-solution  of 


(3;2;a) 

(3.2.b) 

(3^2.c) 


u  €  Wl'^(n)  :  A(u)  =  /+  (F;1)  in 
U  eWi  :B{U^,:=  F  in 
jU  =  Xii  in  C%  . 


Corollary  1'.  Supposeip^i  are  given  as  before  and. assume  {2.11).  For  f,F  as 
above  there  exists  a  unique  solution  of 

(3.3.a)  u  €  Wi'^(Q)':  a  +  (6, 1)  +  4(u)  =  /  +  (F,T)  in 
(3;3.b)  .17  €  Wi:  :-  6  +^(I7)'.=  F  in 

(3.3.C)  yU  =  Xu  in  L^(ii)cT, 

(3.34)  a^€  d(p(u)  in  L«'(Cl)  ,  b  €  d^U)  in  L^'(Q)  . 

fy  , addition,  f  €  L®*  (J2)..  and  F  €  lA'iQ)  if  and  only  if  Am  €  and.  B{U)  € 

L^'"{Q),  and  in  that  case  the  sdlutioii^satishes  almost' everywhere 

a(^x)  €  d<p(u(x)y,  a(x)  -f  A«(ir)  +  \M.x,y)dy  =  f(x)  +  f  F{x,y)dy  ,  a:  €  fi  , 

fn.  do. 


u(s)  =  0 

K^,y)  €  a$(C7(i,y))  ,  .Ka;,y)  +  BU{x,y)  =  F(x,y)  , 
U{x,s)  —  y,{x),. 


s  €  r , 


ye  fi* 


5  6  . 


Remarks.  For  the  very  special  case  of  p  =  g  >  2-a,nd  a(u)  =  u,  b{U)  =  U  in 
the  situation  of  Theorem  1  it  follows  from  Brezis  (-1972)  or  Lions  (1969)  that  the 
Cauchy-Dirichlet  problem  for  (1.1)  is  well-posed  in  the  space  1-^(0,  T;  x 
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:Wp)  with  appropriate initial  data  u(x,0),  r/(a:,i/,0)  and  source  fuirctions 
F(x,y,t).  A  similar  remark  holds  in  the  case  of  Theorem  1'  for  the  matched 
model  wun  _(l. lx').  These  restrictivc’assumptions  will' be  substantially  relaxed 
in  Showalter  and  Walkington  (to  appear). 

Furthermore,  variational  inequalities  may  be  resolved  for  pfqblems  corre^ 
sponding  to  either  the  regularized  or  the  matched  microstructure  model  by 
adding  the  indicator  furiction  of  a  convex  constraint  set  to  the  convex  function 
Thus  one  can  handle  such  problems  with  constfmhts  on  the  global  variable 
Uj.the  local  variables  l^^or  their  difference  Xu  —  iU  on  theinterface. 
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Abstract.  A  model  of  hbniiheair  diffusion  through  a  porous  ihediuih  is 
considered,  where  the  solute  is  adsorbed  through  the  boundaries  of -the 
individual  cells  in  the  prescribed  microstructure,  and  the  flow  within  each 
cell  is  governed'byia'corresponding  porous  medium  equa,tion.  The  resulting 
system  is  ^owri  tb-be  well-pbseddn  appropriate  spacej.  and  certaim 
properties  of  the  solution  are  obtained  for  special  cases. 

1.  Ihtroduction. 

We  begin  with  a  description  of  the  porous  media  system  to  be  studied;  Let  0,  be 
a  bounded  domain  in  which  lies  locally- on  one  side  of  its  boundary T.  Denote 

by  A  the  Laplace  operator  in  and*by^  the  unit  outw^d  normal  on  P.  SuppOse 
that  Q  is  a  given  measurable  subsetof  ft  and  set  ft*  =  {y  €  :  [a;,y]  6  Q}. 

For  each  pOint  cc  €  ft  we  assume  ft*  is  a  domain  in  R"^  with  boundary  P*,  and  these 
satisfy  ^unifofmly  in  x)  the  sanie  hypotheses  as  ft  and  P;  denote  by  Aj,  the  Laplace 
operator  in  the  micro- variable  y  €  ft*  and  by  "n*  the  unit  outward  normal  on  P*. 
Suppose  we  are  given  five  maximal  nionotone  graphs  in  R  denoted  by  a, 

We  shalhcOnsidef  the  singuleir-degenerate  parabolic  system 


(l.l.a) 

■^a{w)-  A{w)  +  J^  VyW■^n,ds  ^f, 

x  €  ft  , 

(l,l.b) 

Vw  •  n  +  B  0  , 

re  €  P  , 

|a(iy)  -  6._yW  9  F  , 

y  €  ft* 

1 

#  e>m- 

- 

VyW  -  na:  +  I3(W)--  7(a(tt>)) 

for  <  >  0  with  initial  conditions;  a(to(a;,0)),  b{W(x,y,0))  prescribed  in  addition  to 
the  sources  /  €  ij^(n'x‘(d,  6o))  and  F  €  x  (0, 66)),  where  Q  =  Hj-gn  Thus 
(1.1. a)  governs  the  flow  on  the  macro-scale  of  the  porous  medium  fi,  i.e.,  through 
the  fissures  in  this  global- domain,  and  (l.l.a)  prescribes  the  flow  on;the  micro-scale 
of  the  indiyidualicell  located:at  the  point  a:  €  -12.  The  ihtegreil-iterm  in  (l.l.a)  is 
the  total  flux^fiowing  across  T®  into  the  cell,  and  this  flux  is  determined  by  (1.1.;0). 

The  monotone  graphs  are  possibly  multi-valued  functions,  so  we  obtain  inclu¬ 
sions  instead  oflequations:  the  corresponding  equation  holds  for  jome  selection  out 
of  the  graph.  In-general,  w(a;).  is  the  density  of.mobil  solute  at  a;;€  ft,  W(a:,t/)  is 
the  density  of;  solute  at  y  efti  adsorbed  in  this  cell,  a(tt))  and‘  Q:(:W^)  represent  cor¬ 
responding  concentrations  ofthe;Solute  present,;^!!),)  is  the.flux  a-cross T  at  a  given 
density  ti),  and-^(W)  —  7(j^)  deterinines  the  tr£insport  of  fliix  acfoss  F*  for  given, 
density  W  on  the  inside  of  ft*  and  concentration  a(u;)  on  the  outside  of  this  cell  nt 
X.  These  last  two  graphs. correspond  to  the  adsorption  isotherms^of  the. media. 

The  problem  (1.1)  will  be  regarded  as  an  abstract  Cauchy  problem  for  the 
evolution  equation 

(1.2)  u'(t) -f  4(u(t))  3 /(<)  ,  f>  0, 

in  the  Banach  space  X  =  (ft)  xi^  {Q).  Recall  from  [3],  [2]  that  OxiAntegral  solution 
of  (1.2)  in  a  Banach  space  X  is  a  u  €  C'([0,oq),X)  such  that  u(i)  €  dom(4)  and 

flLw(<)  -  ^\\  <  |ll^^(5)—  x\\  +  (/(r)  -  y,u(r)  -  x)  dr 

for  each  y  €  A{x)  and  0  <  s  <  t.  The  pmnng  in  the  integral  is  the  semi-scalar- 
product 

(y,x)  =  sup{(y,x*)  :  x*  €  X*  ,  ||a;*||  =  ||x||  =  x*{x)}  . 

A  (possibly  multi-valued)  operator  4  on  X  is  called  accrziive  if 


II -  n2||  <  ||«i  -  U2  +  e(/i  -  /a) II 


for  all  £  >  0;  €  A,  [Baj  /2]  €  A.  If  also  I  A-  A  is  onto  X  then  A  is  called  '^\f 

m-acc'retive.  The  fundamental; result  of  (3]  is  that  for  f  €  i/^([0,6o),X)  and  Uq  € 
doih(A)  thereds  a  unique  integral  solution  u  of  (1.2)  with  u(0)  =  uq.  This  integral 
solution  for  the  appropriate  rh-accretive  operator  m  L^(Q,)  x  L^{Q)  will  be  the 
,  “generalized  solution”  of  (i.l.). 

> 

2.  The  Stationary  Problem. 

■We  begin  with  some  notation.  For  1  <  p  <  oo  we  denote  by  LP(f2)  the  usual 
Lebesgue  space,  and  the  conjugate  exponent  by  p',  so  ^  ^  =  L  is  the 

Banach  space  of  those  functions  whose  derivatives  up  to  order  m  ^belong  to 
C^~(Q,)  is  the  space  of  infinitely  differentiable  functions^with  compact  support  in  ft, 
and  Wo”’^(ft)  ds  the  closure  of  (7^(ft)  in  W™’?(ft).  See;[l]  for  information  oh- these 
Sobolev  spaces.  Specifically,  let:ro.:  W^’^(ft)‘^  TP(F)iandTx  :  -W^>P(ftj;)  — >  Ii*’(rx) 
denote  the  trace  maps  onto  boundary  values. 

An  essential  construction  -L  continuous  direct  sum  or  random  Banach  space 
denote  by  LP(Q,  LP(Qx)y.  We  identify  iiP(fti-£P(ft*))  S  iJ’(Q)  by  the  Fubini-Tonelli 
theorems  with  the  duality 

L{L  dy^dx  =  We  LP{Q)  ,  $  6  Lp'{Q)  . 

Since  W^’^(ftx)  is. continuously  imbedded  in. jbP(fti),  uniformly  for  x  €  ft,  we  can 
identify 

XP(ft,  W'-P(ftx))  =  {W  e  LP{Q) :  W{x)  e  W'’P(ftx.) 

a.e.  X  e  a  ,  and  /  \\Wix)r„,.,  dx  <  oo}  . 

Jsi 

Similarly  we  construct  £^(ft,X^(ri)).  We  shall  assume  the  trace  maps  Tx  :  W^’^(ftx)  — > 
L^(Tx)  3i.ve  uniformly  hounded  and  define  the  distributed  trace  r(W)  in  i^(ft,  (Fx)) 
for  each  W  €  by  r(W)(a;,s)  =  Tx(W(i))(s),  s  €  Fx,  x  €  ft.  Thus 

r  :  ii^(ft,  W'bi(ftj.))  — >  L^(ft,L^(rx))  is  continuous.  Finally,  the  (constant)  imbed¬ 
ding  A  :  L^(ft)  -*  L^{Q,L^(Tx))  given  by  (Au)(.t,s)  =  u{x),  s  €  Fx,  a;  €  ft,  and  its 
formal  dual  A  :  Z»^(ft,  Xr^(rx))  -♦  L^{Q)  given  by  X{U)(x)  =  U(x,s)ds,  x  e  ft, 
will  play  an  essential  role  in  the  porous  media  system. 
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The  linear  Neumann  ;problem 


(2.1)  —/lv}  =  f  in  ,  Vw  -f[—g  in  L^(T) 

'  ^ 

is  an  ingredient  of  each  of  the  bpundaryr value  problems  below.  By  a  solution  we 
mean  % 

(2.1')  w  €  :  [  Vw  Vfp  =  f  f^+  f  gM^)  ,  e  . 

Ja  Jsi  Jr 

By  setting.^?  —  1  it  follows  tha,t  a  necessary  condition-for  the  existence  of  a  solution 
is 

h+  {9=0. 

Ja  Jr 

It  is-%ellrknown  that  this  condition,  is  also  sufficient  for  existence  of  a  solution. 
Solutions- ^e  unique ‘up  to  an  additive  constant,  so  there  is  a  unique  solution  w 
with  w  =  0;  This.solution;satisfies 

lkl|(^  <  const.  ||/||L»(ft)  +  ll^llLi(r). .  / 

For  a.mjpcimal  monotone  graph  a  in.  R  we  say  €  cr(ta)  in  for  a  realr 

valued  function  w  on  if  ^  €  Lp(Q)  and  ^(x)  €  q’(to(x)):  for  a.e.  x  €  fl.  The 
following  result  plays  a  pivotal -roledn  the  following. 


Lemnia  1.  Let  1  <  p  <  00,  /  €  Z^(Q),  g  €  L^{T)  and  w  be  a  solution  of  (2.1). 

Let  cr  be  a  maximal  monotone  graph  in  H  and  0  €  cr(0).  If  ^  €  cr(w)  in  L^’(Q,)  and 
r]  E  cr{Tow)  in  L^'{T),  then 

fi-b  I  grj  >  0  .  „ 

■>  '  -  -  .  __  J^  Jr  <r 

For  a  Lipschitz  a  the  result  is  imniediate.  For  the  more  genersil  case  it  follows  by  ~1 
the  methods  of  [5];  see  [4]  for  deteuls. 

The  stationary  problem  corresponding  to  (1.1)  is  the  following:  for  /  €  L^(f2). 
and  jP  €  Ir^(fl,Ir^(f2*))  =  I-^(Q)  given,  find 

w  €  B^^’^(f2)  ,  u  €  a{w)  in  jb^(ft)  ,  v  €  6(rou;)  in  i^(r)  , 

W  e  (n,  lF^>^(n*))  ,  u  e  a(T'F)..in  L\Q)  ,  y  6  fiirwym  y  {a, LyiTrS) 
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such  that 

(2:2.a5;  u  -  A#  + A^7(A(u))^^ ==  /  in  L^(Q) 

(2.2.b)  V  +  Vw  -  71  =  0  in  L^(F) 

(2.2.0:)  f/~A^W  =  F  in  L^(Q): 

(2.±j3)  ■  V  +  ^yW  •  n*  =  7(%))  in  i\(r*))  . 

Proposition- 1.  Assume  7  is  continuous  and  ;j7(r)[  <  i£:|r|,  r  6  IR.  For  j  =  1,2, 
let  fj  g  and  Fj  6  L^{Q)  be  given,  and  suppose  Uj,Vj,Wj,  Uj,Vj,  Wj  are 

corresponding  solutions  of  (2:2).  Then  ' 

(2;3.a)-  /  |ui  ^  «2|  +  /  #1  -  ^2 1  +  /  IF*  |..  |'r(«i)---  <y(ti2)| 

’  '  Jsi  Jr  Jtt 

Ja  JnJr.i 

<2:3;b)  /  f: 

•(ft  -/pi  JnJ.r„  _  . 

I  *  - 

<//  R-J=2|+/:|Eit|7(u,)-7(t‘4)l, 

JflJSlr  Jn 

laid,  hence, 

(2;4)  /yui_U2l+  f  |t;i-U2|+  /  f  \Ul-U2\  <  [  \f2-f2\+  /  [  \Fi-Fo\^ 
Jn  Jr  JnJQ^  ‘  ./n 

Corollary  1.  The  solutions  satisfy 

(2.5.a)  f  (ui  -ti2)'’'+  / (vi  -02)^  +  /  |F*|(7(ui)  -7(u2))'^ 
jft  ./r  dP- 

<  f(/l  -/2)++  /  /  72 

dn  JnvFx 

{2.6.b)  /  /  {0-i-0-2)++  f  f  {V,-V,)+ 

Ja  Ja.  Jn  Jr. 

<  f  f  (fl-f2)++  /  ,-’.l(7(«l)-7Wf  , 

dn  j!n 


(2.6)  /(«1— +  /(ui  -  V2)'*' +  /  [ 

Jil  Jr  Jil  Jttr 

<  f  iF\-F^ . 

Ja  Ja  Ju, 

Proof.  For  the  two  solutions  subtract  the  corresponding  equations  to  obtain  (2.2) 
with  M  =  ux  —  U2uV  ='V\  —  ‘f}2i  etc.  Note  that  w  is  a^splution-of  a  linear  Neumann 
problem.  Let  p  be  ;the, graph  sgn,  t.e.,  a(r)  =  1  if  r  >  0,  p’(r)  =  — 1  if  r  <  0,  and 
<t(0)  =  [—1,1].  Choose  ^(x)  =  cro(u(x)  +  w(x)),  tj(x)  =  cro(v(s)  +  fow(;s));  note 
,  that  ^  €  ■<7'(u),  ^  €  <t(u>),  since  a  is  rnpnotone,  and.17  €  p(u)^  7/.  €  cr.(Tow)  since  b  is 
monotone.  Thus  Lemma  1  appli^,-and  we  obtain  •(2.3.a);  Similarly,  W'(x,  •)  is  a 
solution  of  a  Newmeim  problem  on  ftx  for  a.e.  x  G  n,.iSo4he  preceding  ^giiment 
'at  a.e..x  G  fl,  followed  by  an  integration  over  ft,  leads  to  (2.3.b).  Also  (2.4);follows 
by  adding  (2.3);  Finally,  the  argum^t  above  with  <r  chosen  to  be,  the  graph  sgn't" 
yields  Corollary  1. 

The  inequality  (2.4)  is  the  fundamental  estimate  in X- (ft)  x  L^(Q)  on  a  solution 
of  the  system  (2.2).  .It  shows  that  the  dynamics  of  (l.'l)icprresponds  to  a  seinigroup 
of  contractions  inithis. space.  Mpreover,:(2l3)  displays  explicitly  the  dependence  of 
the  components  o{.(2. 2)  on  the  excht^ge  or  coupling  terns  in  L^(ft,L^(Fx)). 

Cbrbliary  2.  Suppose  f  €  L^(ft)  and  F  €  L^(Q).  For  j  =  1^2,  let  Vj  € 
:L^(ft,  V-  (Fi))  be  given  and  suppose  uy,  vj^wj  are  corresponding  solutions  of  (2.2.a), 
(2.2;b).  Then 

(2.7.a)  /  |ui  -  «2|  +  /  |Fx|  l7(ui)  -  7(^2)!  <  /  [  \Vi  -Vzl  , 

Jn  Ja  Ja  Jr^ 

Similarly,  let  uj  he  given  for  j  =  1,2  and  suppose  Uj,Vj,Wj  are  corresponding 
solutions  of  (2.2.a),  (2.2.)9).  Then 

(2.7.b)  /  f  \Ui-  U2\  +  [  f:  \Vi-V2\<  f  iFx|  |7(uO-7(«2)|. 

Ja  Ja^  Ja  Jr^  Ja 

Theorem  1.  Assume  the  following  of  the  maximal  monotone  graphs  a,  b,  a,  fi,  7; 
Rg{a  +  6)  =  IR  and  5[doma]  9  0  ; 

Rg{a  +  /?)  =  R  and  j3[doma]  D  Rg  jo  a  \ 
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7  is  LipscKitz'on  Rga:  |7(j'i)  —  7(^2)!  <  ^1^*1  —  j  J’2  €  Rga.  Then  for  each 
/  €  L^(p,),  F  €  L^Q)  there  is  a  solution  of  (2.2). 

'  Proof.  Let  V  €  L^(^,  L^(r*r))  be  given.  From  [4]  it  follows  there  exists  a  solution 
u,t;,rw  of  (2!2.a)  with  V  replaced  by  V,  (2.2.b),  u  £id.{w)  in  and  u  €  b{Tow) 

.  in  L^(f).  Fix  x  €  and  consider  the  boundary- value  problem 


(2.8.a) 

iy(x)  £  ,  U[x)  €  a(Ty(x))  in 

(H, 

y(x)  £  I3[t^W[x))  in  Li(rx) 

-- 

(2.8.b) 

I7(x)  -  Aj,Ty(x)  =  F(x)  in  L^(flx)  , 

- 

(2.8.c) 

y(x)  +  Vyiy(x)-.  nx  =  7((^M)(a:)> 

- 

By  use  of  [r,-s]  £  <^with;^(r)  3  7(v\ti(®))  we  construct  ihe  graphs 
and.;3(<)  =  ${t  r)  r-  7(Ait(x,))sfor  whicl  '^.8)  is  equivalent  to 

,  Ue.oc{W)  ,  V€3{rJW) 

U  -  AyW  =  Ffx)  -r-  s  in 

V  -h  Vj,W  •  =  0  in  - 

for  the  functions  U  '=  U{x)  —  s,  V  =  V’(x)  — 7(Au(x)),  W  =  l^(x)  — jij.and  0  £  5(0), 
0  £ j0(O).  Thus  (2.8)'bas  a  solution  J/(x),y(x),W(x)  by  [4]  for  a.e.  x  €'^  and  from 
the  estimates  above  this  family  of  solutions  is  a  solution  [7,  V,  Ty  of  (2.2.ar),  [2.2.  fi). 
Set  V  =  T[  V );  this  defines  a  self-map  T  bn  (Q,  I?  (Fx)),  and  from  the. estimates 
{2i3)  it  follows  that  T  is  non-expansive.  Moreover,  if  yi,y2  £  L^(fi,L'(rx))  and 
the  corresponding  solutions  of  (2.2)  wdth  y  replaced  by  Vj  in  (2.2,a)  8Je  denoted 
by  Uj,vf,Wj,  Uj,Vj,Wj  for  j  =  1,2,  then  (2.7.a)  with  Vj  =  and  (2.7.b)  hold.  If 
M  =  sup{|ri| :  X  £  ft),  then  we  have  from  the  Lipschitz  condition  on  7 

ThuS'with  (2.7)  we  obtain 


/ 


so  2?  is  a  strict  contraction.  The  fixed- point  V  =  T(V)  in  L^(Q,  L^(Qx))  yields  the 
solution  of  (2.2). 

Define  a  (possibly  multivalued)  operator  A  in  y.  L^  (Q)  a^  follows:  [/,  F];€ 
A[xi,tJ\  if  there  exist  to  €  TV^>^(n),  W  €X*(ft,:^'^*^’^(^x))  ^nd  for  which 

in  L\9)  , 
in  i'(r)  , 
in  L\Q)  , 
in  i^(ft,i^(rx))  . 

Thus  (2.2)  is  equivalent  to  (I  +  A)(u,17]  3  From  (2.4);  ^d  Theorem  1  it 

follows  that  (/ +..4)“^'  is  a  contrzwjtion  on  j&^(f2)  x  L\Q).  By  rescaling  the  five 
graphs  it  follows  easily.'that  the  same  holdsTof  (J  +  e.4)"^  for  every  e  >  0,  so  .4  is 
mTdccrtiivt  on  2/^(fi).:x-X^(Q). 


—Aw  +  X(rf^Xu)  —  V)  =  /  ,  ti  &  d(to) 

V  +  Vto  *  ft  ==  0  ,  V  G  P{tow) 

-AyW  =  F,  Uea{W) 

V  +  Vy W  •  nx  =  7(Au)  ,  V  €  M^W) 
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SUMMARY 

A  new  superconvergent  projection  formula  for  determining  vorticity  boundary  data  in  the  stream- 
functipn^yorticity  method  is  constructed; 


INTRODUCTION 

The  stream-function-vorticity  formulation  is  a  standard  approach  for  numerical 

treatment  of  2D  viscous  flows.  In£thiS;prpccdufe.the  problem  reducesto  solution  of  a  coupled 
pair  of  partial  differential  equations  -  the  vorticity  transport  equation  and  stream-function 
equation.  These  equations:  can-be  discretized  -and  iterativdy  decoupled,  and  then  solved  fpr 
itcf\t|s  approximating  ^  and  w,  Atwell-known  difhculty  in  this.algorithmi is:  the  problem  of 
specifying  vorticity  boundary  data  ^  essentiar  data  for  the  vorticity-tfansjjort  equation.  It  is 
standard  practice  on  rectangular  finite-difTercnce  grids  to  use  one-sided  finite  differences  of  the 
stream-function  iterate  to  compute  an  approximation  to  the  velocity  and  thereby  the  bounaary 
vorticity.'  A  similar  procedure  can  be.used  in  finite-element  methods  but  docs  npt  fit  naturally 
jjn  thisTramework.  Here  we  present  an  alternative  approach  based  on  superconvergent  flux 
ideas  that:  applies  to  both  straight  and  curved  boundaries  and  can  be  used  for  cither  finite- 
element  or  finite-difference  computations. 

FORMULATION 

Recalj  that  in  the  stream-function-^vOrticiiy  method  the  stream  funciibh  -^satisfies  the  Poisson 
equation 

— w  in  ft  (I) 

Here  w  is  the  vorticity  determined  from  the  vorticity  transport  equation  (e.g.  for  steady  Stokes 
flbw_Au-/),  The  strcam-functipn-vorticiiy  equations  arc  frequently  iteratively  decoupled  in 
the  numerical  solution  scheme.  Thc  objcctivc  hcic  is  to  construct  a  proccudre  that  exploits 
supcrconvcrgcncc  ideas  to  develop  a. post-processing  formula  from  (I)  for  approximating  the 
vorticity  on  the  boundary.  This  can  then  be  used  as  data  for  the  vorticity  transport  equation. 
We  introduce  the  familiar  Grccn-rGauss  formula  for  the  Lap|acian  operator 

f  (-Am)i;  dA:  d>'  =  f  Yh  •  Vu  d.Y  d>' -  f  o^  ds  (2) 

Jn  Jn  Jan  dn 
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where  u,  v  are  arbitrary  admissible  functions.  In  particular,  let  us  select  u-yp  satisfying  (I)  so 
that 

f  ujydjcd^=  f  Vxp-Vv  dx  dy  -  f  u  ~  dy  (3) 

Jn  J  n  Jan  on 

Now  (3)  is  an  identity  satisfied  by  the  solution  (tp.u)  for  arbitrary  admissible  v.  In  previous 
superconvergence  studies,  a  similar  construction  has  been  developed  and  the  approximate 
solution  introduced  to  obtain  superconvergent  boundary  flux  approximations  (i.e,  for  dpjdn) 
(e.g.  see  Wheeler^,  Carey ^).  In  the  present  case  we  instead  use  the  known  boundary  data 
^ypjdn  =  u„  where  s  is  the  tangential  direction.  Then  for  known  \p  and  dipldn  in  (3)  we  have  a 
projection  formula  for  vorticity  in  Q.  Now  set  the  approximate  solution  yph  for  yp  on  the 
discretized  domain  Q*  with  v  =  </>,•,  the  piecewise-polynomial  Lagrange  basis  function  assc  :i^ted 
with  node  /  on  the  boundary,  to  get  the  approximate  projection  (for  cj**  w) 

\  dx  dy=  \  VyPH‘V<Pi  dx  dy-  \  u^i  ds  (4) 

J  fl»  •  Jo.  J  an. 

Now  as  /  traverses  the  boundary  nodes  the  integral  on  the  left  involves  only  the  strip  of  elements 
adjacent  to  the  boundary.  Furthermore,  since  <Pi{Xj,  yj)  =  6y,  using  a  Lobatto  (node  point) 
quadrature  on  the  left  simplifies  the  expression  to  yield  an  explicit  superconvergent  extraction 
formula  approximating  the  vorticity  at  boundary  node  /  within  quadrature  accuracy  as 

Qi<j3*  =  f  VyPk'  V<t>i  dx  dy-  [  u^i  ds  (5) 

Jo.  J  ao. 

where  Qi  corresponds  to  the  accumulated  quadrature  weight  at  node  /  from  the  adjacent 
elements. 

Remarks 

1.  For  a  rectilinear  boundary  and  bilinear  elements  the  extraction  formula  (5)  is  equivalent  to 
a  one-sided  second-order  difference  approximation.^  For  higher-degree  elements,  more 
general  boundary  shapes  and  irregular  grids,  (5)  is  still  applicable.  The  scheme  has  been 
applied  in  finite-element  calculations  for  viscous  flow  applications  with  straight  and  curved 
boundary  geometries  as  well  as  moving  surfaces  (Murray*). 

2.  The  scheme  (5)  can  be  applied  with  finite-difference  methods  by  formally  introducing  the 
nodal  interpolant  of  the  finite-difference  solution  as  \pi,. 
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